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ABSTRACT 


This  st\jiy  Is  directed  toward  defining  and  evalu¬ 
ating  prediction  techniques  for  the  influence  of  the  iono¬ 
sphere  on  pseudo-range  observations  from  a  satellite  navi¬ 
gation  system.  The  synchronous  satellite  navigation  system 
for  ;jhich  the  study  vas  intended  is  discussed  and  the 
limitations  Imposed  on  the  correction  technique  by  its 
characteristics  and  the  proposed  navigation  equipment  are 
described.  The  salient  features  of  an  operational  correction 
system  are  defined.  Alternatives  for  the  content,  format  and 
method  of  conveyance  of  the  information  from  vhidi  the  cor¬ 
rections  are  to  be  made  are  discussed.  Three  distinct  algor¬ 
ithms  that  can  be  utilized  to  predict  the  ionospheric  Induced 
time  (or  range)  error  are  defined  and  evaluated.  One  algor¬ 
ithm  is  capable  of  long-term  predictions  while  the  other  two 
are  dedicated  to  near  real-time  predictions.  Data  accumu¬ 
lated  over  several  years  at  several  locations  form  the  basis 
of  the  evaluaticai  of  the  three  algorithms. 
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1.  IHTRODUCnON 


Recent  studies  of  navigation  by  sieans  of  satellites 
have  emphasized  the  passive  approach  in  vhlch  there  is  no 
active  radio  oartlcipatlon  by  the  navigator.  The  Air  Force 
currently  has  vinder  development  a  synchronous  satellite  navi¬ 
gation  system  which  would  be  capable  of  simultaneous  and  near 
instantaneous  determination  of  a  navigator's  position  and 
velocity  f Miller.  197O],  To  navigate, the  user  wouM  utilize 
observations  from  four  geographically  distinct  satellites 
whose  locations  must  be  made  known  to  him.  As  described  in 
more  detail  in  Section  the  observations  on  which  the 
navigation  is  to  be  based  are  the  times  of  reception  of 
signals  that  are  radiated  from  each  sateUlte.  Since  nomi¬ 
nal  times  of  transmission  are  kncwn,  it  is  possible  to 
compute  the  tr^sit  time  of  each  signal  and  consequently  an 
estimate  of  range.  This  technique  is:  generally  referred  to 
as  pseudo-ranging  because  of  the  uncertainty,  in  the  navi¬ 
gator's  estimate  of  the  times  of  tran^ssion.  Errors  in 
making  the  measurements,  predicting  the  ephemerides  and  in 
estimating  the  velocity  of  propagation  of  the  satellite 
signal  at  each  point  alcmg  the  ray  path  will  contribute  to 
reduced  navigation  accuracy. 

For  accurate  navigation  the  influence  of  both  the 
neutral  and  ionized  atmosphere  must  be  accounted  for  in  the 
estimation  of  the  propagation  velocity  of  the  pseudo-ranging 
signal.  The  neutral  atmosphere  of  interest  is  kncwn  as  the 
troposphere  of  which  97  percent  is  below  an  altitude  of 
24  km.  The  vertical  transmission  of  a  signal  through  this 
region  would  Induce  an  apparent  Increase  in  the  range 
measurement  of  from  2,3  to  2,6  meters  if  the  vacuum  speed  of 
light  were  to  be  used  [Honfield.  1971]*  effect  is 

independent  of  the  frequency  of  the  signal  up  to  about 
15  GHz,  For  other  than  vertical  incidence  the  range  error 
introduced  by  the  troposphe.  ;  could  be  appreciably  larger. 
For  exanple,  at  a  zenith  angle  of  85  degrees  a  typical  value 
would  be  in  the  nel§hl>orhood  of  26  meters  [Honfield,  1972]. 
By  utilizing  surface  breather  data  it  is  possible  to  estimate 
the  tropospheric  induced  range  error  for  vertical  incidence 
with  an  rms  error  of  only  a  few  millimeters  [Hopfield.  1971]. 

The  electronic  conponent  of  the  upper  atmosphere 
is  known  as  the  ionosphere.  The  velocity  of  propagation 
for  electromagnetic  radiation  In  this  region  is  different 
than  for  px’opagatlon  in  a  vacuum.  For  propagation  at  the 
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group  velocity,  the  measured  range  would  be  Increased  if  the 
vacuum  velocity  were  used,  which  is  similar  to  the  influence 
of  the  troposphere.  Both  the  group  and  jhase  velocities  can 
be  adequately  approximated  for  propagation  in  an  ionized 
medium  as  a  function  of  the  electrcxi  density  and  transmission 
frequency.  For  propagation  at  a  frequency  at  VHP  or  higher 
the  error  in  range  can  be  approximated  as  being  proportional 
to  the  electron  density  Integrated  alcxig  the  propagation  path 
divided  by  the  square  of  the  transmitter  frequency  [Budden, 

1961].  During  periods  of  maximum  solar  activity,  the  vertical 
range  error  that  could  be  expected  at  lew  latitudes  should  be 
no  larger  than  I50  meters  at  1600  JHz  [Bums  and  Fremouw.  1970] . 
The  1600  MIz  is  used  as  a  reference  because  it  is  the  carrier 
frequency  of  Interest  in  the  succeeding  study.  At  incidence 
other  than  vertical,  the  ionospheric  range  error  can  be 
appreciably  hither. 


There  are  three  alternatives  by  which  the  influence 
of  the  ionosphere  on  the  navigation  accuracy  can  be  minimized. 
A  first  alternative  is  suggested  from  the  fact  that  the 
magnitude  of  the  apparent  range  (or  timing)  error  in  any 
observation  is  inversely  proportional  to  the  square  of  the 
transmitter  frequency.  In  principle  then,  it  is  possible 
to  select  a  frequency  hi^  enough  that  the  errors  will  be 
reduced  to  acceptable  levels.  A  seccxid  alternative  would  be 
to  utilize  the  freqviency  dependence  of  the  induced  error  by 
radiating  two  distinct  coherent  frequencies  from  each  sat¬ 
ellite.  From  the  difference  in  the  apparent  range  to  each 
satellite,  the  major  influence  of  the  ionosphere  could  be 
eliminated.  This  is  the  technique  used  in  the  operational 
Navy  '  avlgatlm  Satellite  System  where  only  one  satellite  at 
a  time  is  observed.  Ihifortunately,  both  of  these  alterna¬ 
tives  have  the  same  deficiency  althou^di  not  to  the  same 
degree,  i.e.,  an  Increase  in  the  ca^lexlty  of  the  equipment 
for  both  the  satellites  and  in  particular  the  navigators. 
Since  any  navigation  system  must  serve  a  broad  class  of 
users,  the  development  of  reliable  and  lew  cost  navigation 
equipment  is  a  most  important  ccxislderatlon.  If  a  single 
frequency,  e.g.,  8OOO  MIz,  were  to  be  used  without  any 
additional  attempt  to  account  for  the  ionosphere,  the  1^ 
meter  vertical  range  error  at  I6OO  Ifiz  referred  to  above 
would  be  reduced  to  6  meters.  However,  several  studies, 
discussed  by  Kteane  [I969],  have  shown  that  operating  at  a 
l^quency  above  L-band  (1000  to  SOOO  Wz)  is  inadvisable 
principally  because  of  power  limitations  and  'Uie  need  fat 
hifidi  antennae.  The  utilization  of  two  frequencies  from 
each  satellite  would  necessitate  a  receiver  capable  of 


receiving  eight  distinct  signals  simultaneously.  However, 
this  approach  becomes  practical  if  a  five-channel  receiver 
is  utilized  in  which  one  channel  is  time  shared  among  the 
different  satellites  for  the  sole  purpose  of  obtaining  the 
ionospheric  refraction  correction.  The  third  alternative 
is  to  correct  for  the  ionospheric  effect  numerically, based 
on  predictions  of  the  characteristics  of  the  ionosphere. 

For  precision  navigation  a  correction  will  have  to  be  made 
by  tne  navigator  for  the  influence  of  the  troposphere.  At 
the  same  time  it  would  also  be  possible  to  computationally 
make  a  correction  for  the  Influence  of  the  ionosphere. 
Unfortunately,  the  electron  density  distribution  of  the 
ionosphere  is  characterlred  by  great  variations  in  magnitude 
which  makes  the  task  of  effecting  a  correction  more  difficult 
than  for  the  troposphere. 

If  it  were  possible  to  predict  accurately  the 
ionospheric  characteristics  and  be  able  to  convey  this 
information  to  the  navigator,  the  third  alternative  would  be 
most  attractive.  It  would  result  in  less  complex  electronics 
for  both  the  navigator  and  the  satellites.  Consetiuently, 
this  paper  is  directed  toward  defining  and  then  evaluating 
prediction  algorithms  which  could  be  -utilized  to  estimate 
the  ionospheric  correction  for  a  high-altitude  pseudo¬ 
ranging  satellite  navigation  system. 
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2.  SUMM/^  OF  BESUUTS 


inae  study  undertook  the  development  of  con5)uta- 
tional  techniques  and  algorithms  that  couM  be  used  in  the 
world-wide  prediction  of  the  ionospheric  induced  time  (or 
range)  error  that  would  be  experienced  by  satellite  pseudo- 
reinging  signals.  Several  alternatives  are  discussed  for  the 
transfer  of  the  ionospheric  correction  to  the  navigators. 
These  include  tables  or  charts  based  on  long-term  predictions 
and  a  series  or  matrix  approach  for  the  dissemination  of  the 
real-time  corrections. 

Hiree  distinct  algorithms  that  predict  the  charac¬ 
teristics  of  the  ionosphere  that  are  needed  to  estimate  the 
time  (or  range)  corrections  have  been  defined  and  evaluated. 
The  first.  Algorithm  I,  provides  long-term  (on  the  order  of 
2  weeks  to  3  months)  predictions  while  the  second  (Algorithm 
n)  and  the  liilrd  (Algorithm  III)  are  capable  of  near  real- 
tiro  predictions. 

The  data  base  provided  for  the  study  consisted  of 
Faraday  rotation  observations  obtained  at  nine  sites  in 
North  America,  one. site  in  Puerto  Blco  and  one  site  in 
Hawaii.  The  time  Intervals  over  which  the  data  was  obtained 
were  during  196^,  a  period  of  low  solar  activity,  and  dxuring 
1968  and  1969,  a  period  of  modei'ate  to  hi^  solar  activity. 
The  ionospheric  Induced  tiro  delays  inferred  from  this  data 
were  used  in  the  real-tiro  predictions  and  as  the  basis  for 
the  evaluation. 

Summarizing  the  results  of  the  study  over  the  data 
base  for  both  daytime  and  ni^ttlro  conditions  indicates 
that  Algorithm  1  is  able  to  predict  in  an  rms  sense  about 
^  percent  of  the  range  error  while  Algorithms  II  and  III 
are  able  to  predict  about  72  percent.  The  maximum  error 
observed  in  each  of  the  algorithms  occurred  for  the  same 
datum  in  ubich  the  vertical  time  delay  Induced  by  the 
ionosphere  was  86  ns  (at  I6OO  VRz),  Algorithms  I,  H  and  III 
predicted  values  of  ^  and  43  ns  too  low.  Because  of  the 
somewhat  limited  extent  of  the  data  base,  restraint  should  be 
exercised  in  attempting  to  interpret  these  values  in  a  global 
sense. 


The  above  figures  represent  the  poredictlon  ability 
on  an  Individual  range  measurement.  However,  the  navigation 
errcac  is  sensitive  only  to  the  uncorrelated  errors  that 
appear  in  four  near-simultaneous  measurements  made  by  a 
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navigator  from  geometrically  distinct  satellites.  For  this 

reason,  the  above  figures  would  represent  a  worse  case 

situation  if,  when  utilized  in  an  error  analysis,  they  are 

assumed  to  be  uncorrelated.  .  i 

* 

V  - 

In  general,  the  results  indicate  an  ability  to 
make  predictions  of  the  ionosjAieric  range  error  which  woulji. 
be  useful  to  some  classes  of  navigators.  Further  Isqpro^ments 
in  the  algorithms  would  expand  the  user  base.  However, ‘at- 
the  present  state-of-the-art,  a  dual  frequency  system  would 
appear  to  be  advisable  for  users  requiring  hl^  precision 
navigation  fixes  at  all  times. 
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3. 


DESCRIPTION  OP  THE  HIC2!  AETITUDE  NAVIGATION 
SATELLITE  SYSTEM 


3.1  INTRODUCTION 

As  discussed  In  Section  3,  a  navic-ti^'n  sat  "iiite 
system  is  being  developed  by  the  Air  Force  -cu  ra:a 
Missile  System  Organization  (SAMSO)  to  enable  military 
users  to  obtain  fast  and  precise  positicxi  and  velocity  infor¬ 
mation.  This  high  altitude  satellite  system  is  the  candidate 
of  the  Air  Force  for  the  Defense  Navigation  Satellite  System 
■which  is  a  proposed  tri-service  satellite  ne'twork  to  meet  the 
navigatiai  requirements  of  a  broad  spectrum  of  military  users. 
These  users  range  from  high  performance  strike  and  reconnais¬ 
sance  aircraft  to  ships,  helicopters,  vehicles  and  foot- 
soldiers.  System  accuracy  is  projected  to  be  on  the  order 
of  tens  of  feet  in  position  and  fractions  of  a  foot  per 
second  in  velocity.  Details  of  the  proposed  Air  Force  system 
are  given  by  Deist,  of  Air  Force  [1970]  and  Miller  [1970]. 
These  unclassified  documents  are  by  nature  abbreviated 
and  do  not  totally  characterize  the  Air  Force  system.  In 
any  event,  the  design  of  such  a  complex  system  is  naturally 
in  a  continual  state  of  evolution.  Consequently,  it  is 
impossible  at  this  time  to  present  a  definitive  description 
of  what  the  final  design  may  be.  However,  for  conq>leteness 
and  to  facilitate  discussion  and  design  of  the  ionospheric 
correction  system, a  navigation  system  is  postulated  here 
which  is  a  composite  of  the  information  given  in  the  above 
referenced  documents.  It  should  be  carefully  noted  that  the 
postulated  navigation  system  should  not  be  interpreted  as 
the  final  Air  Force  candidate.  When  the  design  becomes 
frozen,  any  differences  with  the  postulated  design  should  be 
studied  with  respect  to  their  Influence  on  the  ionospheric 
correction  system. 


3.2  SATELLITE  CONSTELLATIONS 

I 

The  satellite  network  is  to  consist  of  constella¬ 
tions  of  spacecraft  In  which  each  constellation  would  appear 
to  an  observer  on  the  ground  to  be  either  In  a  rotating  ”X” 
configuration  of  five  spacecraft  or  a  rotating  config¬ 
uration  of  four  [Miller.  1970],  Both  of  these  geometrical 
constellations  can  be  achieved  with  a  center  satellite  in  a 
synchronous,  near-circular  and  equatorial  orbit  and  peri¬ 
pheral  satellites  in  Inclined  and  elliptic  orbits  with  the 
same  mean  altitude  as  the  central  spacecraft.  In  each  case. 
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the  three  or  fovur  peripheral  spacecraft  wouM  describe, 
relative  to  an  observer  on  the  earth,  essentially  a  circle 
about  the  centrally  located  satellite.  A  single  constellation 
would  provide  coverage  to  a  region  Just  under  one-third  of 
the  earth* s  surface.  The  coverage  contours  and  satellite 
ground  track  for  a  typical  "X”  configuration  are  given  in 
Figs.  3.1  aud  3*2  which  are  taken  from  Woodford  et  al.  [1969], 
Three  similar  constellations  of  the  ”X''  configuration  properly 
spaced  in  longitude  would  provide  continuous  coverage  anywhere 
north  of  45  degrees  south  latitude. 

3.3  NAVIGATION  SYSTEM  SUPPORT 

The  satellite  navigation  system,  apart  from  the 
user  equipment  and  spacecraft,  is  to  consist  of  from  three 
to  six  ground  stations  constrained  to  possessions  of  the 
United  States,  computational  centers  and  possibly  area  cali¬ 
bration  stations.  The  ground  stations  are  necessary  to 
obtain  tracking  data  so  that  the  ephemeris  of  each  satellite, 
which  is  required  by  the  navigator,  can  be  extrapolated  Into 
the  future  with  sufficient  accuracy.  For  orbit  determlnatirai, 
two-way  (closed  loop)  range  and  range-rate  data  is  to  be 
obtained  at  each  tracking  site  for  all  spacecraft  above  the 
horlzcm.  With  a  knowledge  of  the  positlcais  of  the  tracking 
stations,  a  model  of  the  geopotential  and  models  for  effects 
such  as  the  luni-solar  gravitational  effect  and  solar  radia¬ 
tion  forces,  the  spacecraft  ephemerides  can  be  extrapolated 
in  time.  The  area  calibration  stations  are  to  consist  of 
navigators  at  known  locations.  From  their  measurements  it 
would  be  possible  to  estimate  systematic  errors  in  their 
vicinity.  This  information  would  then  be  made  available  for 
use  by  nearby  navigators. 

3.4  NAVIGATION  OBSERVATIONS 

To  eliminate  active  radio  participation  of  the 
navigator,  each  satellite  is  to  transmit  an  identifiable 
range  code  modulation  on  an  L-band  carrier  (approximately 
1600  Mz) .  These  transmissions  will  also  contain  the 
ephemerides  of  the  spacecraft  and  any  additional  information 
that  may  be  needed  by  the  navigator  for  callbzatlon  or 
effecting  corrections  to  the  observations.  Die  measurement 
that  is  made  is  the  interval  between  the  time  of  reception 
of  the  range  code  signal  and  the  nominal  time  of  transmission 
as  given  by  the  navigator*  s  clock.  The  estimate  of  the  time 
of  reception  is  to  be  accomplished  by  means  of  a  conrelation 
detector. 
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LONGITUDE 


Fig.  3.1  COVERAGE  CONTOURS  AND  SATELLITE  GROUND  TRACK 
FORTHE  "X"  TYPE  CONFIGURATION 
(WOODFORD  ft  al.,  19681 
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Fig.  3.2  COVERAGE  FOR  FIFTEEN  DEPLOYED  SYNCHRONOUS  SATELLITES 
(WOODFORD  at  •!.,  1909] 


The  time  Interval  measurement  can  he  thou^t  of  as 
being  partitioned  into  three  components 


''ej  >  J  =  1  to  4 


(3.1) 


where 


T  .  s  p^/c  ,  is  the  vacuum  propagation  time  from 
the  predicted  position  of  the  jth  spacecraft 
to  the  navigator 

p,  =  range  from  the  predicted  position  of  the  jth 
^  spacecraft  to  the  navigator 

c  s  speed  of  ll^t  in  a  vacuum 


"d  = 


- 


navigator’s  clock  offset 

time  delay  arising  frm  various  error  sources 

the  measurements,  which  are  the  time  intervale 
from  the  navigator’s  clock  epoch 


The  navigator's  clock  offset,  t,  ,  is  the  difference  between 
the  epochs  of  his  clock  and  the^clocks  of  the  spacecraft  if 
the  latter  are  assumed  to  be  s^chronized  among  themselves. 

The  error  sources,  To  >  are  several  in  number,  such  as: 
lack  of  synchrcxiization  among  the  satellite  clocks,  time 
delays  in  the  transmitters  and  receiver,  errors  in  the  pre¬ 
dicted  exiiemerldes  and  the  ionospheric  and  tropospheric 
effects  on  the  propagation  velocity.  If  the  error  so\a:ces 
are  small  enou^  so  that  they  can  be  ignored,  Tej  »  0  ,  and 
the  offset  of  the  navigator’s  clock  is  unknown, then  four 
observations  are  necessary  to  utilize  Eq.  (3.1)  to  effect  a 
position  estimation.  With  the  foia*  observations,  it  is 
possible  to  eliminate  and  obtain  three  equations  in 
terms  of  either  range  differences  or  time  differences.  Con¬ 
sequently,  the  satellite  navigaticm  system  is  conceptually 
similar  to  ordinary  hyperbolic  radio  navigation  in  which  the 
role  of  the  fi:«d  network  of  transmitting  stations  is  replaced 
by  the  satellite  constellation.  Throuf^  normal  doppler  ex¬ 
traction  techniques,  it  is  also  possible  to  measure  the 
range-rate  between  the  user  and  each  spacecraft  from  which 
an  estimate  of  the  velocity  of  the  navigator  can  be  inferred 
[Miller.  19703. 
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3.^  HtlNCIPAL  OBSERmiON/VL  ERRORS 


Tvo  of  the  major  sources  of  navigation  error  are 
the  Inaccuracies  in  the  predicted  eihemerldes  and  the  influ¬ 
ence  of  the  medium  on  the  propagation  velocity.  The  accuracy 
to  which  the  satellite,  ephemerldes  must  he  estimated  is 
within  the  state-of-the-art  [Melton,  1971].  The  influences 
of  the  troposphere  and  the  ionosphere  have  been  discussed  in 
detail  in  Section  1.  However,  it  should  be  noted  that 
because  of  the  hyperbolic  nature  of  the  navigation  solution, 
all  of  the  errors  introduced  by  the  troposphere  and  ionosphere 
do  not  propagate  into  the  solution  for  the  navigator's  posi¬ 
tion.  This  point  is  discussed  in  more  detail  in  Section  9.2.9. 


4.  lONOSHIERIC  EE5ECTS 


1^.1  EP5ECT  OP  THE  lONOSHIERE  ON  SATELUTE  RANGING 

In  the  proposed  hlgh-altltude  navigation  system, 
the  satellite-to-navlgator  geonetric  range  is  to  be  inferred 
from  a  measurement  of  the  satellite-to-navigator  propagation 
time  of  an  identifiable  range  code  modulation  (hereafter, 
ranging  signal)  on  an  L-band  carrier  (at  l600  MIz). 

At  VHP  and  above,  the  time  of  propagation  of  the 
satellite  ranging  signal  between  satellite  and  navigator, 
assuming  group  velocity,  is  shown  in  Appendix  G  to  be  given 
to  first  order  (in  rationalized  mks  units)  as 

T(G,t)  =  <2^  J  N(r,t)  ds  (4.1) 

°  r  G(t) 

where 

T(G,t)  >■  time  of  propagation  of  the  ranging  signal 
between  satellite  and  navigator  along  the 
geonetric  path  G(t) 

p(t)  =  satellite-navigator  geometric  range  at 
time  t 

c  s  speed  of  light 

f  ■  carrier  frequency,  Hz 

G(t)  s  geometric  range  path  at  time  t 

N(r,t)  ■  electron  density  at  r  and  t 

r  •  position  vector  to  a  point  on  G(t)  from 
a  geocentric  coordinate  system 

ds  =  differential  element  of  arc  along  0(t) 

As  shown  by  Eq.  (4.1),  the  effect  of  the  ionosphere  on  the 
ranging  signal  is  to  increase  the  propagation  tine  of  the 
signal  over  its  free  space  value,  p(t)/c  .  If  this  iono¬ 
spheric  induced  prcqpagation  time  error  were  neglected,  the 
effect  would  be  to  increase  the  apparent  range  between  the 
satellite  and  navigator.  The  ionospheric  contribution  to 
the  propagation  tine  is  proportional  to  the  integrated 
electron  content  along  the  range.  Thus  from  a  knowledge  of 
the  electron  density  distribution,  N(r,t)  as  a  function  of 
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space  and  time,  for  given  positions  of  the  spacecraft  and 
navigator,  it  is  possible  to  determine  the  ionospheric  group 
time  delay. 

4.2  FARADAY  ROEATION  OF  SATELLITE  SIGMIS 


Hie  quantity  on  which  the  evaluation  of  the  algor¬ 
ithms  was  to  be  based  in  the  present  study  is  the  equivalent 
vertical  laiospherlc  group  time  delay  obtained  from  syn¬ 
chronous  altltvkde  satellites.  At  the  present  time  a  large 
enough  body  of  direct  measurements  of  the  time  delay  does 
not  exist.  Hcxfever,  the^  does  exist  a  significant  amount 
of  Faraday  rotation  data  from  both  lew  and  high  altitide 
satellites.  iVom  Faraday  rotation  data  it  is  possible  to 
tofer  equivalent  vertical  integrated  electron  content 
(hereafter  EVIC)  data.  In  this  study,  Faraday  rotation  data 
were  used  to  infer  the  EVIC  from  which  equivalent  vertical 
gi‘oup  time  delay  values  were  obtained. 


I 


The  term  Faraday  rotation  refers  to  the  phenonenon 
of  rotation  of  the  plane  of  polarization  of  a  linearly 
polarized  electromagnetic  wave  as  a  result  of  the  blre- 
frlngent  property  of  the  ionosphere  in  the  presence  of  the 
earth’s  magnetic  field. 

At  VRF  and  above,  the  amount  of  Faraday  rotation 
occurring  between  the  satellit«%  and  observer  is  given  to  first 
order  (in  rationalized  mks  units)  as 


n(G,t) 


P  N(r,t)H(r)cos  e(G,r)sec 
f^  ’^G(t) 


x(0,r)dh  (radians) 
(4.2) 


where 


n(G,t)  s  Faraday  rotation  along  G(t)  at  time  t 

H(r)  -  magnetic  field  strength  at  r 

6(G,r)  s  angle  between  G(t)  and  -the  magnetic  field 
vector  (3(?))  at  ? 

x(0,r)  *  vertical  angle  of  G(t)  at  r 

dh  differential  element  of  altitude 

In  order  to  obtain  EVIC  data  from  the  Faraday 
rotation  measurements,  Eq.  4.2  is  rewritten  as 
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where 


and 


then 


n(G,t)  =  K  M(G,t)  f  N(?,t)ah 
G(t) 


(i^.3) 


j*  N(r,t)  H(r)  cos  e(G,r)  sec  x(G,r)dh 

M(0,t)  .  ^ , - 

X  N(r,t)dh 
G(t) 


K 


=  2«97  X  IQ- 


2 


n(G.t) 
K  M(G,t) 


»  J  N(?,t)dh  s  I  (G,t) 
G(t) 


(4.4) 


Note  thatj  strictly  speaking,  I  (G,t)  is  the 
Integral  of  the  vertical  con5>onent  of  N(?,t)  Integrated 
along  the  geometric  path. 


The  equivalent  ionospheric  vertical  group  time 
delay,  Tv(G,t)  ,  is  then  obtained  from  I  (G,t5  using 
Eq.  (4.1j  as 

V(8,t)  .  I^(o,t)  =  4.51  X  10-«  (4.5) 

The  value  of  M(G,t)  is  a  function  of  the  "shape"  of  the 
electron  density  distribution  and  thus  will  e:diibit  both 
diurnal  and  geographical  varlatlcns.  In  order  to  simplify 
calculations,  it  is  common  practice  in  many  quarters  to 
apprcnclma^e  M(G,t)  by  a  constant  value  obtained  by  evalu¬ 
ating  H(r}  cos  9(G,3^)  sec  x(G,^)  at  an  altitude  of  300  to 
450  km  along  the  geometric  »th  [ Yeh  ^d  Gonzales,  i960] . 

The  approximate  value  ol  M(G,t)  obtained  in  this  manner, 
using  a  constant  allitude  for  all  locations  without  resorting 
to  the  explicit  Introductim  of  an  electron  density  distri¬ 
bution,  becmnes  Independent  of  time  but  will  vary  geographically 

Using  either  M(G,t)  or  its  approximation  which 
is  denoted  by  ,  where  h  represents  the  height  used, 
"experimental"  values  of  EVIC  can  be  obtained  from  Eq.  (4.4). 

The  use  of  the  approximation  ignores  the  diurnal  variations 
of  N(G,t)  and  the  EVIC  obtained  using  tL  can  differ  by  as 
much  as  a  factor  of  two  frwn  that  using  M(G,t).  Pigs.  4.1a 
and  4.1b  show  comparison  of  and  M(G,t)  for  the 

Faraday  rotation  observation  stations  at  Arecibo  and  Sagamore 
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Hill,  These  curves  are  representative  of  the  type  of  differ¬ 
ences  to  be  found  at  other  locations.  Figs.  h.2a  and  4,2b 
shew  the  com^ison  of  the  Ty(G,t)  inferred  using  the 
M350  and  M(G,t)  of  Figs.  U.la  and  4.1b  respectively. 

Because  of  the  strong  diurnal  variations  of  M(G,t)  as  well 
as  the  magnitude  of  the  differences  of  the  means,  it  was 
concluded  that  this  approximation  to  the  M- factor  was  not 
adequate  for  the  purposes  of  this  study.  Consequently,  the 
value  of  M(G,t)  was  computed  for  each  data  point  at  the 
time  of  the  data  point. 

It  is  important  to  note  that  because  of  the  lr|"'^ 
fall-off  of  the  gecanagnetic  field,  the  Faraday  rotation  is 
sensitive  to  the  electron  density  in  the  Icwer  ionosphere 
and  relatively  Insensitive  to  the  electron  density  in  the 
higher  ionosphere.  In  order  to  infer  EVIC  data  from  the 
Faraday  data,  as  discussed  above,  a  normalized  profile  must 
be  introduced.  Electron  density  distributions  which  are 
similar  in  the  Icwer  ionosphere  and  depart  widely  in  the 
upper  icaiosphere  would  yield  the  same  measure  of  Faraday 
rotation  but  different  EVIC  data,  A  measure  of  the  signi¬ 
ficance  of  this  is  given  by  the  fact  that  totally  ignoring 
the  upper  ionosphere  in  the  determination  of  EVIC  data  from 
Faraday  rotation  data  can  resixlt  in  errors  of  the  order  of 
20  percent. 

The  Faraday  rotation  data  used  in  this  study  were 
measurements  obtained  at  VHF  (ncsninally  137  geo¬ 

stationary  satellites.  These  data  were  scaled  to  the  frequency 
of  interest  (16OO  MHz)  by  the  l/f^  frequency  dependence 
shown  in  Eq.  (4.2).  Hewever,  it  ^ould  bb  noted  that  due  to 
the  possible  existence  of  sharp  gradients  of  electron  density 
in  the  iemosphere,  this  scaling  procedure  could  result  in 
differences  in  the  inferzed  EVIC  from  that  actually  observed 
at  L-band  [Outer,  1963]  and  [Kelso,  1964],  This  concern  is 
discussed  in  Section  10. 7  es  en  area  for  further  study. 
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Fig.  4.2a  DIURNAL  VARIATION  OF  VERTICAL  TIME  DELAY  BASED  ON  M  (C.t)  AND  M350 
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Fig.  4.2b  DIURNAL  VARIATION  OF  VERTICAL  TIME  DELAY  BASED  ON  M  (G.t)  AND  M3S0 


5.  HIOPOSED  OHKAOJIOML  IONOSPHERIC  CORRECTION  SYSTEM 


5.1  INTRODUCTION 

This  section  is  aevoted  to  delineating  an  opera¬ 
tional  computational  system  by  vhich  the  ionospheric  cor¬ 
rections  can  be  estimated  and  forwarded  to  the  navigators. 

The  computatlcxis  are  discussed  apart  from  the  prediction 
techniques  that  must  be  utilized  to  estimate  the  electron 
density  distribution  of  the  Icxiosphere.  Discussion  of  the 
details  of  the  prediction  techniques  are  given  in  Section  6. 

5.2  CONSTRAIMPS  ON  THE  OPERATIONAL  IONOSPHERIC  CORPIECTION 
SYSTEM 

Biere  are  several  Important  factors  and  constraints 
that  must  be  reflected  in  the  synthesis  of  such  a  system. 
!niese  are: 

5.2.1  Geometrical  Characteristics  of  the  Satellite 
Constellation 

The  geometrical  characteristics . of  the  candidate 
constellations  have  been  discussed  in  Section  3>  Ground 
tracks  and  coverage  contoiurs  for  the  "X”  type  configuration 
are  displayed  in  Figs.  3.1  and  3.2.  In  both  the  "X"  and 
”Y"  constellations,  the  peripheral  spacecraft  rotate  about 
the  geostationary  central  spacecraft  such  that  their  naninal 
positions  repeat  themselves  every  2k  hours.  Departures  from 
the  nominal  orbit  arise  principally  from  the  luni-solar 
gravitatlOTal  effects  and  the  influence  of  solar  radiation 
pressure.  By  periodically  adjusting  the  orbits,  the  magni¬ 
tude  of  the  departxures  could  be  minimized.  This  attribute 
of  a  repeating  ground  track  correlated  with  the  solar  day 
is  emphasized  because  it -is  demcaistrated  in  Section  ^,k  how 
this  can  be  utilized  to  advantage  in  implementing  the 
ionospheric  correction. 

5.2.2  Spacecraft  Hardware  Limitations 

A  description  of  the  spacecraft  hardware  limita¬ 
tions  as  regards  the  ionospheric  correction  is  given  by  Dept. 
of  Air  Force  [1970].  Real-time  transfer  of  information  from 
which  the  correction  can  be  made  is  to  be  accomplished  by 
utilizing  a  portion  of  the  message  content  of  the  L-baxd 
ranging  signal.  This  link  is  also  utilized  to  transfer 
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spacecraft  ephenerides  aM  ranging  signal  tine  synchroniza> 
tlon  data.  For  the  purposes  of  this  study,  a  maximum  contin¬ 
uous  data  rate  of  10  bits  per  second  per  spacecraft  was 
allowed  for  the  ionospheric  correction. 

5.2.3  User  Hardware  Limitations 

Because  of  the  large  number  of  potential  \uters  of 
the  navigation  system,  a  significant  portion  of  the  overall 
cost  will  be  the  navigation  equipnent.  Consequently,  there 
is  a  strong  desire  to  minimize  the  cost  of  the  navigation 
equipment.  The  navigator  is  envisioned  to  have  a  digital 
cooqniter  in  order  to  compute  his  position  from  the  satellite 
observations.  This  device  can  also  be  used  to  cooqpute  the 
ionospheric  (as  well  as  the  tropospheric)  corrections.  The 
significant  limitation  on  the  computer  is  the  maximum  memory 
that  can  be  allocated  for  the  ionospheric  correction  which 
was  tentatively  limited  for  this  study  to  2000  24-bit  words 
[Dept,  of  Air  Force.  19701* 

5.2.4  User  Accuracy  and  Operational  Begulrements 

The  user  population  constitutes  a  broad  spectrum 
of  accuracy  requirements  and  operational  characteristics 
which  InqActs  significantly  in  the  design  of  the  operatl,onal 
correction  scheme.  Because  of  this,  several  levels  of  iono¬ 
spheric  correction  accuracy  should  be  made  available,  each 
being  suited  to  the  particular  needs  of  a  class  of  navigators. 
There  will  be  those  whose  accuracy  requirements  will  be  such 
that  they  will  not  require  any  correction.  Because  of  this, 
the  correction  system  must  be  implemented  in  a  manner  in 
which  it  is  transpaxent  to  these  users.  In  some  instances 
a  user  may  have  svich  lax  operational  requirements  that  he 
may  be  able  to  schedule  his  observations  during  the  evening 
hours.  The  influence  of  the  ionosphere  is  usually  signifi¬ 
cantly  less  in  the  evening  hows  than  during  dayli^^t.  In 
this  case,  the  navigator  msy  need  only  very  approximate 
estimates  of  the  correction,  if  any  at  all.  There  will  also 
be  users  whose  Interest  will  be  in  position  determination 
as  opposed  to  navigation.  For  these,  it  would  be  possible 
to  effect  an  ionospheric  correction  well  after  the  navigation 
observations  are  obtained.  Then,  these  corrections  could  be 
estimated  based  on  ionospheric  data  accxxnulated  before, 
during,  and  after  the  observation  period. 
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The  bulk  of  the  users  will  he  those  Interested  in 
navigation  in  which  it  is  necessary  to  effect  corrections  in 
real  tine.  For  discussion  purposes  it  is  convenient  to 
arbitrarily  categorize  these  users  into  three  groups  of 
distinct  accuracy  requirements  characterized  by  their  desire 
to  utilize;  (a)  long-term  predictions  generated  3  to  6 
months  in  advance,  (b)  short-term  jaredlctions  generated  1  to 
l4  days  in  advance  and  (c)  real  time  estimates.  Whereas  high 
performance  aircraft  mi^t  require  the  best  estimate  possible 
(i.e.,  a  real-time  prediction),  a  merchant  ship  may  be  satis¬ 
fied  with  a  less  accurate  estimate  (i.e.,  long-term  prediction) 
if  the  cost  is  appreciably  less. 

5.2.5  Nature  of  the  Influence  of  the  lonosiihere  on  Navigation 

©le  nature  of  the  influence  of  the  ionosphere 
on  navigation  is  important  in  the  de'/elopnent  of  an  opera¬ 
tional  correction  system.  The  effect  on  the  group  time 
delay  can  be  conveniently  discussed  by  utilizing  Eq,  (3.I) 
which  represents  the  measurement  made  by  a  navigator  as 
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J  =  1  to  4 


To  eliminate  the  effect  of  the  navigator’s  clock  offset  it 
is  necessary  to  difference  the  members  of  a  set  of  four 
observations  obtaining  three  equations  of  the  form 

(Tj-T^)  =  (Tpj-Tpj)  +  (Tgj-Tgi)  ,  .J=l  to  4  and  i=l,2,3  or  4  (5.I) 

In  this  context,  it  is  clear  that  it  is  the  diffeirences  of 
the  errors  which  propagate  into  the  navigation  solution,  not 
the  individual  errors.  If  all  the  t  (k=i,j)  had  the 
same  magnitude,  the  navigation  solution  would  be  unaffected. 

In  general,  this  will  not  be  the  case.  As  discussed  in 
Section  4,  the  ionospheric  Induced  error  is  proportional  to 
the  integrated  electrcm  content  along  the  path  fi’om  the 
navigator  to  the  spacecraft.  The  differences  in  the  iono¬ 
spheric  errors  will  be  a  function  of  the  horizontal  gradient 
of  the  electron  density  distribution  and  the  differences  in  the 
elevation  angles  of  the  spacecraft  with  respect  to  the  navigator. 

5.3  HAVIGAT(»*S  COMPOTATION  OF  lONOSIHERIC  CORRECTION 

!lhere  are  several  alternatives  to  the  computation 
that  the  navigator  can  make  to  infer  the  conrection,  depending 
on  the  form  and  substance  of  the  information  that  is  trans¬ 
mitted  to  him.  !nie  major  constraints  on  effecting  the 
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correction  are  the  maximum  allowed  real-time  data  rate  dis¬ 
cussed  In  Sectlcxi  ^.2.2  and  the  maximum  user  con^uter 
capability  discussed  In  Section  ^.2.3.  Naturally,  the  most 
desirable  cospxtatlon  would  entail  the  minimum  cost  associated 
with  both. 


In  order  to  carry  out  a  navigation  solution,  a 
user  requires  besides  the  observational  data,  the  e^iiemerldes 
of  the  spacecraft.  In  principle  then,  a  user  need  only  be 
additionally  supplied  with  the  electron  density  distribution 
as  a  foiu*  dimensional  function  of  space  and  time.  Given  the 
position  of  the  spacecraft,  the  electron  density  distribution 
and  an  assumed  navigator's  position, the  computation  Is  direct-* 
being  given  to  sufficient  accuracy  by  integration  of  the 
electron  density  al«ig  the  geometrical  range  path  as  discussed 
in  Section  4. 

An  alternative  technique  Introduced  here  would 
minimize  the  users'  computational  burden  relative  to  the 
direct  approach  discussed  above.  Ihe  central  theme  to  this 
alternative  approach  Is  that  at  a  central  computing  facility 
the  range  or  time  corrections  be  computed  for  every  potential 
location  of  a  navigator  for  each  satellite  once  the  electro 
density  distribution  has  been  predicted.  The  maximum  elect on 
density  occurs  at  an  altitude  of  about  350  km.  The  difference 
In  the  geographical  position  of  the  geometrical  path  to  a 
synchronoi’8  satellite  at  this  altitude  for  an  observer  on  the 
surface  of  the  earth  and  one  at  20,000  ft.  varies  from  zero 
for  a  satellite  overhead  to  10  minutes  of  arc  for  a  satellite 
on  the  horizon.  For  an  observer  at  60,000  ft.  the  difference 
increases  to  only  30  minutes  of  arc  for  a  satellite  on  the 
horizon.  It  is  doubtful  that  the  electrcxi  density  distri¬ 
bution  as  modelled  over  the  coverage  region  of  the  satellite 
would  be  sensitive  to  positional  differences  of  these  magni¬ 
tudes.  Ccxisequently,  the  computaticxj  of  the  Ionospheric 
correction  Is  Insensitive  to  the  navigator's  height  so  that 
the  correction  for  each  satellite  can  be  expressed  as  a  three- 
dimensional  function  of  the  navigator's  geographic  location 
and  time.  This  function,  which  will  be  referred  to  as  the 
ionospheric  correction  function,  would  be  supplied  to  each 
navigator  rather  than  a  description  of  the  electron  density 
distribution  from  which  the  range  or  time  correctioi  would 
have  to  be  computed.  To  give  an  intuitive  appreciation  for 
this  approach  it  is  instructive  to  consider  a  hypothetical 
ionosphere  in  which  the  electron  density  distribution  is 
radially  synaetric.  In  this  case,  the  contours  of  constant 
magnitude  correction  would  be  concentric  circles  about  the 
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subsatellite  point  in  which  the  magnitudes  associated  with 
the  contours  increase  with  distance  from  the  subsatellite 
point. 

5.4  TRANSiER  OF  lONOSHIERIC  CORRECTION 

There  are  several  possible  ways  by  which  the 
ionosjiieric  correction  function  can  be  transmitted  to  the 
navigator.  The  alternatives  depend  on  both  his  operaticxial 
and  accuracy  requirements. 

As  discussed  in  Section  5*2.4,  there  are  users  for 
whom  long-term  predictions  (3-6  months)  may  be  adequate.  It 
would  be  possible  to  supply  these  corrections  to  the  navi¬ 
gator  before  a  mission  is  undertaken.  An  estimate  of  the 
volume  of  a  priori  information  that  a  navigator  must  have 
access  to  can  be  obtained  as  follows.  The  satellite  con¬ 
stellations  under  consideration,  as  discussed  in  Section 
5.2.1,  are  such  nhat  the  eiiiemeri-s  for  each  satellite 
nearly  repeats  itself  every  day.  Suppose  that  a  monthly 
median  electron  density  distributioi  given  as  a  function  of 
time  of  day  and  geographic  location  is  estimated  by  the 
conqputational  center.  'Xhen,  the  ionospheric  correction 
functions  based  on  this  electron  density  distribution  for 
one  day  would  be  valid  for  every  other  day  of  that  month.  The 
assvm^ion  that  the  ionospheric  correction  function  need  be 
determined  every  two  hours  in  order  to  provide  adequate 
interpolation  within  the  day,  would  result  in  the  need  of 
12  such  functions  per  spacecraft  per  month.  A  constellaticKi 
of  5  spacecraft  for  a  period  of  6  months  would  require  a 
total  of  360  ionospheric  correction  functions.  Further 
economy  could  be  obtained  by  displaying  only  the  differences 
in  the  correction  for  the  most  desirable  four  satellites  in 
view.  This  would  result  in  a  total  of  2l6  icaiospheric 
correcticxi  functions.  By  having  this  a  priori  information, 
any  navigator  would  have  available  a  backup  capability  to 
a  real-time  scheme.  It  should  be  noted  here  that  it  would 
be  most  efficient  to  include  as  part  of  the  correcticHi  the 
influence  of  the  troposphere.  The  long-term  predictions 
would  be  most  useful  to  those  who  desire  navigation  at 
minimum  cost  and  where  extreme  accuracy  is  not  important. 
These  correction  functions  could  be  provided  in  graphical 
form  either  in  a  book  or  cxi  microfilm,  for  easy  use. 
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Short-tem  estimates^  made  available  1  to  14  days 
in  advance,  can  be  transferred  by  several  means.  For  exau^le, 
overall  corrections  to  the  long-term  estimates  con^lsing 
Just  a  few  numbers  possibly  transferred  along  with  the  space¬ 
craft  e|2iemeris  may  be  adequate. 

It  Is  the  real-time  estimates  of  the  corrections 
that  are  to  be  transmitted  throu^  the  L-band  ranging  signal 
that  are  of  prime  isqnrtsnce  and  will  dictate  the  design  of 
the  operational  correction  system.  It  should  be  noted  that 
in  an  operational  system  it  would  be  a  strai^tforward  matter 
to  include  in  the  corrections  for  the  ionosjiiere  the  influ¬ 
ence  of  the  troposphere,  niere  are  at  least  two  forms  by 
which  it  would  be  possible  to  transfer  the  ionospheric  cor¬ 
rection  function  to  the  navigator  in  real  time.  Ihe  first 
consists  of  representing  the  function  as  a  series  expansion. 
An  estimate  of  the  data  rate  can  be  obtained  by  assuming  that 
the  structure  of  the  correction  contours  closely  resemble  the 
monthly  median  values  of  the  critical  frequency  of  the  F2 
region.  In  sudi  a  case  it  Should  take  no  more  than  ^ 
coefficients  to  represent  the  ionospheric  correction  at 
each  instant  of  time  as  a  function  of  geographic  longitude 
and  latitude  [Davies.  1$6^],  Assuming  an  additional  l8  words 
for  identifiers  and  a  precision  of  seven  bits,  for  a  range  of 
0  to  127,  gives  a  total  of  U90  bits.  For  a  linear  interpola¬ 
tion  or  extrapolation  in  time,  a  total  of  98O  bits  must  be 
transferred.  At  the  maximum  allowed  data  rate  of  10  bits/sec 
discussed  in  Section  ^.2.2,  98  sec  would  be  required  to 
transfer  the  corrections  which  ^ould  be  valid  for  at  least 
an  hour.  The  frequency  at  which  the  correction  information 
is  updated  would  be  dependent  on  the  capriciousness  of  the 
ionosphere  and  shoxxld  be  under  the  control  of  the  ground 
ccmmand  facility.  The  98  seconds  establishes  a  lower  bound 
to  the  updating  interval. 

At  any  one  time,  the  navigator  would  require  about 
UOOO  bits  of  computer  storage  to  be  able  to  correr',  the 
observations  trcm  all  four  satellites.  Ihis  is  considerably 
less  than  the  allowed  maximum  of  48,000  bits  discussed  in 
Section  ^.2.3.  The  ability  to  obtain  powers  of  the  sine  and 
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cosine  functions  of  the  nominal  latitude  and  longitude  and 
the  ability  to  multiply  and  svm  which  would  be  needed  to 
perform  the  computation  should  exist  in  the  software  capa> 
bility  required  to  perform  the  navigation  solution  inde¬ 
pendent  of  the  presence  of  a  numerical  lonosiherlc  correction. 

The  alternative  approach  which  could  be  used  in 
place  of  or  in  addition  to  the  series  representation  would 
be  to  divide  the  viewing  region  of  each  spacecraft  into  a 
matrix  and  to  transmit  individually  the  correction  averaged 
over  each  element  of  the  matrix.  In  such  a  case,  each  navi¬ 
gator  need  only  determine  in  which  region  he  is  located  and 
recover  only  that  particular  number  from  the  satellite  trans¬ 
mission.  The  received  data  would  represent  the  time  error 
or  range  error  itself  so  that  no  processing  of  the  received 
information  would  be  necessary.  An  estimate  of  the  data  rate 
in  this  case  can  be  obtained  in  the  following  manner.  Each 
spacecraft  has  a  viewing  region  of  approximately  40  percent 
of  the  earth’s  surface.  For  elements  of  5  degrees  by  5  degrees, 
1040  words  would  be  required.  If  210  identifier  words  are 
assumed  needed  then  a  total  of  12^  words  mu&t  be  transmitted. 
For  a  word  length  of  seven  bits  and  a  data  rate  of  10  bits 
per  second  it  would  take  875  seconds  (ll*'.6  min)  to  transmit 
the  correction  function  fran  each  spacecraft. 

5«5  IONOSPHERIC  CORRECTION  COMPUTATICKAL  SYSTEM 

A  schematic  of  a  suggested  con^tatlonal  system  is 
given  in  Fig.  3.1.  The  system  is  constructed  in  a  manner 
such  that  long-term,  short-term  and  real-time  estimates  of 
the  corrections  are  obtained.  In  each  case  the  estimates 
are  meant  to  be  either  the  ionospheric  correction  function 
itself  or  adjustments  to  the  function.  It  is  doubtful  that 
real-time  orbit  data  need  be  used  to  Improve  the  ephemeris 
for  the  purpose  of  ccanputing  the  ionospheric  correction. 
Ephemerldes  generated  several  weeks  in  advance  will  probably 
suffice.  In  an  operational  system  it  is  anticipated  that 
the  tropospheric  corrections  would  be  determined  in  a  similar 
manner  and  combined  with  the  ionospheric  correct iois.  Con¬ 
sequently,  a  navigator  would  be  able  to  correct  for  both 
effects  by  making  a  single  computation. 
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6.  GENERAL  APHRQACH  TO  THE  HffiDICTION  ALGCXHTTOB 


6.1  INTRQDUCTKaf 

Ihe  purpose  of  this  section  Is  to  delineate  the 
thought  that  went  Into  defining  the  prediction  algorithms 
that  were  ultimately  evaluated.  As  described  in  the  previous 
sections^ both  the  satellite  ephemerldes  and  the  global 
electron  density  dlstrlbutlon^or  its  eg.ulvalent>  must  be 
predicted  into  the  future.  The  accuracy  to  which  the  sat¬ 
ellite  ephemerldes  must  be  estimated  Is  within  the  state-of- 
the-art,  as  discussed  in  Section  3*5>  The  difficult  aspect 
of  the  Ionospheric  prediction  algorithm  will  be  the  estima¬ 
tion  of  the  Ionospheric  characteristics.  For  this  reason  It 
will  be  Instructive  to  begin  with  a  concise  description  of 
the  state  of  knowledge  of  the  Ionosphere  that  is  relevant  to 
the  estimation  of  the  ionospheric  group  velocity. 

6.2  STATE  OF  KNGWIEDGE  OF  THE  IONOSPHERE 

The  Ionosphere  may  be  defined  as  that  portion  of 
the  earth's  upper  atmosphere  where  Ions  and  electrons  are 
present  In  quantities  sufficient  to  affect  the  propagation 
of  radio  waves.  It  extends  upwards  from  about  ^0  km  with 
no  well-defined  upper  boundary.  The  ionosphere  is  character¬ 
ized  by  several  regions  which  are  designated  as  D,  E,  Fl  and 
F2.  Since  the  early  years  of  ionospheric  research,  theoretical 
models  starting  with  the  well-kncwn  Chapman  model  have  been 
developed  to  attempt  to  understand  the  mechanisms  involved. 

More  con5>lex  theoretical  models  consist  of  numerically  solving 
slmultaneoxisly  the  time-dependent  coupled  chemical,  dynamical 
and  thermodynamical  partial-differential  equations  for  a  multi- 
ccxistitxjent  ionosphere  in  a  mviltlconstituent  neutral  atmos- 
ihere,  e.g.  King  and  Raster  [19711*  Their  soluticms  are 
hindered  by  many  uncertainties  such  as  the  time  and  spatial 
distrlbuticxi  of  the  constituents,  neutral  gas  temperatures, 
ion  and  electron  velocities,  neutral  wind  velocities,  produc¬ 
tion  and  loss  mechanisms,  electric  field  intensity,  geomagnetic 
field  intensity  and  in  the  boundary  conditions.  In  spite  of 
this,  theoretical  models  provide  viseful  tools  for  research 
into  the  physical  concepts  of  ionospheric  behavior.  Unfor¬ 
tunately,  the  departure  of  observations  from  the  model  results 
are  generally  the  nile  rather  than  the  exception. 

Paralleling  the  development  of  theoretical  ionospheric 
models  has  been  the  development  of  what  are  kncwn  as  model 
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icaiospheres  which  are  dedicated  toward  providing  representa¬ 
tive  characteristics  for  given  conditions  and  locations. 

Though  they  draw  heavily  on  theoretical  modeling,  model 
ionospheres  are  visually  eo^irlcal  in  nature.  The  Institute 
of  Teleccmmunlcatlons  Sciences  (ITS)  of  the  Office  of  Tele¬ 
communications,  U.S.  Department  of  Ccmimerce,  has  developed 
extensive  capahllities  in  the  predictions  of  global  iono¬ 
spheric  characteristics.  Their  estimates  are  given  as  functions 
of  the  geographic  location,  season,  time  of  day,  geomagnetic 
characteristics  and  the  12 -month  running  average  of  the 
smoothed  Zurich  sunspot  number  (R12  or  SSN)  [Jones  and  Crallet, 
1962]  and  [Jones  et  al. .  1969].  The  basis  of  these  estimates 
is  the  fitting  of  world-wide  experimental  data  obtained  over 
daily,  seasonal  and  solar  cycles.  For  example,  the  critical 
frequencies  of  the  E2  layer,  foK,  which  define  the  maximum 
electron  densities  are  derived  from  data  that  was  obtained 
from  a  world-wide  network  of  stations  over  a  span  of  several 
years.  The  number  of  different  geographic  locations  for 
which  data  was  available  varied  from  year  to  year,  averaging 
about  160.  Ihe  data  was  reasonably  well  distributed  geo¬ 
graphically,  including  such  areas  as  Evnrope,  Africa,  North 
and  South  Africa,  Asia,  Australia,  the  Pacific,  Russia  and 
mainland  China.  A  typical  global  prediction  of  the  fdF2 
available  from  ITS  is  given  in  graphical  form  in  Fig.  6.1. 

Predictions  of  the  ionospheric  characteristics 
that  are  available  from  ITS  are: 

1.  critical  frequency  of  the  F2  layer,  foF2 

2.  maximum  usable  frequency  factor,  M(3000)  F2 

3.  maximum  usable  frequency,  MUF(o)  F2 

4.  maximum  usable  frequency  at  4000  km,  MUF(4000)  F2 

5.  critical  frequency  of  the  E  layer,  foE 


6.  hei^t  to  semi-thickness  ratio  of  F2  region, 
hmF2/ymi2 


7.  semimonthly  revision  factors  to  correct  the 
MJF(O)  12 

(Oiese  predictions  are  in  use  by  several  communities.  For 
example,  they  are  currently  used  in  the  estimation  of  higii- 
frequency  telecommunication  parameters  which  are  published 
regularly  by  the  Office  of  Naval  Research  for  use  by  the  fleet. 
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Fig.  6.1  DISTRIBUTION  OF  foF2  MEDIAN  (MHz)  [JONES  et  •!.,  1969) 


6.3  ICroCDOLOGY  OF  TUB  ISEDICTION  AliQOBITHMS 

Several  alternative  approaches  exist  in  developing 
t^e  methodology  of  the  prediction  algorithms.  Xn  this  section 
two  distinct  approaches  are  discussed  and  contrasted.  In  the 
first  approach  the  vertical  group  time  delay  itself  forms  the 
hasis  of  the  prediction  in  a  manner  hy  which  it  is  possible 
to  cwnplfitely  avoid  the  e:q>licit  introduction  of  the  electron 
density  distribution.  To  facilitate  discussion,  it  is  appro¬ 
priate  to  introduce  two  terms.  Imagine  the  point  on  the 
position  vector  from  a  navigator  to  a  satellite  at  which  the 
electron  density  is  a  maximum.  Denote  this  point  as  the 
"Ionospheric  point"  and  its  projection  on  the  ground  as  the 
"ionospheric  geographic  location"  of  the  navigator.  Denote, 
also,  the  angle  between  the  intersection  of  the  position  vector 
with  the  geometrical  vertical  at  the  ionospheric  point  as  the 
"ionospheric  zenith  angle".  In  principle  it  would  be  possible 
to  obtain  a  prediction  algorithm  with  the  ionospheric  geo¬ 
graphical  location  of  the  navigator,  the  ionospheric  zenith 
angle  and  time  as  the  independent  variables  and  with  the  season, 
solar  and  magnetic  activity,  etc.  as  parameters.  This  approach 
becomes  more  attractive  by  introducing  the  notion  that  the 
significant  portion  of  the  electron  density  distribution  is 
confined  to  a  narrow  vertical  region.  For  satellites  above 
this  region  it  is  possible  to  approximately  characterize  the 
total  group  time  delay  by  the  vertical  group  time  delay  at  the 
ionospheric  geographic  location  of  the  navigator.  Ilhe  total 
time  delay  for  propagation  -through  this  point  at  other  than 
-vertical  incidence  could  then  be  Inferred;  possibly  by  Just 
multiplying  the  -vertical  time  delay  by  the  secant  of  the  iono¬ 
spheric  zenith  angle.  With  this,  -the  prediction  algorithm 
reduces  to  a  prediction  of  the  vertical  time  delay  with  -the 
ionospheric  geographic  location  of  the  na-vlgator  and  time  as 
the  independent  -variables.  Real-time  predictions  can  be 
achie-ved  by  using  real-time  measurements  of  -the  model  parame-ters 
and  the  group  time  delays  obtained  at  various  sites.  The  basis 
of  the  second  approach  is  -the  prediction  of  the  electron  density 
distribution  from  which  the  group  time  delay  or  equl-valent  range 
error  can  be  determined  in  a  separate  step.  The  time  deley 
would  be  computed  utilizing  Eq.  (4.1)  In  which  N(^,t)  is  given 
by  the  predicted  electron  density  distribution. 

The  distinction  between  the  -two  approaches,  while 
seemingly  subtle,  is  nevertheless  significant.  It  is  our 
belief  that  -this  latter  approach  is  superior  to  the  former. 
Consequently,  the  algorithms  de-veloped  here  predict  the 
electron  density  distribution  as  a  function  of  space  and 
time.  Once  this  has  been  achie-ved  it  is  stral^tforward  to 
compute  the  total  time  delay  for  each  satellite  as  a 
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function  of  time  and  the  actual  geographic  location  of 
potential  navigators  as  discussed  in  detail  in  Sections  5,3 
and  5.^. 

The  decision  to  predict  the  electron  density  distri¬ 
bution  instead  of  the  total  vertical  time  delay  has  important 
consequences  and  is  clearly  in  need  of  additional  discussion. 

It  is  our  contention  that  in  the  implementation  of  the  iono¬ 
spheric  correction  system,  the  electron  density  distribution 
should  be  introduced  for  several  reasons.  If  this  is  true, 
then  rather  than  introduce  it  implicitly,  it  would  be  superior 
to  treat  it  explicitly  in  the  prediction  algorithm.  The 
following  constitutes  several  of  the  advantages  of  this  approach: 

6.3.1  Increased  Accturacy  in  Estimation  of  the  Total  Time  Delay 

Once  the  electron  density  distribution  is  predicted, 
it  is  possible, as  explained  in  Section  5.3, to  determine 
directly  the  total  time  delays  and  equivalent  range  errors 
once  the  position  of  the  spacecraft  has  been  estimated.  This 
procedure  Inherently  takes  into  account  the  horizontal  and 
vertical  gradients  of  the  electron  density  distribution.  If 
only  the  vertical  time  delay  is  predicted  then  it  will  be  nec¬ 
essary  to  infer  from  this  the  total  time  delay  between  the 
observer  and  each  satellite.  The  uncertainty  in  this  process 
has  been  investigated  by  Woodford  and  Dutcher  [197O]  >dio  found 
that  for  an  elevation  angle  of  5  degrees  errors  as  large  as 
±  25^  possible.  To  investigate  this  point  further,  a 
limited  study  was  undertaken.  It  consisted  of  the  computation 
of  the  total  time  delay  for  several  navigators  whose  geographic 
locations  were  different  but  whose  ionospheric  geographic  loca¬ 
tions  were  identical.  Comparisons  were  then  made  of  the  verti¬ 
cal  time  delay  and  the  equivalent  for  oblique  incidence  which 
was  obtained  by  dividing  the  total  oblique  delay  by  the  secant 
of  the  ionospheric  zenith  angle.  The  model  used  for  the  elec¬ 
tron  density  distribution  is  described  in  detail  in  Section  6.4. 
For  an  elevatlcm  angle  of  5  degrees,  errors  as  large  as  44^ 
during  the  evening  and  24^{  dinrlng  the  day  were  frequently  obser¬ 
ved.  For  an  elevation  angle  pf  30  degrees  these  errors  reduce 
to  34^6  and  124  respectively. 

6.3.2  Greater  Theoretical  Background 

A  wealth  of  theoretical  studies  al^ady  exist  which 
can  be  utilized  in  the  estimation  cxT  the  electron  density 
distribution.  In  addition,  there  is  considerable  Interest 
in  understanding  and  in  predicting  the  state  of  the  ionosphere. 
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i.e,,  the  numbers  of  free  electrons  and  the  various  ionic 
species,  from  many  different  viewpoints.  As  discussed  in 
Section  6.2,  besides  the  theoretical  ionospheric  models  that 
have  been  defined,  several  semiemplrlcal  model  ionospheres 
exist  which  characterize  the  ionosphere  as  a  function  of 
geographic  location  and  time  [Bauer.  1971]. 

6.3.3  Greater  Utility 

There  would  be  many  potential  uses,  besides  navi¬ 
gation,  of  estimates  of  the  electron  density  distribution. 
These  include  use  for  over-the-horizon  radar  as  well  as  the 
avoidance  of  detection  by  this  technique  and  the  estimation 
of  high  frequency  communication  parameters.  A  significant 
number  of  the  navigators  may  have  need  for  ionospheric  infor¬ 
mation  for  purposes  other  than  navigation. 

^•3.^  Greater  Experimental  Background 

The  development  of  techniques  to  estimate  the 
electron  density  distribution  has  had  stimuli  other  than 
satellite  navigation  systems.  The  utilization  of  the  semi- 
empirical  model  ionospheres  that  have  already  been  developed 
is  straightforward.  These  models  have  been  developed  on  the 
basis  of  many  different  types  of  data  over  many  decades  over 
all  parts  of  the  earth.  At  the  puresent  time,  there  is  no 
existing  world-wide  model  for  the  vertical  Integrated  electron 
content  which  is  needed  to  determine  the  vertical  time  delays 
directly.  Since  the  launching  of  satellites,  it  has  been 
possible  in  principle  to  obtain  the  integrated  electron 
ccmtent  along  the  path  from  the  observer  to  the  spacecraft. 
Differential  time  delay  on  coherently  related  frequencies 
would  be  the  most  direct  method  but  little  of  this  data  is 


available.  There  has  been  considerable  effort  expended  in 
the  measurement  of  Faraday  rotation  data  from  wnich  it  is 
possible  to  infer  the  integrated  electron  cmtent  by  using 
normalized  models  of  the  electron  density  distributicxi  as 
discussed  in  Section  4.2.  Unfortunately,  the  observations 
that  have  been  obtained  have  been  over  limited  regions  of  the 
earth.  To  devel^  a  world-wide  capability,  it  would  be 
necessary  to  collect  data  over  long  periods  of  time  over  all 
parts  of  the  earth.  Tbis  task,  aside  from  its  magnitude, 
would  be  difficult  becaxise  of  the  inaccessibility  of  certain 
regions  because  of  geographical  or  political  impediments. 


it 


6.3.5  Increased  Confidence  In  the  Validity  of  an  Evaluation 

Care  must  be  taken  in  perfomlng  any  evaluation 
in  order  for  it  to  be  meaningful,  A  most  important  consider¬ 
ation  is  to  ensure  that  the  data  utilized  to  make  predictions 
are  completely  distinct  from  the  data  utilized  to  perform 
the  evaluation.  In  the  prediction  of  the  electron  density 
distribution  this  poses  no  problem  since  the  ionospheric 
models  can  be  based  on  distinctly  different  types  of  data 
than  those  supplied  for  use  in  the  evaluation.  However, 
difficulty  is  encountered  in  performing  a  nonprejudlced 
evaluation  when  estimating  directly  the  vertical  time  delay. 
In  the  absence  of  vertical  integrated  columnar  content  data, 
Faraday  rotation  observations  from  high  altitude  satellites 
can  be  utilized  to  determine  empirical  models  of  the  verti¬ 
cal  time  delay  as  disr.ussed  by  Da  Rosa  [  I969] .  The  body  of 
data,  described  in  Section  7,  that  was  provided  to  perform 
the  evaluation  was  Faraday  rotation  data.  This  data  set 
constitutes  a  significant  portion  of  the  total  body  of  data 
that  is  available  to  obtain  empirical  model  parameters  in 
the  prediction  of  the  vertical  time  delay. 


6,3,6  Increased  Real-Time  UTdating  Capability 


I 


Real-time  updating  of  either  approach  is  possible 
by  using  the  range  residuals  obtained  at  sites  of  known 
location.  To  attempt  to  utilize  Faraday  rotation,  differen¬ 
tial  doppler,  topside  and  bottomslde  sounding  data,  etc. 
would  require  the  introduction  of  the  electron  density  distri¬ 
bution.  Consequently,  it  wouM  be  most  desirable  to  adopt 
an  approach  which  can  directly  utilize  many  types  of  data, 
some  of  which  is  already  being  obtained  for  other  purposes, 
for  both  Increased  acciuracy,  lower  operating  costs  and  for 
a  backup  capability. 

6.4  ALGORITHM  I 
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The  ihilosophy  behind  this  algorithm  was  to  utilize 
the  best  available  long-tezmi  a  priori  predictions  of  ionospheric 
characteristics  in  order  to  predict  an  electron  density  distri¬ 
bution.  Th-n  the  ionospheric  ^oup  time  delay  would  be  com¬ 
puted  as  in  Eq.  (4.l)  where  N(r,t}  is  given  by  the  predicted 
electron  density  distribution.  The  predicted  values  of  group 
time  delay  obtained  in  this  way  would  establish  a  basis  for 
comparison  for  both  short-term  and  real-time  predictions  and 
indicate  a  measure  of  the  extent  to  which  either  is  needed.  To 
assist  in  characterizing  the  electron  density  distribution  a 
vertical  electron  density  profile  was  utilized.  The  model  adopted 
is  that  defined  by  Haydon  and  Lucas  [  I968]  which  is  depicted  in 
Fig,  6.2.  The  vertical  profile  is  characterized  by  eleven 
parameters  each  of  which  can  be  a  function  of  geographic 
location,  time,  solar  activity,  etc.  The  validity  of  the 
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model  has  been  established  by  the  realization  of  good  agree¬ 
ment  between  virtual  hei^t  computations  obtained  using  the 
model  with  experimentally  determined  values  [Haydon  and  Lucas, 
1968].  Typical  comparisons  for  several  different  conditions 
are  given  in  Fig.  6.3.  In  this  algorithm,  of  the  eleven  para¬ 
meters  that  are  required  to  define  the  vertical  electron 
density,  four  are  obtained  through  the  predictions  available 
from  ITS.  Of  these  four,  three  define  the  F2  region  which 
is  dominant  in  the  determination  of  the  integrated  electron 
content  along  the  propagation  path.  The  eleven  parameters 
are  determined  for  Algorithm  I  in  the  following  manner, 
where  the  item  numbers  correspond  to  the  numbers  on  Pig.  6.2: 

(1)  Density  at  5OOO  km  -  an  analytical  function 

of  the  month  of  the  year  and  the  gyrofrequency 
[Smith.  1961]. 

(2)  Density  at  1000  km  -  an  analytical  function 
of  magnetic  dip  and  solar  zenith  angle 
[Chandra  and  Racaswamv.  I967]. 

(3)  Maximum  electron  density  of  the  E2  region  - 
obtained  as  a  function  of  R12  ,  geographic 
location,  month  and  time  of  day  from  predic¬ 
tions  of  foF2  which  are  available  monthly 
from  ITS  [Jones  and  Pallet.  1962], 

(4)  Height  of  the  maximum  of  the  E2  region  - 
obtained  as  a  function  of  R12  ,  geographic 
location,  month  and  time  of  day  from  predic¬ 
tions  of  M(3000)P2  which  are  available 
monthly  from  US  [Jones  and  Gallet.  I962]. 

(^)  Semithickness  of  the  F2  region  -  obtained  as 
a  function  of  R12  t  geographic  location, 
month  and  time  of  day  from  predictions  of 
hmP2/ymF2  which  are  available  mcxithly  from 
ITS  [Jones  and  Gallet.  I962].  ^ 

(6)  Maximum  electron  density  of  the  FI  region  - 
obtained  as  a  function  of  R^  ,  month  and 
time  of  day  from  an  empiricanfomula 
[AVCO.  19631. 

(7)  Height  of  maximum  of  the  FI  region  -  obtained 
as  a  function  of  time  of  day  [AVCO.  I963]. 
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(8)  Semithickness  of  the  PI  region  -  constant  at 
50  km. 

(9)  Maximum  electron  density  of  the  E  region  - 
obtained  as  a  function  of  R]2  >  geographic 
location,  month  and  time  of  day  from  predic¬ 
tions  of  foE  which  are  available  monthly 
from  ITS  f Jones  and  Pallet.  I962]. 

(10)  Height  of  maximum  of  the  E  region  -  constant 
at  130  km. 

(11)  Semi thickness  of  the  E  region  -  constant  at 
20  km. 

6.5  AlfiORITHM  II 

Algorithm  II  was  devoted  to  the  capability  of 
effecting  a  real-time  ionospheric  correction.  The  most 
desirable  approach  would  have  been  to  utilize  one  of  the 
several  rather  sophisticated  theoretical  ionospheric  models 
that  are  available.  With  such  a  model  as  a  basis,  the 
philosophy  of  the  updating  scheme  would  be  rather  straight¬ 
forward.  It  would  be  possible  to  determine  some  of  the  more 
uncertain  model  parameters  as  a  function  of  time  by  "fitting” 
the  model  to  whatever  experimental  ionospheric  data  is  avail¬ 
able,  Prediction  of  the  parameters  and  Integration  of  the 
system  of  nonlinear  time -dependent  partial  differential 
equations  would  then  result  in  the  detenu inat ion  of  the 
electron  density  distribution  at  any  time  in  the  future. 

The  evaluation  of  the  performance  of  any  of  the 
algorithms  developed  in  this  study  was  to  be  achieved  by 
processing  a  large  amount  of  Faraday  rotation  data  obtained 
over  several  years.  This  data  set  and  the  evaluation  plan 
is  defined  in  detail  in  Sections  7  and  8,  Preliminary 
estimates  for  the  cost  of  processing  even  only  a  minimum 
definitive  subset  of  data  Indicated  that  the  computational 
complexities  had  to  be  minimized.  Consequently,  in  order 
to  perform  a  meaningful  evaluation  of  the  algorithm  a  less 
sophisticated  approach  had  to  be  adopted. 

The  a  priori  predictions  of  Algorithm  I,  as 
described  in  the  previous  section,  constitute  in  themselves 
a  rather  sophisticated  model.  The  predicted  electron  density 
distribution  has  considerable  geogrephical  structure  as 

evidenced  by  a  typical  map  of  the  foFS  distribution  (which  is 
proportional  to  the  square  root  of  the  maximum  electron  density) 
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given  In  Fig.  6.1.  Because  of  this,  it  was  concluded  that 
rather  than  atten^rt;  to  model  the  electron  density  distribution 
Itself  it  would  be  more  expedltioxis  to  model  the  differences 
between  the  actual  electron  density  distribution  and  the 
estimates  obtained  from  the  a  priori  predictions  of  Algorithm  I 
The  basis  for  this  approach  is  the  anticipation  that  the  a 
priori  estimates  will  adequately  reflect  the  detailed  structure 
and  that  a  correction  which  is  slowly  varying  in  both  space 
and  time  would  be  adequate. 

The  adjunct  model  developed  for  the  electron  density 
distribution  is  distinguished  by  being  an  analytic  model  con¬ 
sisting  of  distinct  daytime  and  nighttime  functions.  The  day¬ 
time  model  was  obtained  from  a  first  order  perturbation  solu¬ 
tion  to  the  continuity  eqxiation  in  which  the  production  is 
given  by  the  Chapman  function,  the  loss  mechanism  assumed  as 
a  linear  function  of  the  electron  density  and  where  transport 
phenomena  are  neglected.  Under  these  c(»ditlons,  the  electron 
density  at  the  hel^iit  of  maximum  production,  which  is  assumed 
to  be  the  hel^iit  of  the  maximum  density,  is  given  by 

N,^  =  a  cos  X  +  B  tt)  cos  X  cos  X  sin  <«)(t-12)  +  y  (6.1) 

Cl  ns 

where 

a,P,Y  =  parameters,  m 

X  -  solar  zenith  angle  at  "ionospheric  point" 

X^  =  "ionospheric  latitude"  of  navigator 

X  =  latitude  of  sun 

8 

tt»  =  earth's  rotaticai  rate,  rad/hour 

t  =  local  time,  hour 

=  maximum  daylight  electron  uenslty,  m"^ 

nils  expression  for  the  maximum  electron  density  exhibits 
some  of  the  more  typical  characteristics  that  have  been 
observed.  The  first  term  alone  would  give  an  electron 
density  that  is  a  maximum  at  the  subsolar  point  and  decreases 
with  increasing  solar  zenith  angle.  The  second  term  intro¬ 
duces  a  phase  (or  time)  lag  in  the  appearance  of  the  noirimum 
which  conforms  with  observations.  One  effect  that  this  ex¬ 
pression  does  not  Introduce  is  diffusion  of  the  electrons 
along  the  geomagnetic  field  lines.  This  phenomenal  is  dis¬ 
cussed  by  Rastogi  [  19^91  who  demcxistrated  that  it  is  signifi- 


-38- 


can't  in  the  'vicinity  of  the  geomagnetic  eg,uator  but  can  be 
neglected  at  the  middle  latitudes.  The  magnitude  of  the 
effect  is  dxcwn  to  be  correlated  with  both  magnetic  dip  angle 
and  solar  activity.  In  order  to  account  for  this  phenonenon 
the  following  approxlmaticm  was  used  which  reasonably  fits 
the  experimental  data  given  by  Rastogi  [ I959] ; 


“a  =  "a 


(6.a) 


where 


K(D)  = 


,  1D|  iTT/4 


(1-0.305  cos  6  D)  ,  |d1  <  n/U 


D  =  magnetic  dip  at  a  height  of  350  km>rad 
a  maximum  daylight  electron  density,  m”^ 

The  ni^ttime  model  that  was  utilized  represents 
a  solvttlon  to  the  continuity  equation  which  consists  of  a 
linear  loss  mechanism  and  no  production  mechanism.  The 
maximum  nighttime  electron  density  is  then  Just  siiiq>ly 


=  N,(t  „)  e 
n  d'  ss 


-6(t-t  „) 

+  Njl-e  ®®  )  (6.3) 


where 


=  ni^ttime  maximum  electron  density,  m  » 

=  parame-ter  to  be  de'fcermlned,  m 
N,(t  )  =  maximum  electrcxi  density  at  sunset,  m”^ 

Cl  S8 

t  =  local  time,  hour 

t  =  time  of  sunset,  hovir 
ss 

6  »  effective  recombination  coefficient,  12/day 

This  solution  represents  an  exponential  decay  of  the  maximum 
density  from  sunset  to  sunrise. 

Now  that  the  maximum  electron  density  as  a  function 
of  time  and  geograihical  location  has  been  specified,  it  is 
necessary  to  define  the  'vertical  profile  as  a  function  of 
this  single  variable.  Let 


I  i 


N(h) 


r  w(h) 

L“n  S(») 


where  h  Is  the  height  above  the  earth's  surface  and  w(h) 
and  g(h)  represent  normalized  vertical  profiles.  There 
are  many  alternatives  available  for  both  w(h)  and  g(h) 
but,  as  will  be  shown  later  in  Section  8. 5. 2, they  will  not 
have  to  be  explicitly  defined  for  this  particular  study. 


In  summary,  an  analytic  model  for  the  electron 
density  distribution  has  been  specified  as  a  function  of 
geographic  location,  local  time  and  magnetic  dip.  The  model 
is  a  linear  function  of  the  parameters  which  are 

to  be  determined  frcxn  real-time  observaticais  and  other 
measures  of  the  ionosphere.  Because  of  the  D.inearlty  of  the 
model,  it  is  possible  to  utilize  these  functions  to  model 
directly  the  difference  between  the  observed  icaiosphere  and 
the  a  priori  predictions  of  Algorithm  I. 

6.6  ALGORITHM  III 


To  further  explojre  the  potential  of  a  real-time 
icxiospheric  correction  an  alternative  approach  to  Algorithm  II 
was  explored.  In  Algorithm  II  the  basis  of  the  correction 
was  the  development  of  a  theoretical  ionospheric  model.  In 
contrast,  the  basis  of  Algorithm  III  is  the  utilization  of 
an  en^irlcally  developed  model  ionosphere.  At  any  geographic 
location  the  vertical  electron  density  profile  is  taken  to 
be  that  utilized  in  Algorithm  I  and  depicted  in  Fig.  6.2. 

As  described  in  Section  6.4,  five  of  the  eleven  parameters 
used  to  define  the  profile  for  any  specific  geographic 
position  are  functions  of  the  predicted  •  Most  signi¬ 
ficant  are  the  three  parameters  defining  the  F2  regic»i 
which  contributes  overwhelmingly  to  the  total  group  time 
delay.  These  are  functions  of  geographic  location,  time, 
season,  and  R22  •  World-wide  predictions  of  each  of  the 
parameters  can  be  obtained  from  ITS  as  discussed  in  Section 
6.3.  With  this  background  it  is  now  possible  to  describe 
the  central  theme  of  Algorithm  HI. 

The  R22  used  in  the  ITS  predictions  represent 
the  predicted  12-month  running  average  of  the  smoothed 
Zurich  sunspot  number  determined  from  optical  observations 
of  the  sun.  Suppose  this  interpretation  is  dianged  so  that 
the  R22  to  simply  represent  a  number,  or  parameter. 
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on  which  the  electron  density  distribution  is  dependent. 
Given  ionosiheric  data  for  some  interval  of  time,  then,  it 
would  be  possible  to  determine  an  effective  Rj2  (designated 
by  R)  such  that  the  modeling  electron  density  distribution 
best  represents  the  observed  data.  Then  by  using  this  R  , 
or  a  value  predicted  from  a  sequence  of  such  numbers,  a 
real-time  prediction  of  the  electron  density  distribution 
can  be  obtained. 

Inherent  to  this  approach  are  the  spatial  and  time 
correlations  which  have  been  determined  experimentally  and 
reside  in  the  sets  of  coefficients  which  form  the  basis  of 
the  ITS  predictions.  One  concern  in  this  approach  is  that 
the  ITS  model  ionosphere  represents  monthly  median  values 
and  the  same  correlation  may  not  exist  for  the  Instantaneous 
electron  density  distribution .  Utilizing  data  over  a  period 
of  one  to  several  days  should  minimize  any  difficulty 
Induced  by  this  Inconsistency.  The  Implementation  of  this 
approach  is  discussed  in  detail  in  Section  8.^.3. 


r 


7.  EVALUATION  DATA 


7.1  INTRODUCTION 

The  experimental  data  provided  for  use  in  the  study 
were  Faraday  rotation  observations  obtained  from  VHF  (nominally 
137  mz)  transmissions  from  near-synchronous  satellites.  From 
this  type  of  data  j.t  is  possible  to  "infer"  an  "experimental" 
vertical  ionospheric  group  time  delay,  data 

obtained  in  this  manner  was  to  be  used  for  updating  and  evalu¬ 
ation  of  the  prediction  algorithms.  In  Section  4.2  theoretical 
expressions  for  converting  the  experimental  Faraday  rotation 
data  to  T»  data  are  given  and  the  approximations  in  these 
e 

expressions,  which  would  result  in  errors  in  t  ,  are  dis- 

cussed.  In  this  section,  the  errors  in  resulting  frcmi 

experimental  errors  in  the  Faraday  rotation® and  the  defici¬ 
encies  of  the  data  base  as  regards  its  use  in  verification 
of  algorithms  for  use  in  a  navigation  context,  are  discussed. 

7.2  EXPERIMENTAL  ERRORS 

The  experimental  errcars  in  the  Faraday  rotation 
data  result  from  uncertainties  in 

a)  the  initial  polarization  of  the  satellite 
radiated  signal 

b)  the  msasurement  of  the  polarization  of  the 
received  signal 

c)  resolution  of  the  nw  uncertainty  inherent  to 
any  2>olarlzation  measurement. 

Together,  a)  and  b)  have  been  estimated  by  the  experimenters 
to  amount  to  about  ±  1^°  which  at  I6OO  HRz  would  range  frm 
^  0.2  to  ±  0.8  ns  depending  on  station  location  and  time  of 
day.  The  time  of  day  dependence  arises  throu^  the  diurnal 
variation  of  M(6,t}  ,  as  discussed  in  Section  4.2.  As 
regazds  c),  it  is  usually  assumed  that  the  nn  ambiguity 
has  been  successfully  resolved.  However,  the  existence  of 
negative  values  of  rotation  in  the  data  make  the  assun^lon 
suspect.  An  error  of  ±  tt  in  resolving  this  ambiguity  can 
lead  to  errors  in  the  range  ±2  to  :t4  nsec  at  I6OO  tHz 
for  daylight  periods, depending  on  station  location.  In  this 
study  it  has  been  assumed  that  the  values  of  the  Faraday 
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rotation  data  that  have  been  supplied  are  correct.  Thus 
the  Inferred  Ty^  data  are  subject  to  the  experimental 

error  in  the  Faraday  rotation  measurement. 

7.3  IIMITmONS  OF  TEIE  UNSA  BASE 

I 

One  of  the  deficiencies  of  the ‘data  base  in  pro¬ 
viding  a  meaningful  evaluation  of  ionospheric  prediction 
algorithms  in  a  navigation  environment  is  its  limited  geo¬ 
graphic  and  geomagnetic  coverage.  Table  7.I  gives  a  list 
of  the  data  provUed  for  use  in  the  study.  Table  7.2  gives 
the  geographic  coordinates  and  geomagnetic  dip  angles  of 
each  of  the  measurement  sites.  Fig.  7*l'is  a  world  map 
depicting  the  measurement  sites.  Since  the  satellites  were 
in  eq.uatorial  orbits  with  slow  drift  rates,  the  data  was 
confined  to  a  limited  geographic  area.  With  the  exception 
of  Honolulu,  Cold  Bay,  and  Areclbo,  the  remaining  stations 
are  in  the  continental  thiited  States  and  at  Edmonton.  It 
is  already  known  that  there  exists  a  hi^  degree  of  coire- 
latlon  in  the  behavior  (magnitude  and  local  tenoral 
dependence)  of  the  data  at  the  mid- latitude  stations. 
Moreover,  as  is  clear  from  Fig.  7*1  Table  7.2  for  all 
stations  except  Honolulu,  the  geomagnetic  dip  angles  of 
the  ionospheric  points  fall  in  the  range  of  50°N  to  72®Nj 
Honolulu  is  at  37*5°N.  On  the  basis  of  the  behavior  of  the 
ionospheric  electron  density  as  a  function  of  geomagnetic 
dip,  one  would  expect  that  the  magnitudes  and  Ivcal  temporal 

dependence  of  at  all  stations  except  Honolulu  would 

'^e 

be  similar.  It  would  be  expected  that  the  magnitudes  of  Ty 
at  Honolulu  would  be  hi^er  than  all  other  stations.  Ihese  ® 
expectations  were  borne  out  as  evidenced  by  the  cumulative 
fi«quency  distributions  of  the  vertical  time  delay,  Ty  , 
discussed  later.  ® 

Another  problem  with  the  data  was  that  the  satel¬ 
lites  were  always  south  of  che  measurement  stations.  In  an 
actual  navigation  system,  a  navigator  would  in  general  be 
looking  at  satellites  north,  east,  south,  and  west  frexn  his 
location.  The  hi^  degree  of  correlation  esdiibited  by  the 
data  as  discussed  above  would  then  tend  to  obscure  the 
validity  of  the  evaluation. 


OIABIE  7.1 

DTJBmONS,  OBSERVATION  STATIONS,  SATELLITES  AND  SOURCES  OF 
FARADAY  ROTATION  DATA.  (AFCRL,  SSL,  UH  ani  UI  represent 
Air  Force  Cambridge  Research  Lahoratories,  Stanford  Elec> 
tronics  Laboratories,  IMiversity  of  Hawaii  and  University 
of  Illinois,  respectively.) 


Duration 

Location 

Satellite 

Source 

1  Jan  1965 

Stanford.,  Ca. 

SINCOM3 

SSL 

to  31  Dec  1965 

Honolulu,  Hi. 

SYNC0M3 

SEL 

1  Dec  1967 

Arecibo,  P.R. 

ATSl 

SEL 

to  18  Apr  1968 

Honolulu,  Hi. 

ATSl 

SEL  and  UH 

Sagamore  Hill,  Ma. 

ATS3 

AFCRL 

Stanford,  Ca. 

ATSl 

SEL 

Stanford,  Ca. 

ATS3 

SEL 

Urbana,  II. 

ATS3 

UI 

1  Jan  1968 

Honoltilu,  Hi. 

ATSl 

UH 

to  31  Dec  1968 

Sagamore  Hill,  Ma. 

ATS3 

AFCRL 

Stanford,  Ca. 

ATSl 

SEL 

IDec  1968 

Arecibo,  P.R. 

ATS3 

SEL 

to  13  Nov  1969 

Clark  Lake,  Ca. 

ATSl 

SEL 

Cold  Bay,  Ak. 

ATSl 

UI 

Edmonton,  Can. 

ATSl 

SEL 

Fort  Collins,  Co. 

ATSl 

SEL 

Honolulu,  Hi. 

ATSl 

UH 

Rosman,  N.C. 

ATS3 

SEL 

Sagamore  Hill,  Ma. 

ATS3 

AFCRL 

Stanford,  Ca. 

ATSl 

SEL 

Stanford,  Ca. 

ATS3 

SEL 

Urbana,  11. 

ATS3 

UI 

TABLE  7.2 


GEOMAGNETIC  DIP  AND  GEOGRAPHIC  COORDINATES  OF  THE  "lONOSHlERIC 
GEOGRALHIC  LOCATION"  OF  THE  FARADAY  ROTATION  MEASUREMENT  STATIONS 

(Epoch  Janiiary  15,  1969) 


Geoioagiietic  Dip  Latitude  longitude 
Station  Satellite  degrees  north  degrees  north  degrees  east 


Edmonton 

ATSl 

71*3 

48.1 

239.5 

Sagamore  Hill 

ATS3 

70.6 

39.3 

289.4 

Urhana 

ATS3 

68.7 

37.0 

273.8 

Rosman 

ATS3 

65.2 

32.6 

278.4 

Cold  Bay 

ATSl 

64.0 

49.8 

199.5 

Fort  Collins 

ATSl 

63.6 

3T.0 

248.5 

Stanford 

ATSl 

58.4 

34.6 

234.8 

Clark  Lake 

ATSl 

55.6 

30.9 

240.1 

Arecibo 

ATS3 

50.0 

17.3 

.  293.1 

Honolulu 

ATSl 

37.5 

19.9 

202.7 
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THE  INCLINATION  (DIP)  OF  THE  EARTH'S  MAGNETIC  FIELD  FOR  EPOCH  1965.0 
GIVING  THE  MAGNETIC  DIP  POLES  AND  MAGNETIC  DIP  EQUATOR:  CONTOURS 
OF  CONSTANT  I  lOEG).  {AFTER  LEATON.  1971.) 


Fig.  7.1  WORLD  MAP  SHOWING  GEOGRAPHIC  AND  MAGNETIC  DIP  COORDINATES  OF  THE  "IONOSPHERIC 
GEOGRAPHIC  LOCATION"  OF  THE  FARADAY  ROTATION  MEASUREMENT  STATIONS  (EPOCH 
JANUARY  15. 1969) 


8.  EVAIUASION  IHIIOSOIHY 


8.1  EVALUATION  CONDITIONS 

In  order  to  simulate  the  application  of  the  algor- 
ithm  in  a  navigation  system  the  Earaday  rotation  data  pro¬ 
vided  for  the  study  and  listed  in  Table  7.I  was  subsetted 
into  five  so-called  "evaluation  conditions”  [Dept,  of  Air 
Force .  19701*  five  evaluation  conditions  are  described 
in  Table  8.1.  In  an  evaluation  condition,  the  stations 
listed  as  "observation  stations”  correspond  to  the  system 
ground  stations  whose  data  is  to  be  used  in  the  updating 
process.  The  stations  listed  as  "evaluation  stations" 
correspond  to  navigators.  Their  data  is  to  be  used  for 
ccmparlson  with  the  predictions  and  thus  forms  the  basis  of 
the  evaluation  of  the  algorithms. 

8.2  CRITERIA  lOR  EVALUATING  THE  AIGCRUHMS 

In  Section  6.3.5  attention  was  directed  to  the  care 
that  must  be  exercised  in  performing  the  evaluation  of  the 
prediction  algorithms.  The  algorithms  developed  in  this 
study  are  truly  prediction  algorithm  for  use  in  a  navigation 
system.  This  point  is  emphasized  by  two  in^rtant  criteria 
that  have  been  observed  in  their  eval\zationj 

1)  none  of  the  experimental  data  provided  for 
evaluation  of  the  prediction  algorithms  at 
stations  designated  as  evaluation  stations 
(navigators)  were  used  in  any  way  to  make  the 
predictions  of  the  time  delays  at  an  evaluaticxi 
station 

2)  the  experimental  data  from  the  stations  desig¬ 
nated  as  observatlcxi  stations  (system  ground 
stations)  were  used  in  a  prediction  sense 
rather  than  a  filtering  or  smoothing  sense 
(the  terms  prediction,  filtering  and  smoothing 
being  used  in  their  technical  estimation  theory 
sense). 

The  first  criterion  indicates  that  it  was  not 
appropriate  to  derive  any  correlations  between  observation 
and  evaluation  stations  and  then  vtse  these  correlatlcxis  in 
the  evaluation  of  the  algorithms.  The  in^lication  of  such 
an  approach  in  a  navigation  system  is  knowledge  of  the 
correlations  between  observation  stations  and  any  geographic 
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Pnciiv  p«i  Mart 


!EABI£  8.1 

EVAUIATIOH  CONDinC^ 

EVALUATION  CONDITION  1:  Day  3^-6?  to  Day  109-68 

Observation  Stations 

Areclbo  ATS3 
Honolulu  ATSl 

Evaluatloi  Stations 

Stanford  ATSl 
Stanford  ATS3 
Urbana  ATS3 
Sagamore  Hill  ATS3 

EVALUATICK  CONDITION  2:  Day  342-6?  to  Day  109-68 

Observation  Stations 

Sagamore  Hill  ATS3 
Areclbo  ATS3 
Honolulu  ATSl 

Evaluation  Stations 

Stanford  ATSl 
Stanford  ATS3 
Urbana  ATS3 

EVALUATION  CONDITION  3:  Day  1-65  to  Day  365-65 

Observation  Station 
Stanford  SINC0M3 
Evaluation  Station 
Honolulu  SINC0M3 

EVALUATION  CONDITION  4:  Day  1-68  to  Day  366-68 

Observation  Station 

Stanford  ATSl 

Evaluation  Statirn 

Sagamore  Hill  ATS3 
Honolulu  ATSl 
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TABIE  8.1  (cont*d) 


EVALUATION  COtJDITION  5:  Day  336-68  to  Day  317-69 

Observation  Stations 

Edmonton  ATSl 
Sagamore  Hill  ATS3 
Honolulu  ATSl 


Evaluation  Stations 


location  and  the  temporal  variations  (diurnal,  seasonal, 
solar  cycle)  of  these  correlations. 

The  second  criterion  indicates  that  in  the  real¬ 
time  updating  (Algorithms  II  and  III)>  the  observation  station 
data  is  used  only  up  to  some  time  prior  to  the  time  of 
Interest  at  evaluation  stations  (prediction)  rather  than 
up.  Including,  and  beyond  (filtering  and  smoothing)  the  time 
of  use  at  an  evaluation  station.  Thus  the  algorithms  truly 
represent  the  navigation  situation. 

8.3  UIDATING  DATA 

As  was  indicated  in  the  general  discussion  in 
Section  6.3}  e  major  virtue  of  the  estimation  of  the  electron 
density  dlstributloi  approach  that  was  adopted  in  this  study, 
is  the  diverse  body  of  data  vhich  could  be  utilized  for  the 
near  real-time  updating  of  Algorithms  n  and  III.  Such  data 
as 

1)  Integrated  electrcai  content 

2)  ground  based  vertical  and  oblique  soundings 

3)  topside  vertical  and  oblique  soundings 

4)  geomagnetic  activity 

3)  ionospheric  refraction  data  available  from 
other  systems  (e.g. ,  Navy  Navigation  Satellite 
System) 

6)  Faraday  rotation 

could  easily  be  utilized  for  purposes  of  updating. 

An  updating  scheme  using  some  or  all  of  the  above 
data  naturally  partitions  itself  into  three  types  of  updating: 

1)  global 

2)  regional 

3)  local 

In  global  updating, data  from  a  specified  geographic 
region  would  be  used  to  update  and  predict  the  ionospheric 
time  delay  at  other  geographic  locations,  both  within  and 
outside  the  region.  In  an  operational  system  this  must 
correspond  to  having  the  system  ground  stations  confined  to 
the  territory  of  the  Ihiited  States  and  using  data  collected 
by  these  stations  for  world-wide  predictions. 

In  a  regional  updating,  where  the  interest  is  con¬ 
fined  to  the  cover  region  of  a  single  satellite  constellation^ 


data  collected  at  system  ground  stations  within  the  region 
would  be  used  to  update  the  algorithms  in  order  to  predict 
the  ionospheric  time  delays  within  the  entire  region. 

A  local  updating  would  be  similar  to  the  regimal 
updating  in  which  both  the  system  ground  stations  and  navl> 
gators  are  confined  to  a  limited  aarea  within  the  region. 

Each  of  the  algorithms  proposed  in  this  study  has 
inherent  to  it  a  global  capability.  However,  in  performing 
the  study,  only  the  integrated  electron  content,  inferred 
from  the  provided  data  base,  was  used  as  a  data  source  for 
updating.  Thus, given  the  limitations  of  the  data  base  as 
discussed  in  Section  7,  the  evaluation  of  the  algorithms  was 
necessarily  local. 

8.4  VEIGHTIHG 

A  considerable  effort  was  expended  in  considering 
whether  or  not  a  weighting  scheme  should  be  implemented  as 
pert  of  the  algorithms.  When  one  considers  wei^ting  it  is 
for  the  reasons  of: 


1)  data  of  distinct  differences  in  q,uality 
(precision  and  accuracy) 

2)  deficiencies  in  the  model  with  regard  to  the 
intended  use  of  the  model. 

In  this  study  all  the  data  furnished  was  assumed 
to  be  of  equal  quality.  This  assumption  is  somewhat  suspect 
due  to  the  existence  of  negative  Faraday  rotation  values  in 
the  Honolulu  data.  However,  without  the  proper  means  of 
assessing  the  quality  of  all  the  data  provided,  it  is  the 
only  reasonable  assumption. 

All  models  suffer  from  some  deficiencies  and  one 
might  want  to  weight  the  observation  station  data  so  as  to 
improve  the  predictions  at  some  specified  sites.  This  was 
particularly  true  of  the  Honolxilu  data  which  was  of  a  dis¬ 
tinctly  different  character  when  compared  with  the  data  ftrom 
all  other  stations.  One  mi^t  wish  to  deweigh'li  the  Honolulu 
data  when  Honolulu  appeared  as  an  observation  station  along 
with  other  observatlcxi  stations. 

The  development  of  weighting  schemes  is  not  a 
trivial  matter  of  geographic  or  relative  distance  wei^ting 
but  involves  complex  considerations.  This  is  made  apparent 


. . . . . . . . . . . . 


In  evaluation  condition  vhen  one  notes  that  the  character 
of  the  ionos|iheric  time  delay  data  at  Arecibo  is  much  more 
like  that  at  Cold  Bay  than  is  the  H^olulu  data,  (see  Fig.  7.I) 
Yet  a  sinqple  relative  distance  \relghtlng  scheme  would  tend  to 
inappropriately  develght  Arecibo  relative  to  Honolulu.  From 
our  knowledge  of  the  correlation  of  ionospheric  electron  densl* 
ties  with  geomagnetic  dip.  Section  6.^,  it  is  apparent  th?t 
any  weleiitlng  scheme  would  necessarily  involve  geomagnetic 
parameters. 


Limited  comparisons  of  the  Algorithm  II  capability 
with  and  without  Honolulu  observation  data  were  made  for 
January  1969  from  evaluation  condition  Ihe  results  are 
given  in  Table  8.2.  13ie  results  indicate  that  if  Honolulu 
data  were  not  used,  the  xms  of  the  residuals  averaged  over 
the  evaluation  stations  are  about  12  percent  less  than  when 
Honolulu  data  were  used. 


In  order  to  avoid  creating  a  false  is^iiresslon  by 
developing  a  wel^iiting  scheme  peculiar  to  the  given  regional 
character  of  the  data  base,  it  was  decided  to  give  all  obser¬ 
vation  station  data  equal  wel^t*  This  was  consistent  with 
the  desire  to  Interpret  the  results  of  the  evaluation  in  a 
world-wide  context. 


8.:^  EVALUATION  DETAILS 


8.5.1  Algorithm  I 


Ihe  basic  features  of  Algorithm  I  were  given  in 
Sectim  6.4.  Therein  it  was  pointed  out  that  the  underlying 
phlloeophy  behind  this  algorithm  was  the  use  of  the  best 
available  long-tezm  estimates  of  ionospheric  (diaracterlstics 
needed  to  define  the  electron  density  dlstrlbixtion.  One  of 
the  best  long-texn  predictions  of  ionospheric  characteristics 
presently  available  are  provided  by  ITS  as  discussed  in 
Section  6.2.  Of  tbe  predictions  available.  Algorithm  I  uses 


1] 


E-layer,  maximum  electron  density 
F-layer,  maximum  electron  density 
F-layer,  hei^t  to  semi-thickness  ratio 
F-layer,  M(3000)  factor  from  which  hei^t  can 
be  inferred 


with  which  to  chtaln  four  of  the  eleven  required  parameters. 
The  long-term  predictions  of  the  above  characteristics  are 
prepared  by  ITS  five  months  prior  to  the  month  of  ^diction. 
Use  of  these  predictions  in  Algorithm  I  would  enable  one  to 
provide  a  priwi  estimates  of  world-wide  ionospheric  time 
delay  five  mcmths  in  advance  of  use.  Sucdi  estimates  would 
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TABI£  8.2 


COMPARISON  OF  AlfiOBiTHM  II  IREDICTION  CAPABILm  WITH  AND 
WITHOUT  HONOLULU  OBSERVATION  DATA 


January  1969  Evaluation  Condition  3 

RMS  of  the  Residuals 
(ns) 


Evaluation  Station 

Number 
of  Points 

With  Honolulu  Without  Hoi 
Data  Data 

Stanford-ATSl 

2021 

3.42 

3.69 

Clark  LakerATSl 

477 

3.47 

3.14 

Port  Comna-ATSl 

2084 

3.62 

2.81 

Urhana-ATS3 

2137 

3.86 

3.10 

Rosman*^ATS3 

2039 

3.30 

3.22 

Areciho-ATS3 

1847 

5.47 

4.13 

CoM  Bay^ATSl 

2105 

2.21 

2.26 

Mean  over  all  stations  — 

3.62 

3.19 

be  useful  to  certain  classes  of  navigators,  to  vhon  they  couM 
be  provided  In  advance  as  discussed  In  Section  ^.4. 

In  addition  to  the  long-term  predictions,  ITS  also 
provides  the  so-called  Semi-monthly  Revlslcm  Factors  as 
described  In  Zacharlsen.  Ostrotr  and  Huang  [  I969I .  These 
factors  apply  only  to  the  F-layer  critical  frequency  and  can 
be  used  to  provide  corrections  to  the  long-term  estimates  of 
the  F-layer  maximum  electron  density.  They  are  prepared  by 
subjective  methods  of  ccmparlson  of  the  long-term  predictions 
of  the  MtIF  (3CXX))  F2  with  experimentally  determined  values  of 
MUF  (3000)  is  from  a  small  set  of  lonosonde  stations.  The 
cooqiarlsons  are  based  on  experimental  data  available  one  to 
two  weeks  prior  to  the  half -month  of  application. 

The  Implementation  of  Algorithm  I  has  Included  the 
use  of  these  revision  factors.  IMs  was  in  keeping  with  the 
desire  to  provide  the  best  possible  estimates  of  the  ionospheric 
characteristics.  Since  the  revision  factors  aare  available 
about  two  weeks  prior  to  use.  Algorithm  I  as  in^lements  is, 
strictly  speaking,  a  two-week  prediction.  Hcwever,  comparlstais 
of  Algorithm  I  predictions  with  and  without  the  revision  factors 
shewed  little  substantial  differences. 

To  Investigate  the  significance  of  the  semi-monthly 
revision  factors,  a  compilation  of  their  statistics  for  a 
representative  saitq>le  of  the  data  supplied,  with  which  to 
perform  the  evaluation,  is  given  in  Table  6.3.  (The  indicated 
months  actually  represent  those  which  formed  the  subset  of  data 
ultimately  used  in  the  evaluation  as  dlscxissed  in  Section  9*1*) 
The  mean  and  both  zms*  are  obtained  by  summing  over  time  through¬ 
out  the  day  and  over  space  over  the  entire  globe.  The  rms*  of 
the  semi-monthly  revision  factor  over  all  the  months  is  six 
percent,  l.e.,  a  six  percent  correction  to  the  critical  frequency 
of  the  F2  region.  What  is  of  interest  however,  is  the  influence 
on  the  vertical  ionospheric  time  delay  which  is  proportional  to 
the  Integrated  electron  density.  It  is  not  unreasonable  to 
assume  that  the  vertical  time  delay  is  proportional  to  the 
maximum  electron  density  which  is  itself  proportional  to  the 
square  of  the  critical  firequency  of  the  F2  region.  For  this 
reason  the  rms  of  the  deviation  from  unity  of  the  square  of 
the  semi-monthly  revision  factor  is  given  for  each  month.  Ihe 
rms  is  12  percent  over  all  the  data. 

Based  on  this  discussion,  it  is  possible  to  make 
a  statement  as  to  the  maximum  benefit  incurred  by  utilizing 
the  semi-monthly  revlsicxi  factors;  or  alternately  to  estimate 
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TABIE  8.3 


STATISTICS  OF  SEMI-MONTHLY  REVISION  FACTORS 
Data  Points  Per  Month 


Month 

Year  Mean 

RMS  from  Unity 

RMS  of  the  Squso 

January 

65 

0.99 

0.06 

0.12 

April 

65 

1.00 

0.08 

0.16 

July 

65 

1.01 

0.09 

0.19 

January 

68 

1.03 

0.05 

0.11 

February 

68 

1.05 

0.07 

0.15 

Iferch 

68 

1.05 

0.06 

0.13 

April 

68 

1.03 

0.05 

0.11 

October 

68 

0.99 

0.05 

0.09 

January 

69 

1.02 

0,06 

0.10 

April 

69 

1.00 

O.OU 

0.07 

July 

69 

1.00 

0.02 

0.01^ 

October 

69 

0.98 

0.04 

0.07 

from  Unity 


the  maxhnum  degradation  to  be  expected  were  Algorithm  I 
utilized  without  them.  The  degradation  would  have  been  no 
greater  than  12  percent.  Consequently,  to  a  precision  better 
than  12  percent,  It  Is  possible  to  Interpret  the  results  of 
Algorithm  I  as  a  five  month  prediction. 


8.5.2  Algorithm  II 

Algorithm  II  constitutes  a  real-time  updating  of 
the  a  priori  estimates  provided  by  Algorithm  I.  The  details 
of  the  algorithm  are  given  In  Section  6.5.  Briefly,  the 
algorithm  consists  of  a  parameterized  analytic  model  which 
describes  the  electron  density  distribution  as  a  function 
of  geographic  location,  local  time,  and  magnetic  dip. 

Using  Eq.  (6.4)  the  vertical  Ionospheric  time 
delay  at  local  hour  t  Is  given  as 


T  (t)  *  C  J  N,(t)  w(h)dh  ,  for  daytime  (8.1a) 
d 

and 

(t)  «  C  J  N^(t)  g(h)dh  ,  for  nighttime  (8,1b) 


where  C  =  1.34  x  10" f  rationalized  mks  units  where  f 
is  the  carrier  frequency  of  the  ranging  signal  in  Hz, 
Inserting  Eqs.  (6.2)  and  (6.3)  Into  the  above  yields 


-58- 


In  the  Implenientation  of  the  algorithm,  a  simpli¬ 
fication  was  introduced  which  allowed  the  use  of  the  Algor¬ 
ithm  I  results  and  thereby  removed  the  need  for  additional 
ray  tracing.  The  simplification  is  based  on  the  fact  that 
the  major  contribution  to  comes  from  the  F2- layer  of 
the  ionosphere.  It  is  then  possible  to  write  that 

f 


T^(X,0,t,R)  «  C  \ax(\,0,t,R)  J  W(\,e,t,h)  dh  (8.5) 


where 


C  =  1.3^  X  10“'^/f^  mks  units  where  f  is  the 

ceurrier  frequency  of  the  ranging  signal  in  Hz. 

X  =  ionospheric  geographic  latitude  of  the  obser¬ 
vation  station  (in  fitting)  or  the  evaluation 
station  (in  predicting) 

0  =  ionospheric  geographic  lonrtldue  of  the  obser¬ 
vation  station  (in  fitting)  or  the  evaluation 
station  (in  predicting) 

t  =  UT  time 

R  =  effective  value  of  R^^ 

N^v(X,0,t,R)  =  maximum  electron  density  of  the  F2  layer  at 
X,0,t,  and  R  ,  given  explicitly  as  a 
function  of  R 

W(X,0,t,h)  =  normalized  electron  density  profile  at 
X,0,  and  t;  independent  of  R 

Strictly  speaking,  the  normalized  vertical  larofile,  W  , 
should  be  a  function  of  R  .  The  independence  of  R  is 
introduced  here  only  as  a  simplification,  as  discussed  above. 


so  that  the  value  of  J*  W(X,0,t,h)dh  may  be  obtained  from 
the  Algorithm  I  results.  The  value  of  J  W(X,0,t,h)dt  is 

(8.6) 


then  given  as 

T  (X,0,t,R,,>  ) 


where  R 


12^a 


is  the  value  of  R22  used  in  making  the 


Algorithm  I  predictions. 
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given 


vhere 
taken 

predictions  of  the  critical  frequency  of  the  F- layer,  fc?8 
[Jones  and  Obitts.  1970]* 

The  fitting  function  Eq.  8.7  was  then  least  sqviares 
fitted  to  the  experimental  t  (i.e.,  j  )  values  obtained  at 

the  observation  stations.  The  fitted  value  of  R  thus  ob¬ 
tained,  when  used  in  Eq.  (8.7)  at  the  evaluation  stations, 
gives  the  updated  prediction  at  the  evaluation  station. 

The  fitting  span  in  this  algorithm,  unlike  the 
case  in  Algorithm  II,  has  no  defined  minimum  length.  However, 
to  provide  a  meaningful  estimate  of  R  ,  it  was  felt  that  the 
fitting  span  should  be  at  least  2k  US  hours.  As  in  Algorithm 
II  this  was  taken  as  the  span  length  so  as  to  make  the  model 
responsive  to  daily  variations.  Longer  fitting  spans  would 
make  the  model  more  sluggish  in  its  response.  The  time  betvreen 
the  end  of  the  f-'ttlng  span  and  the  beginning  of  the  prediction 
time  was,  as  in  Algorithm  II,  taken  as  l/2  hour*  A  prediction 
span  of  3  hours,  as  in  Algorithm  H,  was  found  to  be  suitable. 
Shorter  spans  did  not  show  any  improvement  in  the  results. 


From  Eqs.  (8.9)  and  (8.6)  the  fitting  function  is 
as 

r  T„(A.,©,t,R,« )  ^ 

(8.Y) 

I  max  (Lfe^a  | 

R  is  the  fitting  parameter  and  the  functional  form 
for  Nmnv(X,9,t,R)  is  that  utilized  by  ITS  in  their 
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iMBliHiBlMM 


%  (t)  =  (Of  C  J*  w(h)  dh}  K(D)  cos  X 
cl 

{P  C  j*  ■w(h)  dhj  K(D)u)  cos  X^cos  ^gSln  tB(t-12) 

+  {V  C  J  w(h)  dh}  K(D)  (8.2a) 

atvi 

C  J  g(h)  dh) 

,  .  -6(t-t  ) 

+  {N^  C  J  g(h)  dh]  (1-e  )  (8.2b) 

Itt^oslblon  of  the  continuity  criteria 


T  (t  ) 


restricts 


%  <‘s.) 

n 


“  I  e<'')  “  ^  \  (t..) 

d 


In  the  study,  updating  was  effected  using  only 
vertical  time  delay  data,  Ty  .  In  this  case  there  is  no 
^^d  define  the  profiles  w(h)  and  g(h)  . 


a'  =  a  c  J  w(h)  dh 

P'  =  P  C  J*  w(h)  dh 

v'  =  Y  C  J*  w(h)  dh 

=  “«  C  J  g(h)  dh 

then  E(is.  (8,2)  yield 

T  (t)  =  a'  K(D)cos  X  +  P 

^d 


(8.3) 


K(D)  u)  cos  X^cos  XgSin  u)(t«»12) 


+  y'  K(D) 
for  the  daytime 


(8.4a) 
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and 


T  (t)  = 


-6(t-t  ) 

(KJ  +  K  (l-e  )  (8.Ub) 


for  the  nighttime. 


Thus  a  fit  ifc  made  directly  for  Of',  p',  and 
N  without  explicit  definition  of  f(h)  and  g(h)  .  In 
the  implementation  of  the  fitting  process,  the  daytime 
function  Eq.  (8.4a)  would  be  fit  first  so  r.s  to  define  t  (t  ) 
which  is  required  for  the  nighttime  function  given  by  ^d 
Eq.  (8.4b).  The  fitting  functions  t  (t)  and  T  (t)  as 

^d  \ 

given  by  Eqs.  (8.4)  are  linear  in  the  fitting  parameters 
^  t  ^  t  y  t  and  .  Thus  it  is  possible  to  fit  thir  model 
directly  to  the  residuals  of  the  observation  station  ditn 
obtained  in  Algorithm  I. 

The  minimum  length  of  the  fitting  time  span  is 
dictated  by  the  model  itself.  The  model  is  a  function  of 
local  time.  Thus,  in  order  to  obtain  both  day  and  ni^t 
estimates  at  any  observation  station,  allowing  for  the 
varying  longitudes  of  the  stations,  the  minimum  fitting  time 
span  is  24  hours  or  less.  The  length  of  the  span  that  was 
chosen  for  the  study  was  taken  as  24  hours  to  be  consistent 
with  our  attempt  to  track  the  daily  variations  expected  in 
the  residuals.  A  suitable  irediction  span  (that  is,  the 
time  jpan  over  which  the  fitted  parameters  are  used  for  up¬ 
dating)  for  purposes  of  making  the  model  responsive  to 
fluctuations  in  the  residuals  was  found  to  be  3  hours.  The 
time  between  the  end  of  the  fitting  span  and  the  beginning 
of  the  prediction  span  was  taken  to  be  l/2  hour.  In  an 
operational  system  this  would  correspond  to  the  time  required 
for  a  central  ground  facility  to 

1)  gather  data  from  the  system  ground  stations 

2)  perform  the  canputations 

3)  inject  corrections  into  the  satellite  memory 


Algorithm  III  is  an  alternative  approach  to  a  real¬ 
time  updating  scheme.  The  fundamentals  of  the  algorithm  are 
given  in  Section  6.6. 
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EVAIMTION  BESUliTS 


INTRODUCTION 


In  Section  8.1,  the  general  outline  of  the  evalua¬ 
tion  of  the  algorithms  in  terms  of  evaluation  conditions 
was  discussed  and  the  particulars  involved  in  the  In^lemen- 
tations  of  each  of  the  algorithms  were  given  in  Sections 
8.5.1j  8.5.2,  8.5.3.  The  basis  of  the  evaluation  of  each 
algorithm  was  the  difference  in  nanoseconds  referenced  to 
1600  MIIz  between  the  predicted  vertical  group  time  delay  and 
the  group  time  delay  as  inferred  frcmi  the  experimental  Fara¬ 
day  rotation  data,  .  This  difference,  denoted  as  6t  , 
is  termed  the  residual?  ^ 

Due  to  constraints  imposed  on  the  study  and  the 
expense  of  numerical  Integrations  inherent  to  the  algorithms, 
it  was  necessary  to  pick  frcan  each  evaluation  condition  a 
representative  data  rate  and  subset  of  days  on  which  to  base 
the  evaluation.  The  data  rate  was  taken  as  one  data  point 
per  20  minute  interval.  The  representative  subset  of  days 
for  each  evaluation  condition  is  called  the  subsetted  evalu¬ 
ation  condition  and  is  given  in  Table  9.1*  In  future  refer¬ 
ences,  the  term  evaluation  condition  will  mean  the  subsetted 
evaluation  condition. 

The  results  p:'esented  for  each  of  the  evaluation 
conditions  to  the  extent  permitted  by  the  Faraday  rotation 
data  are: 


Plots  of  the  monthly  rms  of 


6t^  (for 


each  algorithm)  for  each  one  hour  (UT)  interval, 
for  each  station  (both  evaluation  and  observa¬ 
tion),  for  each  month  of  the  evaluation  condi¬ 
tion,  together  with  the  rms  value  of  and 

e 

6t,^  over  the  month  and  the  number  of  residuals 
used  in  computing  these  values  (Appendices  A 
and  B). 

Plots  of  the  cumulative  frequency  distributions 

of  the  Ty  and  6t„  (for  each  algorithm)  over 
e  » 

the  entire  span  of  months  of  the  evalxiaticxi 
condition  together  with  the  rms  value  of  Ty 
and  fiTy  over  the  entire  span  and  the  ® 
number  of  residuals  used  in  conq^nzting  these 
values  (Appendix  C). 
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Evaluation  Condition 


CONDITIONS 


Subset  of  Data  Used  In  Studs 


Day  1  to  91  -  1968 
(January,  February,  March) 


Day  1  to  91  -  1968 
(Janiwry,  February,  March) 


Days  1  to  31,  91  to  220, 

182  to  212  -  1965 

(January,  April,  July) 

Days  1  to  3I,  92  to  121, 

183  to  213,  275  to  305  -  1968 
(January,  April,  July,  October) 

Days  1  to  31,  91  to  120, 

182  to  222,  274  to  30Jf  -  1969 . 
(January,  April,  July,  October) 


3)  Tables  of  the  values  of  the  correlation  coeffi¬ 
cient  of  the  6t  residuals  (for  each  algorithm) 
for  each  of  the  station  pairs  designated  in 
Table  9.2,  for  each  hour  (IB?)  Interval  over  the 
month  of  the  evaluation  condition,  together  with 
the  number  of  residuals  used  in  commting  each 
correlation  coefficient  (Appendix  D). 

4)  Tables  of  the  values  of  the  correlation  coeffi¬ 
cient  of  the  6t  residuals  (for  each  algorithm) 
for  each  of  the^station  pairs  designated  in 
Table  9.2  over  the  entire  span  of  months  of  the 
evaluation  condition  and  the  number  of  residuals 
used  in  computing  these  values  (Appendix  E) . 

5)  Tables  of  the  cumulative  frequency  distribution 
of  the  daily  correlation  coefficient  of  the  6t 
residuals  (for  each  algorithm)  for  each  of  the  ^ 
station  pairs  designated  in  Table  9.2  over  the 
entire  span  of  months  of  the  evaluation  condition 
and  the  number  of  residuals  used  (Appendix  F). 

On  each  of  the  cumulative  frequency  plots,  Appendix  C, 
a  table  has  been  added  which  summarises  the  plot  by  giving 
percentage  of  and  6t  values  which  fall  in  the  four 

intervals,  ±  3  ns,  ±  9  ns,  ±  18  ns,  and  ±  27  ns.  (Note  that 
1  ns  “•  1  foot.)  The  plots  are  given  on  a  probability  grid 
with  the  ordinate  range  0.01  percent  to  99*99  percent  and 
the  abscissa  range  -70  ns  to  +110  ns  in  1  ns  steps.  To 
determine  0  percent  subtract  1  ns  from  the  left-most  abscissa 
point  of  the  curve  and  to  determine  100  percent  add  1  ns  to 
the  right-most  abscissa  point  of  the  curve.  The  residual 
pairs  used  in  computing  the  correlation  coefficients  were 
sepeurated  by  a  time  interval  of  not  more  than  19  minutes. 

A  given  residual  was  not  used  in  more  than  one  residual  pair 
for  a  given  correlation  coefficient  evaluation  pair. 

9.2  DISCUSSION 

®ie  results  of  the  study  are  given  in  detail  in 
the  figures  and  tables  located  in  the  appendices.  In  this 
sectiem  we  wish  to  bring  out  sane  of  the  salient  features 
of  the  results. 

Table  9,3  summarizes  the  rms  of  the  6t  residuals 
over  evaluation  conditicxis  for  the  evaluation  stations. 


Table  9.4  summarizes  the  rms  of  the 


residuals  over 


TABI£  9*2 


C(»BEMriON  COEFFIClEm;  EVALUmON  PAIRS 


Evaluation  Condition  Station  Pair 

Stanford  ATS3  -  Stanford  AUKl 

1  Stanford  ATS3  -  Urbana  ArS3 

Ibrbana  ATS3  -  Sagamore  Hill  ATS3 

2  Stanford  ATS3  -  Stanford  ATSl 

Stanford  AlTSl  -  Stanford  ATS3 

Stanford  ATSl  -  Clark  Lake  ATSl 
Stanford  ATSl  -  Port  Collins  ATSl 
Rosman  ATS3  -  Urbana  ATS3 
Areclbo  ATS3  -  CoM  Bay  ATSl 


5 


S 1 


St-inford  AISl  1 
Stanford  ATS3  1 
Urb  ina  ATS3  1 
S>ii;rimore  Hill  ATS3  1 
Stanford  AT31  2 
Stanford  AT33  2 
Urbana  AT33  2 
Honolulu  syiC0H3  3 
Sagamore  Hill  AJS3  U 
Honolulu  AT31  It 
Stanford  ATSl  ^ 
Stanford  ATS3  $ 
Clark  I/ike  AI31  9 
Cold  Bay  ATSl  5 
Fort  Collins  ATSl  $ 
Urbana  ATS3  5 
Rosmnn  AT33  5 
Arecibo  AT33  ^ 


T  (ns) 
e 


TAW£  9.3 

suHURY  OF  me  VAUEs  am  bvawation 

CCHOmONS  FOR  EVAUIATIOH  STATIONS 
STy  Alg.  l(na)  8t^  Alg.  ll(ns)  Alg.  m(ns) 

(ST^Ty^) 


Hiadber  of 
Residuals  Used  In 
Computation  of 

roe _ 


ALL 

ALL  I'j.ce. 

<1.60 

•'1.27 

it.3« 

97i69 

(O.W) 

(0.27) 

(0.2a) 

^  I 

~X  a 

s 


67 


TAB1£  9.h 


SUWUBY  OF  R>6  VAIASS  OVER 
EVAUUATKm  CONDlTlOieS  FOR  OBSERVATION  STATIOtS 


.'Stntion 

E.C. 

Ty^(ns) 

Alg.  l(na) 

ST^  Alg.  ll(ns) 

e 

6t^  Alg.  lll(nB) 

Niiid)er  of 

Hoslduals  Uaed  In 
Computation  of 

RW 

Areclbo  ATS3 

1 

19.35 

8.13 

(0.1.2) 

'..53 

(0.23) 

If. 89 
(0.25) 

5180 

Honolulu  AT31 

1 

28.13 

U.73 

(o.w) 

7.88 

(0.28) 

10.70 

(0.38) 

5632 

Sagamore  Kill  AXS3 

2 

15.30 

7.12 

(0.».7) 

1..01 

(0.26) 

h.hj 

(0.29) 

5607 

Areclbo  AT83 

2 

19.<*8 

8.1'. 

(0.U2) 

I..I7 

(0.21) 

i..75^ 

(0.21.) 

529't 

Honolulu  ATSl 

2 

28. 3<* 

11.81 

(0.1.2) 

7.79 

(0.27) 

ss 

5717 

Stanford  SnC0H3 

3 

^.CZ 

l.«*7 

(0.26) 

1.32 

(0.23) 

1.1.3 

(0.25) 

4401 

Stanford  ATSl 

U 

lb.06 

U.79^ 

(0.3'.) 

2.59 

(0.18) 

3.11 

(0.22) 

6305 

Sagamore  Hill  AIS3 

5 

12.99 

6.18 

(0.U8) 

1..01 

(0.31) 

4.57^ 

(0.35) 

8538 

Edmonton  ATSl 

5 

13.90 

'..70, 

(0.55) 

5.10 

(0.37) 

5.06 

(0.36) 

7582 

Honolulu  ATSl 

5 

2)t.l<l 

9.78, 

(O.bO) 

5.33. 

(0.22) 

7.68 

(0.31) 

8236 

AU. 

ALL 

19.53 

8.29 

(o.w) 

5.09 

(0.26) 

6.43 

(0.33) 

62492 

evaluation  conditions  for  the  observation  stations.  These 

tables  also  include  a  column  for  the  r--®  values  to 

e 

provide  a  basis  for  measuring  the  perfori.»^^ice  of  each  of  the 
algorithms.  The  is  the  experimentally  determined 

e 

vertical  time  delay  Induced  by  the  ionosphere,  llie 

is  the  vertical  time  delay  remaining  after  making  a  correction. 

Thus  6t  /t  ,  which  have  been  included  parenthetically  under 

V  Vg 

the  time  value  for  each  algorithm,  give  a  measure  of  the 
performance  of  the  algorithm. 

In  general,  as  indicated  in  Table  9.3,  the  pre¬ 
diction  abilities  of  Algorithms  II  and  III  are  about  the 
same,  the  relative  prediction  ability  varying  from  case  tu 
case.  In  the  case  of  Honolulu,  Algorithm  II  is  superior  to 
Algorithm  IH  as  can  be  seen  in  both  Tables  9.3  and  In 

Evaluation  Condition  3,  a  period  of  low  solar  activity,  the 
difference  is  negligible. 

Summarized,  over  all  evaluation  stations  and  over 
all  evaluation  conditions,  the  rms  value  of  T  is  15,68  ns. 
The  rms  of  the  iresiduals  6t  for  Algorithms  If  II,  and 
III,  are  respectively  6.6O  ns,  4.27  ns,  and  4.38  ns.  Using 
these  rms  values,  the,  ratio  fof  each  algorithm  is 

respectively  0.42,  0,27,  0,28.  These  results  indicate  that 
at  the  present  time  long-term  predictions  (3-6  months)  of 
the  rms  vertical  time  delay  can  be  made  to  within  42  percent 
and  real-time  predictions  to  within  27  or  28  percent.  As 
pointed  out  in  Section  8.3,  the  only  data  used  in  the  present 
study  for  updating  in  Algorithms  II  and  III  were  the  t 

^e 

values  at  the  stations  designated  as  observation  stations. 

Overall,  the  long-term  predictions  of  Algorithm  I 
tend  to  be  Irw  (i.e..  Algorithm  I  tended  to  under  predict 
T„  ).  ®ie  Algorithm  II  and  III  residuals  are  more  normally 
^e 

distributed  about  zero  than  are  the  Algorithm  I  residuals, 

A  more  detailed  picture  of  the  distribution  of  the  errors 
is  given  in  Appendix  C. 

The  extreme  values  of  T  and  fir  observed  in 

the  study  occxirred  at  Honolulu  in  Evaluation  Conditions  1 
and  2,  The  values  were  86  ns  in  ,  -54  ns  in 

Alg.  I,  -50  ns  in  6t^  Alg.  II,  and  -43  ns  in  Alg,  III. 


-69- 


10.  lUTURE  EFFORTS 


10.1  mTRQDUCnON 

The  study  described  above  represents  a  first 
attempt  to  synthesize  a  prediction  technique  to  estimate 
the  lonosjherlc  induced  group  time  delay.  Because  of  this, 

Its  scope  had  to  be  constrained  so  that  several  ioqportant 
areas  of  interest  couM  not  be  explored.  The  evaluation 
of  the  two  real-time  algorithms  described  in  Section  8 
indicated  that,  while  both  were  acceptable  ftom  a  statis¬ 
tical  viewpoint,  periods  of  time  did  exist  when  further 
improvements  would  be  necessary  in  order  to  meet  the  navi¬ 
gation  accuracy  requirements.  In  addition,  several  funda¬ 
mental  questions  have  arisen  during  the  development  of  the 
algorithms  that  need  be  answered  to  substantiate  the  vali¬ 
dity  of  the  prediction  methods.  Consequently,  it  is  suggest¬ 
ed  that  extension  or  continuance  of  this  study  be  directed 
to  the  following  items. 

10.2  ADDinOML  TEST  FOR  REGIONAL  PREDICTIONS 

In  any  future  studies,  efforts  should  be  made  to  ob¬ 
tain  an  expanded  data  base  giving  adequate  coverage  in  longi¬ 
tude,  latitude  and  gecanagnetic  dip  along  with  solar  activity. 

In  p^icular,  it  would  be  useful  to  obtain  data  at  a  given 
location  from  two  or  more  satellites  so  as  to  obtain  a  aata 
base  which  repres'inted  a  navigation  situation.  Due  tr.  in.ierent 
limitations  of  Faraday  rotation  data  in  providing  the  cxp-p-ded 
geographic  coverage  Just  mentioned,  other  sources  of  c-da  will 
be  required.  One  such  source  presently  available  would 
U.S.  Navy  Navigation  Satellite  ^tem.  Ionospheric  reiractlon 
data  which  is  presently  available  as  a  by-product  of  both  the 
tracking  and  navigation  process  could  be  used  to  help  establish 
a  base  to  permit  a  world-wide  evaluation.  The  only  limitation 
of  the  data  would  be  that  it  would  provide  only  a  measia^  of 
the  ionospheric  time  delay  to  1000  km  (the  nominal  altitude  of 
a  satellite  in  the  system).  However,  this  would  measure  about 
80  percent  of  the  time  delay  to  synchronous  altitude  satellites. 
One  advantage  of  the  ionospheric  data  obtained  fron  the  U.S. 
Na%'y  Navigation  Satellite  System  would  be  the  measurement  of 
the  ionospheric  induced  range  error  over  a  large  region  of  the 
sky  during  a  transit  of  the  satellite  over  a  station.  This 
would  provide  a  measure  of  the  correlation  of  the  ionospheric 
range  error  from  near-simultaneous  observations  which  is  the 
quantity  of  interest  because  of  the  hyperbolic  nature  of  the 
pseudo-ranging  system. 
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10.3  GLOBAL  IREDICTIONS 


It  would  be  Instructive  to  broaden  the  region  over 
which  the  evaluation  Is  made  to  determine  the  true  global 
capability  of  the  current  technology.  Each  of  the  algorithms 
proposed  In  the  study  has  Inherent  to  It  a  global  prediction 
capability.  A  priori  estimates  of  the  group  time  delay 
obtained  throu^  Algorithm.  I,  which  does  not  require  real¬ 
time  data,  can  be  made  for  anywhere  in  the  world  subject 
only  to  specification  of  the  predicted  solar  activity  and 
the  locations  of  the  navigator  and  navigation  satellite. 

The  data  base  on  which  the  a  prior  estimates  are  based  con¬ 
sists  of  data  from  all  p?rts  of  the  world  Including  such 
areas  as  Africa,  Europe,  Russia  and  mainland  China.  The 
real-time  predictions.  Algorithm  II  and  III,  can  also  be 
used  In  their  current  state  to  estimate  the  group  time  delay 
anywhere  in  the  world.  The  Information  required  to  effect 
a  real-time  prediction  is  the  position  of  the  navigator  and 
the  navigation  satellite  and  real-time  data  taken  anywhere 
in  the  world. 

It  is  litq)ortant  to  note  that  all  the  algorithms 
utilized  in  this  study  are  already  defined  in  a  global 
sense  so  that  the  only  data  that  would  be  required  is  that 
on  which  the  evaluation  is  to  be  based. 

10.4  IMEROVEMBNTS  IN  MODEL  lONOSIHERBS 

The  model  ionospheres  utilized  in  the  study  to 
represent  the  electron  density  as  a  function  of  space  and 
time  represent  only  a  first  attempt  as  discussed  in  Section  6. 
More  detailed  models  do  exist  which  could  be  adapted  for 
utilization  in  the  prediction  algorithms  [King  and  Ruster^ 
1971].  Improvements  in  the  prediction  accuracy  are  to  be 
expected  if  more  detailed  model  ionospheres  are  exploited. 

10.5  UTIIIZATION  OP  OTHER  DATA  IN  THE  KffiDICTIONS 

In  the  completed  study,  the  data  on  which  the 
real-time  predictions  were  based  were  solely  the  group  time 
delay  inferred  from  Faraday  data  measured  at  the  observation 
stations.  There  is  a  significant  number  of  other  types  of 
data  that  could  be  used  in  the  prediction  algorithms.  Some 
of  these  are:  geomagnetic  and  solar  indices,  satellite 
observations  such  as  the  UV  intensity  and  ionosonde  data. 
These  should  be  explored  for  their  effectiveness  in  the 
predictions  of  the  group  time  delay. 


-72- 


10.6  VALIBm  OF  GROUP  TIME  DELAYS  INIEBRED  HIOM  FARADAY 
ROTATION  DATA 

As  discussed  in  Section  4,  the  group  time  delays 
■were  inferred  from  Faraday  rotation  data,  IMs  data, 
because  of  the  l/r3  (where  r  is  geocentric  distance) 
fall-off  of  the  geomagnetic  field,  is  sensitive  to  the 
lower  ionosphere  and  insensitive  to  the  higher  ionosphere. 
Consequently,  group  time  delays  inferred  from  Faraday 
observations  could  be  in  significant  error  because  of  the 
uncertainties  in  the  upper  atmosphere.  It  is  important 
that  the  group  time  delays  as  inferred  from  Faraday  data  be 
compered  with  simultaneous  measurements  of  the  group  time 
delay  itself.  An  experimental  study  directed  to  this  end 
would  establish  a  confidence  in  the  Faraday  data  which  does 
not  currently  exist. 

10.7  BFIECT  OF  MEASUREMENT  FREQUENCY  IN  FARADAY  ROTATION 
OBSERVATIONS 

Ibe  Faraday  data  utilized  in  the  study  was 
obtained  from  satellite  frequencies  in  the  nei^borhood  of 
137  (VRF)  ,  Ihe  data  was  then  scaled  to  the  l600  MRz 
frequency  of  Interest  (L-band)  by  use  of  Eq,  (4.2),  While 
this  formula  is  appropriate  for  a  near-homogenous  medium 
there  is  some  question  as  to  its  applicability  in  a  medium 
with  sharp  gradients  in  the  electron  density  distribution 
[Quier.  I963]  and  [Kelso.  1964].  Experiments  should  be 
undertaken  to  confirm  the  validity  of  the  extrapolation  of 
the  data  Arom  the  VHF  band  to  the  L-band. 

10.8  SCINTILLATIONS 

Scintillations  are  characterized  by  amplitude 
and  phase  variation  in  radio  waves  that  pass  throu^  the 
ionosphere.  This  phenomenon  arises  from  irregularities  in 
the  electrcm  density  dlstrlbxxtion.  There  is  some  concern 
as  to  the  depth  of  the  scintillations  at  L-band  and  whether 
they  would  Impact  significantly  in  the  navigation  measure¬ 
ments.  Consequently,  a  study  should  be  made  to  investigate 
the  effect  of  scintillations. 
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10.9  '.ilLOCITY  OP  PROPAOmON 

^^ere  are  several  velocities  of  propagation 
associated  with  the  transmission  of  electromagnetic  energy; 
the  phase  velocity,  the  group  velocity  and  the  signal 
velocity.  In  the  measurement  of  limited  wave  motion  the 
velocity  of  interest  is  generally  the  signal  velocity  which 
is  the  velocity  with  which  the  main  part  of  the  wave  motion 
propagates  in  a  dispersive  medium.  The  signal  velocity  is 
practically  the  same  as  the  group  velocity  whenever  the 
wave  motion  proceeds  without  strong  absorption.  As  dis¬ 
cussed  in  Section  4,  the  assunqition  that  the  signal  velocity 
is  equal  to  the  group  velocity  is  inherent  to  this  study. 

This  assumption  should  be  studied  to  investigate  the  accuracy 
of  the  approximation  for  the  L-band  ranging  system  proposed 
for  the  High  Altitude  Navigation  Satellite  System. 
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Fig.  A.28  VERTICAL  TIME  DELAY  MID  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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Fig.  VERTICAL  TIME  DELAY  AND  RESIDUALS4IOURLY  RMS  OVER  MONTH 
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NANOSECONDS  tat  1600  MHz) 


EVALUATION  CONDITION  2 
JANUARY  im 
HONOLULU- ATS  1 
NUMBER  OF  POINTS  2003 


SYMBOL 


\  AL6.I 
ALG.II 
ALG.III 


MONTHLY  RMS 
(m) 

28.78 

12.51 

7.54 

11.49 


.  .  • 

**•*■!?■ 


*  2  4  !  •  •  4  i  *  ,  • 
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Fifl.  A.30  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  R..:SOVER  MONTH 


NANU^CONDS  hit  1600  MH2) 
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EVALUATION  CONDITION  2 
FEBRUARY  1968 
SAGAMORE -ATS 3 
NUMBER  OF  POINTS  1790 


SYMBOL 


ALG.1 
ALG.II 
Sf  ALG.III 


MONTHLY  RMS 
(nt) 

15.41 

7.38 

3.23 

443 


♦  ♦  *  *  ♦ 


♦  • 

*  ••  • 
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Fig.  A.31  VER7ICAL  TIME  DELAY  AND  RESIDUAL&HDURLY  RMS  OVER  MONTH 
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EVALUATION  CONDITION  2 
FEBRUARY  1968 
ARECIBO-ATS3 
NUMBER  OF  POINTS  1979 

SYMBOL  ITEM  MONTHLY  RMS 

ins) 

♦  \  19.19 

•  Sfv  ALG.I  8.16 

A  \  ALG.II  3.99 
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Fifl.  A.32  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 
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EVALUATION  CONDITION  3 
JANUARY  1965 
HONOLULU  -  SYNCOM  3 
NUMBER  OF  POINTS  1377 


SYMBOL 

ITEM 

MONTHLY  RMS 
(n>) 

♦ 

7.96 

• 

St^  ALG.I 

2.79 

A 

ALG.II 
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ALG.III 
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Fig  A.37  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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Ut  1600  MHz) 
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Fig.  A.42  VERTICAL  TIME  DELAY  AND  RESIOUAL&4(OURLY  RMS  OVER  MONTH 
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NANOSECONDS  tat  1600  MHz) 


EVALUATION  CONDITION  4 
JANUARY  1868 
SAGAMORE  HILL  -  ATS  3 
NUMBER  OE  POINTS  1220 


SYMBOL 


MONTHLY  RMS 
int) 

15.60 

730 


fif^  ALG.  I  730 

ALG.  II  4.50 

Sf  ALG.  Ill  a97 


♦  ♦  ♦ 

♦  ♦ 


•  •  • 
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■  ^  ♦  ♦  ♦  ♦  ♦ 
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UNIVERSAL  TIME  (howt) 
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Fig.  A.43  VERTICAL  TIME  DELAY  AND  RESIDUALSNOURLY  RMS  OVER  MONTH 
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NANOSaCONOS  lat  t600  MHz) 


Fig.  A.45  VERTICAL  TIME  DELAY  AND  RESIDUALS4IOURLY  RMS  OVER  MONTH 


Fig.  A.46  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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Fifl.  A.47  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 


-126- 


NANOSECONDS  (at  1600  MHz) 


UNIVERSAL  TIME  (hours) 


Fig.  A.49  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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EVALUATION  CONDITION  4 
OCTOBER  1968 
HONOLULU-ATS1 
NUMBER  OF  POINTS  1941 


SYMBOL 

ITEM 

MONTHLY 

im) 

♦ 

29A4 

• 

Sf,  ALQ.I 

13.18 

A 

ALG.II 

&66 

■ 

8r„  AL6.III 

10.88 

s  •  • 

^  ■  • 


A  A 


•  ♦ 


♦  1 
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Fig.  A.50  VERTICAL  TIME  DELAY  AND  RESIDUAL&HOURLY  RMS  OVER  MONTH 
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Fig.  A.51  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 


NANOSECONDS  (at  1600  MHz) 


Fig.  A.54  VERTICAL  TIME  DELAY  AND  REStDUALS4K)URLY  RMS  OVER  MONTH 
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Fig.  A.55  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 


EVALUATION  CONDITION  5 
JANUARY  1969 
COLO  BAY  -  ATS  1 
NUMBER  OF  POINTS  2105 


SYMBOL 

ITEM 

MONTHLY  RMS 
(nt) 

♦ 

9.17 

• 

6tv  ALG.1 
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A 

ALG.II 
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Fig.  A.57  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 
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EVALUATION  CONDITION  B 
JANUARY  1969 
URBANA-ATS3 
NUMBER  OF  POINTS  2137 
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Fig.  A.59  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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Fig.  A.60  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 
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Fig.  A.61  VERTICAL  TIME  DELAY  AND  RESIDUALS4IOURLY  RMS  OVER  MONTH 
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tet  1600  MHz) 


EVALUATION  CONDITION  S 

APRIL  1989 

STANFORD- ATS  1 

NUMBER  OF  POINTS  2015 
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ITEM 

MONTHLY  RMS 
(nt) 

♦ 
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• 
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A 

ALG.II 
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EVALUATION  CONDITION  5 
APRIL  1969 
COLO  BAY- ATS  1 
NUMBER  OF  POINTS  2035 


SYMBOL 

ITEM 

MONTHLY  RMS 
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• 

AL6.I 

7.38 

A 

ALG.  II 

S.32 

G 

ify,  ALG.  Ill 
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Fla.  A.65  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 
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Fig.  A.69  VERTICAL  TIME  DELAY  AND  RESIDUALS4I0URLY  RMS  OVER  MONTH 
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Flo.  A.74  VERTICAL  TIME  DELAY  RESIDUALS+IOURLY  RMS  OVER  MONTH 
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EVALUATION  CONDITION  5 
OCTOBER  1966 


FORT  COLLINS  -  ATS  1 

NUMBER  OF  POINTS  1816 

■ 

SYMBOL 

ITEM 

MONTHLY  RMS 
(nt) 

♦ 

’■v. 

15.00 

• 

ALG.I 

6.57 

A 

ALG.II 
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- 

■ 

ALG.  Ill 
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Fig.  A.77  VERTICAL  TIME  DELAY  AND  RESIDUALS4iOURLY  RMS  OVER  MONTH 
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Fig.  A.81  VERTICAL  TIME  DELAY  AND  RESIDUALS4(0URLY  RMS  OVER  MONTH 
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Fig.  A.84  VERTICAL  TIME  DELAY  AND  RESIDUALS410URLY  RMS  OVER  MONTH 
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Fig.  A.91  VERTICAL  TIME  DELAY  AND  RESIDUAL&HDURLY  RMS  OVER  MONTH 
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TABD3  B.l 

Number  of  Data  Uccd  in  each  UT  Hoxxr  Interval 


Evaltiation  Condition  1 
Stanford  -  ATSl 


Hour  Interval 

January  1968 

Number  of  Data 
February  3968 

March  3.968 

10,1) 

76 

86 

80 

[1,2) 

81 

87 

81 

[2,3) 

77 

87 

81 

13,4) 

78 

86 

81 

[4,5) 

79 

85 

80 

[5,6) 

80 

84 

81 

[6,7) 

77 

84 

78 

[7,8) 

79 

84 

75 

[8,9) 

78- 

82 

75 

19,10) 

73 

81 

76 

[10,11) 

79 

78 

68 

[11,32) 

79 

80 

72 

[32,13) 

79 

80 

72 

[13,14) 

74 

80' 

74 

[14,15) 

•79 

86 

76 

[15,16) 

81 

87 

76 

[16,17) 

80 

87 

79 

[17,18) 

79 

87 

77 

[18,19) 

79 

86 

78 

[19,20) 

84 

87 

78 

[20,21) 

8? 

86 

81 

[21,22) 

&i 

87 

80 

I22,2j,) 

7y 

'’■< 

i'.o 

l23,2Ji) 

75 

86 

81 
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TABIil!  B.2 

Nuinber  of  Data  Used  In  each  UT  Hour  Interval 


,0 


■  Evaluation  Ccidxtion  1 
Stanford  -  ATS3 


Hour  Interval 

% 

January  1968 

Number  of  Data 
February’ I968 

March  I968 

[0,1) 

82 

80 

84 

11,2) 

83 

83 

82 

[2,3) 

85 

82 

81 

[3,^) 

87 

82 

80 

14,5) 

86 

76 

77 

15,6) 

84 

■  73 

75 

16,7) 

84 

69 

75 

I7v^^ 

87 

71 

71 

ia,9) 

84- 

74 

7li 

[9,10) 

86 

72 

77 

[10,11) 

84 

70 

75 

[11,12) 

82 

75 

77 

[32,13) 

83 

67 

78 

[13,iM 

86 

78- 

78 

[14,15) 

•77 

rr 

84 

[15; 16) 

82 

76 

84 

[16,17) 

88 

81 

83 

[17,1a) 

90 

84 

81) 

[18,19) 

87 

80 

84 

[19..20) 

86 

79 

82 

[20,21) 

86 

73 

8r. 

(21,22) 

86 

76 

81 

122,23) 

1"; 

r 

[23,2t) 

86 

•rr 

•  » 

83 
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TABIE  B.3 

Number  of  Bata  Used  in  each  UT  Hour  Interval 


Evaluation  Condition  1 
Urhana  -  AT33 


Hour  Interval 

January  I968 

Number  of  Data 
February  I968 

March  I968 

10,1) 

74 

87 

48 

ll,2) 

74 

87 

54 

12,3) 

75 

87 

54 

l?,4) 

75 

87 

54 

75 

87 

54 

15,6) 

75 

87 

54 

[6,7) 

75 

81 

57 

[7,8) 

75 

81 

57 

[8,9) 

75 

61 

57 

[9>10) 

75 

81 

57 

[10,11) 

75 

81 

57 

[11,12) 

75 

81 

57 

[12,13) 

75 

79 

57 

[13,1.4) 

75 

78- 

57 

[14,15) 

•72 

78 

57 

[15,16) 

73 

81 

60 

[16,17) 

76 

84 

60 

[37,18) 

74 

Qh 

60 

[18,19) 

72 

87 

57 

[19,20) 

70 

86 

52 

r  20,21) 

69 

84 

53 

[21,22) 

71 

86 

51 

[ '  ''  ) 

'I'y 

.  # 

5i 

[23,24) 

1 

M 

1 

87 

51 

TABUS  B  .4  . 

Number  of  Data  Used  In  each  UT  Hour  Interval 


Evaluation  Condition  1 
Sagamore  Hill-ATS3 


I 


I 


Hour  Interval 

Janviary  I968 

Number  of  Data 
February  196O 

March  .1,968 

[0,1) 

.  '^5 

75 

76 

[1,2) 

72 

75 

75 

[2,3) 

75 

75 

75 

[3,0 

75 

75 

75 

t't,5) 

75 

75 

75 

[5,6) 

75 

75 

75 

[6,7) 

75 

75 

75 

17,8) 

75 

75 

75 

(8,9) 

75- 

75 

75 

[9,10) 

75 

75 

75 

[10,11) 

75 

75 

74 

[11,12) 

75 

75 

72 

[12,13) 

75 

75 

73 

[13,14) 

75, 

-  .. 

'rr 

[14,15) 

•75 

75 

78 

[15,16) 

75 

75 

78 

[16,17) 

74 

78 

[17,18) 

72 

75 

78 

[18,19) 

71 

75 

78 

[19,20) 

72 

75 

78 

[20,21) 

72 

75 

78 

[21,22) 

r2 

74 

78 

U^-,23) 

72 

72 

7c 

[23,24) 

’  72 
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72 

78 
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TABIB  B.5 

Number  of  Data  Used  in  each  UT  Hour  Interval 


Evaluation  Ccud.’tion  1 
Arecibo  -  A'x’S3 


Hour  Interval 

January  1968 

Kum’ •  of  Data 
Febr— ;  y  3.968 

March  3968 

[0,1) 

84 

-- 

54 

11,2) 

83 

82 

54 

[2,3) 

79 

82 

54 

[3,M 

77 

83 

53 

[4,5) 

74 

84 

54 

[5,6) 

75 

•  80 

54 

[6,7) 

75 

80 

52 

[7,8) 

76 

77 

5i* 

[8,9) 

78- 

81 

54 

[9,10) 

78 

77 

5.3 

[10,11) 

80 

75 

54 

[11,12) 

81 

78 

52 

[12,13) 

83 

76 

49 

[13,14) 

81 

80 

51 

[14,15) 

•87 

81 

51 

[15,16) 

84 

87 

51 

[16,17) 

87 

87 

51 

[17,18) 

85 

87 

48 

[18,19) 

89 

86 

47 

[ 19,20) 

85 

87 

48 

[20,21) 

83 

86 

48 

[21,22) 

85 

37 

51 

[r-'jK'i) 

87 

'’I 

50 

i23,2J|) 

87 

87 

51 
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TABUS  B.6 

Nunber  of  Data  Uued  lu  each  UT  Hour  Interval 


Evaluation  Condition  1 
Honolulu  -  ATSl 


Hour . Interval 


January  1968 


Number  of  Data 
February  1968 


March  3.968 
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TAB1£  B.8 

.  Number  of  Data  Used  In  each  UT  Hour  Interval 

.Evaluation  Condition  2 
Stanford  -  ATfl3 


Hovu:  Interval 


January  1968 


Number  of  Data 
February  1968 


March  X968 


I 


Evaluation  Condition  2 
Urbana  -  ATS3 


Hour  Interval 

January  1968 

Number  of  Data 
February  3968 

March  1968 

[0,1) 

74 

87 

51 

[1,2) 

74 

87 

57 

[2,3) 

75 

87 

57 

[3,4) 

75 

87 

57 

[4,5) 

75 

87 

57 

[5,6) 

75 

87 

57 

[6,7) 

75 

81 

60 

[7,8) 

75 

81 

60 

[8,9) 

75- 

81 

60 

[9,10) 

75 

81 

60 

[10,11) 

75 

81 

60 

[11,12) 

75 

81 

60 

[12,13) 

75 

79 

60 

[13,14) 

75 

78 

60 

[14,15) 

.72 

78 

) 

[15,16) 

73 

81 

63 

[16,17) 

76 

84 

63 

[17,18) 

74 

8J» 

63 

[18,19) 

72 

87 

60 

[19,20) 

70 

86 

53 

[20,21) 

69 

84 

51 

[21,22) 

71 

86 

51 

[22,23) 

75 

87 

51 

[23,24) 

77 

87  • 

53 
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TABI£  B.IO 

Number  of  Data  Used  in  each  UT  Hour  Interval 


. Evalmtion  Condition  2 
,  Sagamore  Hill  -  ATSS 


Hour  Interval 


% 

January  1968 

February  1968 

March  I968 

[0,1) 

.75 

75 

85 

[1,2) 

72 

75 

84 

[2,3) 

75 

75 

84 

[3,M 

75 

75 

84 

[4,5) 

75 

75 

84 

[5,6) 

75 

•  75 

84 

[6,7) 

75  • 

75 

84 

[7,8) 

75 

75 

84 

[8,9) 

75 

75 

84 

[9,10) 

75 

75 

84 

[10,11) 

75 

72 

83 

[11,12) 

75 

75 

81 

[12,13) 

• 

75 

75 

be 

[13,14) 

75 

75  .. 

86 

[14,15) 

•75 

75 

87 

[15,16) 

.'5 

75 

87 

[16,17) 

74 

75 

87 

[17,18) 

72 

75 

87 

[18,19) 

71 

75 

87 

[19,20) 

72 

75 

87 

[20,21) 

72 

75 

87 

[21,22) 

72 

74 

87 

[22, 2i) 

72 

72 

87 

[23,2J4) 

72 

72 

87 

A.  Ik.*  •ittS.!^''- 


TABUS  B.ll 

Nuinber  of  Data  Used’  In  each  UT  Hour  Interval 


Evaluation  Ccvidition  2 
Areclbo  -•ATS3 


Hour  Interval 


Number  of  Data 


«* 

January  I968 

February  I968 

March  I968 

10,1) 

84 

82 

57 

[1,2) 

83 

82 

57 

[2,3) 

79 

82 

57 

13, 1») 

77 

83 

56 

14,5) 

74 

84 

57 

[5,6) 

75 

80 

57 

[6,7) 

75 

80 

55 

[7,^^) 

76 

77 

57 

[8,9) 

78  ■ 

81 

57 

[9,10) 

78 

77 

56 

[10,11) 

80 

75 

60 

[11,12) 

81 

78 

58 

[12,13) 

83 

76 

55 

[13,14) 

81 

80 

57 

[14,15) 

87 

81 

57 

[15,16) 

84 

87 

57 

[16,17) 

87 

87 

57 

[17,18) 

85 

87 

54 

[18,19) 

89 

86 

53 

[19,20) 

85 

87 

54 

[20,21) 

83 

86 

5J1 

[21,22) 

85 

87 

57 

122,23) 

87 

t  r 

[23,21i) 

87 

87  J  ■ 

5T 

-lb2- 


~  I  nnaa 


p-OTBi.w>. 


TABIE  B.12 

Number  of  Data  Used  in  each  UT  Hour  Interval 

Evaluation  Ccndition  2 
Honolulu  -  Al'Sl 


Hour  Interval 


Number  of  Data 

January  1968  February  I968 


^farch  1968 


TABUS  B.13 


Number  of  Data  Used  in  each  UT  Hour  Interval 


Evaluation  Condition  3 
Honolulu-Symcom3 


Hour  Interval 

January  I965 

Number  of  Data 
April  1965 

July  1965 

[0,1) 

56 

54 

66 

[1,2) 

55 

52 

66 

[2,3) 

53 

52 

65 

[3,4) 

54 

54 

66 

U,5) 

53 

54 

66 

[5,6) 

55 

53 

66 

[6,7) 

56 

53 

65 

[7,8) 

59 

51 

66 

[8,9) 

•  59 

49 

66 

[9,10) 

59 

48 

66 

[10,11) 

59 

53 

66 

[11,12) 

60 

51 

66 

[32,13) 

60 

45 

66 

[13,14) 

60 

45  •> 

66 

[14,15) 

60 

47 

66 

[15,16) 

60 

51 

65 

[16,17) 

60 

51 

66 

[17,3.8) 

'6o 

51 

66 

[18,19) 

60 

51 

66 

[19,20) 

^  58 

51 

63 

[20,23.) 

57 

48 

63 

[21,22) 

56 

48 

62 

>’• 

64 

[23,24) 

fiJj 

49 

66 

-l8it. 


TABLE  B.U 


Hour  Interval 

% 

Number  of  Data  Ucc'i  in  each  UT  llovir  Interval  ' 

Evaluation  CoriiUion  3 
^  .  StMford-SyncomS 

*• 

Number  of  Data 
January  1965  A-pril  1965 

July  1965 

10,1) 

55 

63 

65 

U,2) 

60 

63 

65 

[2,3) 

60 

63 

66 

13, M 

60 

63 

64 

U.S) 

60 

63 

66 

[536) 

59 

63 

66 

16,7) 

60 

63 

66 

17,8) 

60 

63 

63 

(8,9) 

•60 

63 

63 

19,10) 

60 

63 

63 

[10,11) 

60 

63 

63 

[11,12) 

60 

63 

64 

[12,13) 

60 

54 

61 

[13,14) 

60 

49  •• 

•• 

■58 

[14,15) 

60 

53 

64 

[15,16) 

60 

54 

66 

[16,17) 

59 

5T 

65 

[17,18) 

57 

62 

66 

[18,19) 

54 

63 

65 

[19,20) 

59 

63 

63 

[20,21) 

57 

63 

62 

[21,22)  • 

57 

62  ■ 

63 

[22,23) 

57 

60 

65 

[23,2'i) 

53 

-185- 

60 

66 

/f 


TABtS  B.15 

Number  of  Data  Used  in  each  UT  Hour  Interval 


Evaluation  Condit; on  h 
Sag^ore  Hill  -  ATS3 


Hour  Interval 


Nambcr  of  Data 
January  3968  April  I968 _ July  3.968 


October  I968 


Co,i) 

51 

54 

63 

78 

ll,2) 

48 

57 

62 

78 

12,3) 

51 

57 

60 

78 

l3,n) 

51  • 

57 

.  63 

78 

U,5) 

51 

57 

63 

78 

15,6) 

51 

57 

63 

78 

[6,7) 

51 

57 

63 

78 

17,8) 

51 

57 

63 

78 

18,9) 

51 

57 

63 

78 

l9,10) 

51 

57 

...  60 

78 

[10,11) 

51 

57 

63 

78 

[11,12) 

51 

«  • 

56 

63 

78 

[32,13) 

51 

.  57 

63 

78 

[13,11^) 

51 

57 

■63 

71 

[14,15) 

51 

57 

63 

73 

[15,16) 

51 

57 

63 

75 

[16,17) 

51 

57 

63  • 

76 

[17,18) 

51 

57 

63 

78 

[  18,19) 

50 

57 

64 

78 

[19,20) 

.  :.51 

57 

66 

78 

[20,21) 

51 

57 

66 

78 

[21,2?) 

51 

57 

66 

78 

[22,23) 

51 

57 

66 

78 

[23,2li) 

51 

57 

CG 

78 

-186- 

TABI£  B.16 


Number  of  Data  Used  in  each  UT  Hour  Intejrval 


Evalviatlon  Conditaon  4 
Honolulu  -  ATSl 


Hour  Interval  Nanbor  of  Data 


January  I968 

April  I960 

July  1968 

October  -1968 

10,1) 

83 

56 

66 

80 

11,2) 

61 

57 

66 

81 

.  12,3) 

62 

57 

66 

81 

I3,h) 

62  • 

57 

66 

81 

*14,5) 

60 

56 

66 

81 

15,6) 

60 

42 

66 

81 

16,V) 

59 

38 

■  66 

81 

* 

[7,8) 

59 

34 

66 

81 

[8,9) 

59 

27 

66 

81 

[9,10) 

57 

28 

...  66 

81 

[10,11) 

57 

33 

63 

81 

[11,12) 

57 

•  • 

39 

65 

81 

112,13) 

57 

.  34 

65 

81 

[13,14) 

57 

39 

•66 

81 

[14,15) 

55 

39 

66  ■ 

81 

.  [15,16) 

54 

29 

65 

81 

[16,17) 

53 

16 

66  . 

80 

[17,18) 

52 

4l 

66 

01 

[18,19) 

56 

53 

66 

81 

[19,20) 

'60 

54 

66 

81 

[20,21) 

63 

56 

66 

81 

[23,22) 

54 

57 

66 

80 

[22,23) 

5^» 

57 

66 

81 

[23,24) 

57  . 

57 

66' 

81 

-187- 

TA1JIJ2  B.1T 

Kurnbcr  of  J)ata  Uccd  in  each  OT  lloiu*  Interval 

Evaluation  Condition  4 
Stanford  -  Al'Sl 

Houi*  Interval  N’liibcr  of  Data 

_  _ January  19^^  April  1<)CS _ Ju3.y  October  19'>8 


10,1) 

53 

57 

66 

80 

U,2) 

60 

56 

66 

80 

[2,3) 

60 

56 

65 

81 

< 

13, ij) 

60  ■ 

57 

.  66 

■k 

H,5) 

63 

57 

66 

81 

[5,6) 

62 

56 

66 

79 

[6,7) 

61 

57 

66 

81 

[7,8) 

59 

57 

66 

81 

[8,9) 

60 

57 

65 

81 

[9,10) 

57 

56 

..  65 

81 

[10,11) 

63 

56 

65 

81 

111,12) 

62 

57 

64 

81 

[12,13) 

61 

■  55 

66 

81 

[13,1H) 

56 

56 

■65. 

•• 

81 

[li»,15) 

61 

56 

66 

81 

[15,16) 

63 

57 

66 

81 

[16,17) 

61 

56 

66  . 

81 

[17,18) 

59 

5*< 

66 

81 

[18,19) 

62 

56 

66 

81 

[19,20) 

.63 

57 

65 

81 

(20,21) 

62 

56 

64 

81 

[21,22)  •  . 

60 

57 

66 

81 

[22,23) 

60 

56 

66 

Gl 

(23,2*1) 

5*^ 

5*1 

-188- 

66' 

81 

TABIB  B.18 

Iluniber  of  Data  Used  in  each  UT  Hour  Interval 

Evaluation  Condi^-ion  5 
,  Stanford  -  ATSi 

Hour  Interval  Nuniber  of  Data 

January  19^»9  April  1969  July  iq6< 


TAWJi  B.20 

Nvunlicr  of  Data  Uaed  in  each  UT  Hour  Interval 


Evaluation  Condition  5 
Clark  Lake  -  ATSl 


Hour  Interval 


January  I969 


Muraber  of  Data 
ATOlll9r'9  July  1Q5<; 


October  1.Q6Q 


Co,i) 

19 

87 

42 

42 

11,2) 

18 

85 

39 

42 

12,3) 

18 

83 

40 

42 

13, h) 

18  • 

75 

4l 

39 

14,5) 

18 

73 

41 

37 

15,6) 

18 

64 

39 

38 

16,7) 

18 

66 

30 

39 

17,8) 

18 

69 

30 

36 

18,9) 

21 

58 

21 

26 

19,10) 

21 

61 

21 

2'i 

110,11) 

21 

55 

22 

21 

Ul,.12) 

21 

58 

24 

23 

102,13) 

21 

56 

25 

29 

113,14) 

21 

60 

■■  39  • 

37 

114,15) 

21 

63 

32 

42 

115,16) 

■  21 

'  64 

35 

42 

116,17) 

21 

70 

40 

41 

117,18) 

18 

77 

44 

42 

118,19) 

20 

79 

4l 

J'O 

119,20) 

21 

8o-  • 

41 

45 

120,21) 

21 

77 

42 

45 

C 21,22) 

21 

80 

42 

45 

122,23) 

21 

82 

41 

40 

21 


1»2 


5!ABW3  B.21 


Nuniber  of  Data  Deed  in  each  UT  Hour  Interval 


Evaluation  Condition  ^ 
Cold  Bay  -  ATSl 


Hour  Interval 

January  1969 

Number  of  Data 
April  1969 

July  I960 

10,1) 

83 

85 

63 

11,2) 

87 

84 

66 

12,3) 

87 

83 

66 

13,4) 

87  • 

81 

66 

14,5) 

87 

81 

66 

15,6) 

37 

81 

66 

(6,7) 

81 

66 

17,8) 

90 

81 

66 

18,9) 

90 

81 

66 

19,10) 

90. 

.  66 

110,11) 

90 

84 

.  66 

111,12) 

90 

85 

66 

112,13) 

90 

87 

60 

[13,14) 

90 

87 

to 

114,15) 

90 

87 

to’ 

Il5,i6) 

■90 

87 

to 

(16,17) 

90 

87 

60 

(17,18) 

90 

87 

to 

(18,19) 

90 

87 

to 

( 19,20) 

\  00 
VO 

87 

59 

(2*0,21) 

85- 

8r 

57 

(21,22) 

81 

87 

•57 

(22,23) 

81 

37 

57 

(23,24) 

81 

87 

57* 

TADJJ3  B.22 

Huiiibcr  of  Data  Used  in  each  UT  Ilour 


Bvaluation  Condition  5 
Fort  Collins  -  ATSl 


Interval 


Hour  Interval 

% 

January  I969 

Numher 
April  1969 

of  Data 

July  1969 

October  1969 

10,1) 

86 

24 

66 

74 

tl,2) 

87 

24 

45 

72  . 

12,3) 

87 

24 

40 

72 

13,M 

87  • 

26 

54 

72 

U,5) 

87 

27 

68 

72 

15,6) 

87 

27 

75 

72 

[6,7) 

87 

27 

75 

72 

17,8) 

87 

27 

81 

75 

18,9) 

87 

27 

84 

75 

r<),3.p) 

87 

27 

.  84 

75 

110,11) 

87 

27 

84 

71 

[11,12) 

•  • 

27 

81 

69 

[12,13) 

87 

27 

77 

69 

[13,1*^) 

87 

27 

■■  ..  • 

72 

[14,15) 

87 

27 

76 

76 

[15,16) 

87 

24 

76 

81 

[16,17) 

87 

24 

79 

81 

[17,18) 

87 

24 

80 

81 

[18,19) 

84 

24 

77 

81 

[19,20) 

.S'? 

'24  • 

78' 

81 

[20,21) 

87 

21) 

75 

81 

[21,22) 

87 

24 

77 

79 

(22,23) 

87 

22 

81 

• 

8') 

[23,2>i) 

87 

21 

*81 

79 

A 

*8 

I 


i  M 

i4 

I  1 

3  M 


a  I 


Hour  Interval 


TABUS  B.23 

Nunibcr  of  Data  Uccd  in  each  UT  Hour  Interval 


Evaluation  Condition  5 
■  Urhana  -  ATS3 


January  19^*9 


Number  of  Data 
April  1969  Jul.y  196^ 


October  196^ 


10,1) 

90 

87 

90 

90 

[1,2) 

90 

87 

90 

90 

[2,3) 

90 

87 

90 

90 

t3,i») 

90  • 

87 

90 

90 

U,5) 

90 

87 

90 

90 

[5,6) 

90 

87 

90 

90 

[6,7) 

90 

87 

90 

90 

[7,8) 

90 

87 

90 

90 

[8,9) 

90 

87 

90 

90 

19,10) 

90 

87 

90 

90 

[10,11) 

90 

87 

90 

90 

[11,32) 

90 

87 

90 

90 

[32,13) 

90 

87 

90 

90 

[13,14) 

90 

87 

••  '90 

.  .  90. 

[14,15) 

90 

87 

90 

90 

[I5,i6) 

89 

87 

90 

87 

[16,17) 

86 

85 

90 

87 

[17,18) 

8lt 

84 

90 

87 

[18,19) 

86 

84 

90 

87 

[19,20) 

86 

84-  ■ 

90 

90 

[20,21) 

87 

87 

90 

90 

[21,22)  .  . 

89 

87 

90 

90 

122,23) 

90 

87 

90 

90 

[23,24)  " 

90 

87 

■90 

90 

-194-  _ 

m 


m 


H 


li 


a 


n 


I 


f 


TABUS  B.24 

MumBcr  of  Data  Uced  In  each  UT  Hour  Interval 


Evaluation  Condition  5 
Rosinan  -  A'i‘S3 


Hour  Interval 


January  I9C19 


Numher 
April  1969 


of  Data 
July  19ft 


Octotcr  196*^ 


[0,1) 

87 

39 

80 

32 

11,2) 

86 

39 

80 

27 

t2,3) 

86 

38 

80 

27 

[3,4) 

84  • 

36 

79 

32 

14,5) 

81 

38 

83 

29 

(5,6) 

80 

39 

83 

29 

16,7) 

85 

37 

80 

29 

(7,8) 

87 

37 

87 

30 

(8,9)- 

82 

•  38 

82 

29 

(9,10) 

85. 

36 

82 

30 

(10,11) 

87 

34 

88 

30 

(11,12) 

86 

38 

87 

32 

(32,13) 

90 

36 

87 

35 

ri3,l4) 

89 

36 

■■  k  . 

•  33 

(14,15) 

78 

36 

82 

33 

(15,16) 

75 

36 

82 

32 

(16,17) 

83 

36  • 

88 

36 

(17,18) 

83  *' 

34 

86 

36 

(18,19) 

86 

35 

..87 

36 

(19,20) 

33 

83 

32 

(80,21) 

90 

34 

82 

36 

(21,28)  . 

87 

32 

•79 

33 

(82,23) 

66 

36 

67 

• 

35 

(23,24) 

87 

36 

79 

35 

.¥  i 

a  I 

53 


-195- 


SABl^  B.26 


Nuiriber  of  Data  Used  in  each  DT  Hour  Interval 


Evalxiatlon  Condition  5 
Sagamore  Hill  -  ATS3 


Hour  Interval  Number  of  Data 

_ January  I969  April  1969  July  October  ‘ 


[0,1) 

90 

87 

90 

90 

[1,2) 

90 

87 

90 

90 

[2,3) 

89 

87 

90 

90 

13,4) 

90  • 

8'^ 

90 

90 

14,5) 

90 

87 

90 

90 

1 5,6) 

90 

’87  • 

90 

90 

16,7) 

87 

87 

90 

90 

17,8) 

90 

87 

90 

90 

18,9) 

90 

87 

90 

90 

19,10) 

90. 

87 

.  •  ■  90 

90 

[10,11) 

90 

87 

.  90 

90 

[11,12) 

90 

87 

90 

90 

[12,13) 

90 

87 

90 

90 

[13,14) 

90 

86 

••  90 

90 

[14,15) 

89 

87 

89 

90 

[15,i6) 

•  88 

87 

87 

90 

[16,17) 

87 

87 

90 

90 

[17,18) 

87 

87 

90 

90 

[18,19) 

87 

85 

90 

90 

[19,20) 

■■  87 

87- 

90 

90 

[20,21) 

.88 

87 

90 

89 

[21,22)  . 

89 

87 

90 

90 

[22,23) 

90 

87 

90 

90 

[  23,2*1 ) 

90 

87 

■96 

90 

197' 


Number  of  Data  Used  in  each  UT  Hour  Interval 


Evalu£ition  Condi'^iun  5 
H.onolulu  -  ATSl 


*  1 


Hour  Interval 


January  I969 


Number  of  Data 
April  1969  July  19^ 


Octobo'i 


10,1) 

81 

86 

90 

90 

[1,2) 

80 

87 

90 

90 

12,3) 

81 

85 

90 

90 

13,J0 

83  • 

a 

90 

90 

t‘o5) 

86 

81 

90 

89 

[5,6) 

87 

81 

90 

90 

[6,7) 

87 

81 

90 

90 

[7,8) 

87 

80 

87 

89 

[8,9) 

87 

80 

89 

89 

[5,10) 

87. 

79 

88 

89 

[10,11) 

87 

73 

.  85 

85 

[11,12) 

73 

83 

88 

[32,13) 

87 

76 

87 

88 

fl3,l4) 

87 

79 

••  k  .. 

89 

[lit,l5) 

87 

81 

86 

90 

[15,i6) 

■87 

79 

87 

89 

[i6-,17) 

87 

80 

81 

87 

[17,18) 

87 

78 

78 

87 

[18,19) 

87 

81 

81 

89 

[19,20) 

86 

a  • 

89 

90 

[20,21) 

84 

35 

90 

90 

[21,22) 

8^4 

87 

90 

90 

12;.', 2.0 

bh 

}SiJ 

89 

[23,2l|) 

84 

87 

88 

89 

*  I 
I  1 

^  A 


i  I 


-4  ? 

■4 


Pnceding  pigt  link 
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NANOSECONDS  (at  1600  MHz) 

Fki.  C.3  VERTICAL  TIME  DELAY  AND  RESIDUALS-CUMULATI VE  FREQUENCY  DISTRIBUTION 


eVALUATION  CONDITION  1 
ARECIBO-ATS3 


S  cs 

§u.z 

oa  -s 

I  -  “JS  # 


ui  ~ 

§0^1 
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= 
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giSSS 
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2  5s 
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CONDITION  2 


EVALUATION  CONDITION  2 
HONOLULU  -  ATS  1 
JANUARY,  FEBRUARY.  MARCH  1968 
NUMBER  CF  POINTS  5717 


Fig.  C.13  VERTICAL  TIME  DELAY  AND  RESIDUALS-CUMULATIVE  FREQUENCY  DISTRIBUTION 


EVALUATION  CONDITION  4 
SAGAMORE  HILL  -  ATS  3 


FIb.  C.18  VERTICAL  TIME  DELAY  AND  RESIDUALS-CUMULATIVE  FREQUENCY  DISTRIBUTION 
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TABI£  D.x 

Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  l 
January  1968 

Correlation  Pair:  Stanford  ATSl  -  Stanford  ATS3 


UP  Hour  Interval 


6Ty  Alg.  I 


6t^  Alg.  II 


6t^  Air,,  in 


1  Ot  1) 

0.88529 

0,98449 

C. 

924C5 

1  If  2) 

0.95<‘a0 

0.98271 

0. 

95806 

I  2f  3) 

0.95240 

O.vBi 17 

C. 

95277 

1  3«  4) 

0.95407 

0.98192 

0. 

96210 

1  4,  5» 

0.B95B9 

0-97265 

c. 

96^- 60 

i  5f  6) 

0.87755 

0.97935 

V  • 

98C54 

1  6f  7) 

0.86595 

0. 97945 

0. 

99146 

1  7,  B) 

0.30337 

0.972i8 

c. 

8'>ol  8 

1  8f  9) 

0.80144 

6.95738 

0. 

89144 

i  9fl0) 

0.768«6 

0.93321 

c. 

R6C92 

llOfll) 

0.79554 

0.95397 

0. 

f.72n2 

illf 12) 

0.72577 

0.92657 

c. 

B^r>24 

112,13) 

0.72629 

0.93257 

c. 

81735 

113,14) 

0.77943 

0,95753 

0.39214 

114,  15) 

0.59343 

0.99536 

C 

84^77 

113,15) 

0.93525 

0.71103 

c 

94379 

116,17) 

0.32353 

0.91705 

•  c 

938f.5 

117, IB) 

0.37617 

0.58189 

0 

94329 

113,19) 

0.93372 

0,90323 

r 

93957 

119,20) 

0.95581 

0.95652 

‘-62C4 

l?0,21) 

0.93210 

0. Q5542 

f 

■ri  7  L'  7 

121,22) 

0.90339 

0.94098 

C 

97312 

122,23) 

0.57370 

0.95569 

9:;C  37 

123,2-) 

0.79533 

0.95701 

c 

90”  9? 

Nominal  Valxie  in  Degrees  for  Month 

Stanford  ATSl 

Stanford  ATS3 

Elevation 

57.5 

38.C 

Azimuth 

221.5 

139.5 

Eff.  Latitude 

34.6 

34.6 

« 

Eff,  Longitude 

234.8 

L 

240,7 

-232. 


Nuniber  of 
Residual  Pairs 


71 

7^ 

74 

75 
75 
75 
75 
75 

74 
*•  > 

75 
73 

75 
70 

7? 

76 
79 

78 
PI 

79 

77 
7f> 
70 


5a 


TABIE  D.2 

i’orr*  lation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  1 
January  1968 

Correlation  Pair:  Stanford  ATS3  -  Urbana  ATS3 


OT  Hour  Interval  6t„  Alg.  I  6t..  Alg.  II 


1  Of 

1) 

1  1* 

2) 

1  2f 

3) 

1  3t 

4) 

1  4, 

5) 

I  5t 

6) 

1  6, 

7) 

1  7« 

8) 

1  fif 

91 

1  9, 

10) 

no, 

U) 

111, 

12) 

1  12, 

13) 

113, 

14) 

114, 

15) 

13  5, 

16) 

116, 

17) 

117, 

18) 

118, 

19) 

119, 

20) 

120, 

21 ) 

I2li 

2?) 

122, 

23) 

123, 

?4) 

0.65786 
0.85277 
0.86183 
0.76662 
0.71567 
0.77153 
0.64089 
0.64631 
0.5<.672 
0.35266 
0.10366 
0.06496 
0.40377 
-0.D3113 
0.34020 
0.14035 
0.51968 
0.70404 
0.79030 
0.70932 
0.70357 
0.620o5 
0.5  3  9  8  7 
0.53429 


0.96797 
0.83296 
0.88260 
0.85610 
0.33159 
3.3-;777 
0.37513 
0.90355 
0.86269 
0.66151 
0.61 l64 
0.60554 
0.84770 
0.33933 
0.17591 
0.17923 
0. 67560 
0.70374 
0.80406 
0.31915 
0.35736 
0.9762'. 
0.39979 
0*36105 


Number  of 

6t„  Alg.  Ill  Residual  Pairs 


C.71:b7 
C. 88168 
C.9369? 
C. 90914 
0.89165 
0. 9:202 

0.63178 
C.63C50 
r. 57142 
f. 33737 
<‘.27391 
C.  271  82 
C.6l*3€5 
0.52456 
C. 41551 
C. 19539 
0.44874 
0.77694 

C.7'»r92 
<■.63069 
0.57*'41 
0.45197 
0.40797 
0.4f 174 


nominal  Value  in  Degrees  for  Month 
Stanford  ATS3  Urbana  AT'S3 


Elevation 
Azimuth 
Eff.  Latitude 
Eff.  Longitude 


Stanford  ATS3 

38.0 

139.5 

34.6 

240.7 


U3,2 

190.2 

37.0 

271.1 


TABIE  D. 

Correlation  Coefl’icionts  of  Vertical  Time  Delay  Residvials 
Evaluation  Condition  1 
Jantiary  I968 

Correlation  Pair:  Urbana  AIIS3  -  Sagamore  Hill  ATS3 


Number  of 

UP  Hour  Interval  6t^  Alg.  I  6t^  Alg.  II  6t^  Alg.  Ill  Residual  Pairs 


1  0«  I  ) 

0.74112 

0.3959O 

0.73224 

71 

1  It  2) 

0.35433 

0.9197.. 

0.9C76? 

66 

1  2t  3) 

0.9034;. 

0.93597 

0.93346 

72 

1  3t  4) 

0.39478 

0.93069 

0.9291-: 

72 

1  4t  5) 

0.93672 

0.95397 

C. 95230 

72 

i  5t  6) 

0.89572 

0.92717 

0.93255 

72 

1  6t  7) 

O.B  7370 

0.91917 

0. 36619 

72 

1  7,  8) 

0.64025 

0.39935 

''.81998 

72 

1  B,  9) 

0.80731 

0.88079 

C.75'‘:6 

72 

1  9flO) 

0.82023 

0.75601 

C. 7051 7 

72 

1 lOt  11) 

0.610S7 

0.64939 

C.4395? 

72 

|11»12) 

0.55^62 

0.74448 

C.5372Q 

72 

|12tl3) 

0.24660 

0.49711 

0.56745 

7? 

1 13f 14) 

0.64433 

0.61122 

0.64512 

72 

|14tl3) 

0.70254 

0.59184 

0.6  09.1 5 

59 

|15t 16) 

0.85907 

0 • 54365 

r,  3C95r 

69 

116,17) 

0.77290 

0.72017 

.C.fc.51't5 

6° 

117,18) 

0.82333 

0.90t4? 

0,91546 

65 

1  18,  19) 

0.95107 

0.93390 

C.96C06 

62 

119,20) 

0.94  342 

0.94357 

C.95''74 

61 

|?0,?1 ) 

0.39211 

0.931 61 

C.  H6594 

60 

121  ,22) 

0.94354 

0.96513 

C.9?C04 

52 

122,23) 

0. 52726 

0.V0359 

C.7>9n 

66 

123,24) 

0.31-,30 

0. 83323 

0.  74!  j 

6? 

Nominal  Value  in  Degrees  for  Month 

Urbana  ATS3 

Sag^unorc  Hill  ATS3 

. 

Elevation 

I13.2 

:5.i 

Azimuth 

190.2 

213.3 

Eff.  Latitude 

37.0 

39.2 

• 

Eff.  Longitude 

271.1 

286.3 

^^^^SWBaBIWBMWE»iggSawmnu.-» 


TABIE  D.l^ 

Con  elation  Coefficients-  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  1 
Petruai^y  I968 

Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 


Ur  Hour  Interval 


6t^  Alg.  I 


6t^  Alg.  II 


6t  Alg.  Ill 

V 


Nunibe-r  of 
Residual  Pairs 


0, 

’  1 

C  69 

C.Q7313 

C.955''i5 

1, 

21 

0. 966?  6 

0.97C08 

C. 96969 

2 , 

31 

0.96.^  *57 

0.96582 

0.97241 

?» 

A) 

0,9?9^d 

C,Q42Ari 

C.9399? 

At 

f.) 

0.87526 

C. 91880 

0.92^  '••8 

5t 

0.9C89? 

0.92774 

6 1 

71 

O.fifGAA 

C. 93139 

0.90563 

•'t 

3) 

0.7?r  16 

0.93KB 

C.B^OOl 

“f 

9) 

0.70976 

r,04,Pft8 

<•  .351*56 

"•f 

101 

C .  7  3 .3  0  < 

0.'*.435B 

C .86162 

1 

i  -  f 

11  1 

C.A.fj?6 

C.022?" 

C. 76009 

lit 

1.?) 

O.A.:/.07 

0.«164e 

C.76O0L 

I2t 

1?1 

C.  3^*7  36 

C.9GA66 

0.63=61 

I2t 

141 

0.-331  ?<? 

0.91783 

C,63t  75 

!  A, 

151 

0.73517 

0.638-1 

!  5, 

161 

0.7*'6A7 

0.RR552 

0.34054 

■ft. 

171 

0.73?'^: 

0.P73C5 

0.685  32 

P.7,  ).«>) 
lift  l<-) 


.  ) 

,/.•? ) 

I 

1,7..) 


7  A ''.'9 

^,7 

G.  -Tl 

C.  i-ic‘5., 

^ 7?.-  •» 
<.<■ 


r.909CA 
0.r6164 
O.PB.?'»9 
G. 85913 
C.»=?776 
r.,  Q7A30 
r.«^2574 


(  .917'  3 

C.9Uf-5 

c,e.^:7: 

t.OlAjji 
f.  ,9-  2Ll 


Nominal  Value  in  Degrees  for  Month 


Elevation 
Azimuth 
Eff.  Latitude 
Eff.  Longitude 


Stanford  ATS3 


32.6 

329.7 

3'v.5 

241.9 


Stanford  ATSl 


TABIE  D.*; 


Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  1 
February  1968 

Correlation  Pair:  Stanford  ATS3  -  Urbana  ATS3 

Number  of 

WT  Hour  Interval  6t  Alg.  I  6t  Alg.  II  6t  ,  Alg.  Ill  Residual  Pairs 


I  0,  1) 
I  If  ?) 
I  2» 

I  3f  h) 
I  A,  51 
I  5.  b» 
I  6,  7) 
I  7, 

I  B, 

I 

IK,  11) 
lUfi?) 
P.2f  1?) 

IK,  K) 
llAfl*:) 
|l!5f  K) 
n.6,17) 
fl7,lf ) 
ll-’-fV^) 
119,7'^) 
l?.C,?l  ) 
|:if ?/) 

12?  ,  ’4} 


Elevation 
Azimuth 
Eff.  Latitude 
Eff.  Loncituie 


0.37531 

C. 851 AR 

C 

0.P9421 

f^.7?6A9 

c 

^  .BPS*'? 

C. 73827 

c 

O.9701A 

C. 80050 

c 

C.?.S75^ 

0.95092 

c 

r  .917?'; 

0.79147 

f 

C.76?3'» 

0.77AB2 

c 

C.h70«i: 

0.78681 

c 

0.B2919 

f 

'•.3A612 

C.6R677 

c 

'^.13776 

c 

0. ?57o4 

C.P1P21 

t 

r,  7ei.o3 

0.33518 

t 

1 

0.17982 

c 

C.A98  7? 

n 

C .47°n4 

0.72752 

r. 

0.7'-, 46? 

0.78513 

C 

C.51557 

0.P3573 

c 

C.8271A 

( 

C 

n.7UA7 

c 

C.  H''*)  jC 

C.f.op?'' 

f. 

O.T-lpf*'’ 

C..77C7C 

r 

t'.  ^9?.  11 

C. 87057 

r- 

C. 87511 

r 

Nominal  Value  in  Degrees  for  Month 

Stanford  ATS3 

Urbana  ATS3 

32.6 

43.6 

129.7 

176.  !• 

34.5 

37.0 

241.9 

272.0 

flpmi^  JAM JBJjl 


TABIE  D.<i 

Corralation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  l 
February  1968 

Correlation  Pair;  Ifrbana  AIS3  -  Sagamore  Hill  ATS3 


UT  Hour  Interval  6Ty  Alg.  I  6t^  Alg.  II 


Number  of 

6t^  Alg,  III  Residual  Pairs 


i 

i  1  Ct  1) 

0.81611 

C.71212 

C. 73514 

1 

1  1  If  ?) 

^■?.7qii! 

0.7555? 

0.74649 

1 

1  1  2. 

r.75«5?iO 

0.727P5 

0.70'  83 

1  1  1  •'t  h) 

C.732^6 

0,79169 

0.7^923 

e  1  1  A,  *?) 

G.6A257 

0.71771 

0,59947 

1  \  15,  (>) 

0.51648 

C. 63639 

C. 50296 

1  8  1  A.  71 

ri.4'“?7i 

C,62C51 

0.46148 

g 

1  1  7f  fl) 

C.35^ ’3 

C.63629 

0 • 35663 

■i. 

P  1  fit  '>1 

0.74061 

C.^2983 

t 

}  1  9  f  IC) 

C.5572A 

0.76^49 

C.6-:t57 

t 

1  ilOtlU 

0.455<»7 

0.78503 

C. 55174 

i  1 1 1  f  1 2  ) 

O.51''?0 

C. 65229 

C. 57669 

w 

i  ii2fi’) 

0.AA1C3 

C. 75338 

C. 5 3926 

1  Il2fl<^) 

C. 731 38 

0,83??*^ 

0,55411 

[  jlAfl?) 

C.82506 

0.8?97C 

C.7CU5 

il5fU) 

0.7«<'-<t? 

0.7411? 

C. 73354 

: 

lUtlTJ 

0.704^6 

0,69295 

0.75561 

4 

117f!f') 

C.6C22-5 

0.66230 

0.7825'^ 

i 

nsfi*?) 

0.646fl7 

0.69A76 

C.Pi5K 

1, 

0.8  8H7<, 

0.82118 

0  •  8  6“  3 1 

\ 

i?'*f  2M 

:.P«5’7 

''.87126 

0,86538 

|21f2?» 

0.‘-.2.?.76 

C. 890 Oh 

0- .  8  3243 

IZPf?.') 

o.v:  •'o? 

''.87766 

0,84535 

4 

I2?f,>) 

C. An?64 

0.7537A 

C.8->b'^3 

f 

Nominal  Value 

in  Degrees  for  Month 

Urbana  ATS3 

Sagamoi^  Hill  ATS3 

Elevation 

43.6 

38.4 

Azimuth 

176.5 

201.8 

Eff.  Latitude 

37.0 

39.3 

Eff,  liongitude 

272.0 

L - - - - 

287.4 

-237- 


# 

75 

1 

1 

! 

75 

i 

75 

i 

75 

1 

75 

1 

1 

75 

% 

-i 

69 

1 

-5^ 

69 

1 

69 

$ 

69 

% 

66 

j 

69 

•  .( 

ts 

s 

67 

j 

1 

66 

J 

! 

1 

66 

% 

69 

1 

72 

1 

72 

"5 

1 

75 

i 

s 

16 

74 

i 

72 

j 

-^1 

'1 

7? 

5 

1 

72 

f 

It 

72 

i 

TABIED.7 

Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  l 
.  March  1968 

Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 

Number  of 

OT  Hour  Interval  Alg.  I  6t^  Alg.  II  6t^  Alg'.  Ill  Residual  Pairs 


1  0, 

1) 

1  11 

2) 

1  2» 

31 

1  3t 

4) 

1  4, 

5) 

1  5t 

6) 

t  6« 

7) 

1  7, 

6) 

1  8t 

9) 

1  9, 

10) 

1  1C« 

11) 

ni» 

12) 

!12f 

13) 

n3f 

14) 

114. 

15) 

115. 

16) 

116, 

17) 

117, 

IP) 

P6, 

19) 

119, 

29) 

I2C, 

?l) 

l?li 

,22) 

1  2?  1 

,2?) 

123, 

»?4) 

0.93608 

C. 93661 

0.91735 

0.65795 

0, 68041 

0.36975 

0.43615 

0.46895 

C. 454  64 

0.33603 

C, 34561 

0. 38696 

C.C83C9 

0.14592 

C. 57637 

'>.72459 

n,85*’-7e 

0.89608 

0.9CC68 

r.37<59C 

0.83275 

0.91265 

0.91655 

9,90767 


0.98627 

0.98638 

0.98509 

0.96036 

0.93813 

0.93716 

0.93522 

0.94898 

0.95616 

C. 92551 

0.92C39 

0.91891 

0.88C09 

0.83217 

0.42733 

0.75371 

0.66170 

0.92415 

0.93515 

0.93627 

0.90928 

0.96289 

0.95615 

0.96756 


C. 93199 
C. 93192 
C. 91625 
C. 84358 
C. 76336 
C. 76447 
r.68318 
0.63481 
C.616‘:4 
C. 53380 
0.63247 
(  .6C589 
C.561C5 
C. 62639 
C. 64163 
C. 71479 
0.61981 
G. 89245 
‘;.90127 
C.8SC77 
8t;494 
C.«l<-44 
r, 93184 
C. 93575 


80 

79 

78 

77 
73 
7? 
69 

65 

66 
69 
61 
66 
66 
68 
76 

76 

78 

77 

78 
75 
77 
77 
8^ 
86 


Nominal  Value  in  Deerces  for  Month 


Stanford  AT53 

Stanford  ATSl 

Elevation 

31.3 

37  7 

Azimuth 

127.8 

221.2 

Eff.  Latitude 

34.5 

34.6 

Eff.  Longitude 

242.2 

_ 

234.8 

-238- 


R- 


TABIE  D.fi 

Corralation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  3 
March  I968 

Correlation  Pair;  Stanford  ATS3  -  Urhana  ATS3 


OT  Hoiar  Interval 


6Ty  Alg.  I 


6T^  Alg.  II 


6t^  Alg.  Ill 


i: 

'i 

lie 

5. 

i.-: 

?• 

111 

112 

1 

113 

Il« 

115 

t 

lu 

1 

117 

[f 

118 

139 

> 

I2C 

1 

121 

1 

12? 

U 

123 

0.P9340 

0.81515 

0.735B8 

0.77331 

0.92?92 

0,84487 

0,80044 

0,83545 

0.88125 

0.82848 

0,85184 

0.59614 

-0.32500 

-0,13853 

0.44943 

0.62388 
0.66374 
0.71657 
0.73901 
0,76166 
0.8''216 
0.85535 
C. 86292 
0.87C13 


nominal  Value  in  Degrees  for  Month 


l> 

C. 34541 

2) 

0.32413 

3) 

C.39C00 

4) 

0.4?5t8 

5) 

C. 31085 

6) 

0.08Q59 

7) 

-0.07281 

8) 

-0.14806 

9) 

0.15C78 

ICl 

0.22214 

111 

0.2CC82 

121 

-0,03879 

131 

-0.25293 

14) 

-0.19550 

15) 

0.23805 

16) 

0.48169 

17) 

C.5453C 

18) 

0,58787 

19) 

O.6057I 

20) 

0.52549 

?l) 

0,60593 

22) 

0.54394 

23) 

0.64642 

24) 

0.58190 
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C. 27344 
0.12522 
C. 20349 
C.19C53 
f ,25690 
0.19811 
C. 11574 
C.  0  5310 

0.39304 
0.60475 
C. 58462 
0,28925 
0.03869 
0.C7129 
16CP8 
34929 
4C794 
0,-9628 
39951 
27781 
32607 
51842 

c. 54554 

C.45529 


0. 
0  . 
0. 


C. 

C. 

C. 

0:. 


Stanford  ATS3 

Urbana  ATS3 

Elevation 

31.3 

43..> 

Azimuth 

127.8 

173.  C 

Eff,  Latitude 

34.5 

37.0 

Eff.  Longitude 

242.2 

272.2 

Nunber  of 
Residual  Pairs 


•*8 

53 

54 
54 
50 
'  9 

5:; 

46 

52 

52 

50 

5? 

50 

51 
57 
60 

59 

60 
57 

50 
46 
48 

51 
50 


I 


I 


TABIE  D.9 

Corr<2latlon  Coefficients,  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  1 
March  I968 

Correlation  Pair;  Urhana  ATS3  -  Sagamore  Hill  ATS3 


Number  of 

UP  Hour  Interval  6t^  ->''lg.  I  Alg.  II  6t^  Alg.  Ill  Residual  Pairs 


I 

I 


i 


I 

i 


8 


''f  1 1 

C. 78087 

0.90723 

C.  83617 

It  21 

C. 78188 

0.91531 

C.e5CC3 

2t  31 

0.65501 

0.90782 

C.701t4 

3t  *^1 

0.56C97 

0.90128 

C,  5*  (  ‘-P 

4 »  5) 

0.52638 

0.90165 

C. 53693 

5 1  61 

C. 51403 

0.91094 

0.52354 

6 »  M 

C. 38719 

0.89076 

C.40P'V8 

7 1  81 

C. 39609 

0.90595 

f .393C3 

8 1  ^1 

0.38569 

0*89671 

C.37f  22 

1C) 

C. 41275 

0.85532 

C. 49713 

1 v  t  ill 

0.54557 

0.76788 

C. 62677 

1 1 1 1  2 1 

C. 58506 

0.45435 

C.6276C 

!  2 1 1 3 1 

0.47667 

0.51753 

C.C5928 

13 1 1 41 

0.52780 

0.24408 

C.  174  61 

14,15) 

C.5m9 

0.44047 

C.5l«s6 

15  f 1 6) 

0.61  •»49 

0.54913 

C,4ifi67 

U  1  1  7  1 

0.62763 

0.64603 

r. 42383 

17,18) 

0.71454 

0.77492 

0.562  3'" 

1 8  f  1 ')  1 

0.74900 

0.90289 

C. 55553 

1<?  t  ?  C  1 

C.7B557 

0,94425 

f  .5  J8  9'' 

2  C"  1 2 1  1 

0.75P04 

0.96156 

C . 446C6 

21 t??) 

0.78299 

0.97308 

C  *  65645 

22t23) 

0. 80288 

C. 97033 

r  ,i7lOf.. 

2?f ?4) 

C.P63n8 

0.05863 

<“.73598 

Nominal  Value 

in  Degrees  for  Month 

• 

Urhana  ATS3 

Sagamore  Hill  ATS3 

Elevation 

39.0 

Azimuth 

173.6 

199.f' 

Eff,  Latitude 

37.0 

39.^ 

Eff.  Longitude 

272,2 

287.7 

45 
43 
48 
48 
43 
43 
51 
51 
51 
51 

<8 

50 

51 
*.4 
5^ 

51 

46 
45 
'.5 
45 
45 
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TABI£  D.tO 

Corr>ilation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  2 
January  1968 


Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 


Number  of 

Ui'  Hour  Interval  Alg.  I  6t^  Alg.  II  Alg.  Ill  Residual  Pairs 


1  0.  1) 

0.83529 

0.97550 

C. 90903 

1  it  2) 

0.95480 

0.97442 

C.  94575 

1  2t  3} 

0.96240 

0.97558 

0.93431 

1  3,  4) 

0.95407 

0.97317 

0.94696 

1  4,  5) 

0.89639 

0.96613 

0.95310 

1  5,  5) 

0.87755 

0.97536 

C.97533 

1  6,  71 

0.36595 

0.97113 

0.9402C- 

1  7,  8) 

0.80597 

0.  95775 

0. 87C77 

1  8,  9) 

0.50144 

0.95=t63 

0.864C5 

1  9,10) 

0.76836 

0.93954 

0. 84PC9 

1 10,  11) 

0.78554 

0.95918 

0.86845 

|Utl2) 

0. 72377 

0.93353 

0. 84100 

1  12, 13) 

0.72629 

0.93454 

C. 79605 

1 13, 14) 

0.77943 

0.96031 

0.63362 

114,15) 

0.59343 

0.90599 

C.3271? 

115, 16) 

0.93526 

0.64535 

C.  9364-=. 

115,17) 

0.82358 

0.83096 

0.92311 

117,13) 

0,87617 

0,36031 

C. 93663 

118, 19) 

0.93872 

0.39224 

0,93532 

119,20) 

0,95631 

0.954  38 

C.96r43 

120,21) 

0.93210 

0.94930 

0,92244 

121,22) 

0.90389 

0.92361 

0.90316 

122,  23) 

0.87370 

0.93966 

r. 37341 

123,24) 

0,79838 

0,93555 

O.RS^n 

nominal  Value  in  Degrees  for  Month 


, 

Stanford  ATS3 

Stanford  ATSl 

Elevation 

38.0 

37.) 

Azimuth 

139.5 

221.5 

Eff.  Latitude 

34.6 

34.5 

Eff.  Longitude 

240.7 

234.8 

71 
74 

74 

75 

76 
75 

75 

76 

74 

72 

75 

73 

75 
fO 

69 
73 
79 
79 

78 
PI 

79 

77 

76 

70 


I 

i 


I 


•* 

3 


j 

■» 

I 

I 


1 

I 
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TABIE  D.ll 

Corre  lation  Coefficients-  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  2 
February  I968 

Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 


Number  of 

UP  Hour  Interval  Alg.  I  Alg.  II  6t^  Alg.  Ill  Residual  Pairs 


i  0, 

1) 

0.96C68 

0.96361 

C  .  9  H  6  7  4 

79 

1  If 

2) 

0.96836 

e.96C0! 

0.96297 

63 

i  2f 

3) 

0.96957' 

0.95519 

G. 96672 

62 

1  3, 

4) 

0.92948 

C.Q2533 

C.925r-9" 

8? 

1  4, 

5) 

0.87526 

C.890B2 

r .90875 

73 

i  5t 

6) 

0.90816 

0.87648 

C..92023 

70 

I  6t 

7) 

0.95046 

0.98929 

C. 6-746 

65 

1  7f 

8) 

0.72016 

C. 89996 

C  .  P  3  5  8 

69 

i  B, 

9) 

0.70978 

0.93339 

0,84725 

70 

1  9,10) 

0.73302 

0, 93184 

C.E5P.7C 

63 

110,11) 

0.4P536 

C. 89956 

0 . 7  6 1  i  :» 

64 

111,12) 

0.4C497 

0.67998 

0.7  49  5  6 

69 

112,13) 

0.347H6 

C. 87659 

C. 62266 

6? 

113, lA) 

0.33182 

0.8B621 

G..£  13^-6 

71 

11A,15) 

0.32990 

C. 67588 

0.62267 

75 

115,16) 

0.77647 

C. 87595 

C.5  3£:i5 

76 

116,17) 

0.78250 

9.S5351 

C  .  d  r,i,  7  c. 

81 

117, IB) 

0.8744Q 

0.90132 

0.91402 

64 

118,19) 

C. 90087 

0.84583 

0.57262 

79 

|19,2U) 

0.91822 

0. 8691 3 

C. 3930  b 

79 

|?0,21) 

0.85950 

C.85r'l4 

0,P2f.*.5 

73 

121,22) 

0.85477 

0.P7'^q7 

C  •  0  9  6  7  9 

76 

122,23) 

0.89367 

C,9C736 

C,92i --5 

76 

123,24) 

0.90466 

0,90736 

f.  9301 

7t. 

nominal  Value  in  Degrees  for  Month 

Stanford  ATS3 

Stanford  ATSl 

Elevation 

32.6 

37.^ 

Azimuth 

029.7 

221.3 

Eff.  Latitude 

34.5 

34.^^ 

Eff.  Longitude 

241.9 

234.  ^ 
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TABLE  D.12 

Corre  lation  Coefficients-  of  Vertical  Time  Delay  Residuals 
Ev6ilua1;ion  Condition  2 
March  I968 

Correlation  Pair:  Stanford  ATS3  -  Stanfwd  ATSl 


OT  Hour  Interval 


AIe.  I 


6t^  Alg.  II 


6T^Alg.  Ill 


Number  of 
Residual  Pairs 


0,  n 

0.95124 

If  2) 

0.94834 

2t  3) 

0.92799 

?f  4) 

C.B7438 

4,  5) 

0.69P15 

5,  6) 

''.38115 

6f  7) 

0.44390 

7,  8) 

0.4  7674 

Bf  9) 

0,48438 

9,10 

0.41334 

1C, in 

0,37989 

11,1?) 

0.31999 

12,13) 

0,13518 

13,14) 

C.IP7J 1 

)4,15) 

C. 59781 

15,16) 

0.72325 

U,17) 

C.  856  30 

17,18) 

0.39971 

16,19) 

0.91414 

19, ?C) 

0,9.??41 

2C,?1  ) 

G. 86206 

?l,?2} 

0,92654 

22, ?3) 

0, 920^5 

23,54j 

C. 92605 

0.9 7794 
C. 97361 
C.9367C 
0. 69«!9 
0.9C269 
C.9C635 
0.926C0 
0.93449 
0.39568 
0.87630 
0.37603 
0.06345 
0.78764 
0.38355 
0.70459 
0.07C05 
0.91594 
0. 97453 
0.91016 
0.  8  7443 
C.946?'3 
•',9  nr  39 
0.95161 


C. 94088 
U. 93341 
C. 91255 
C. 84349 
0.73363 
C. 76402 
0.720  63 
0.68822 
0.63437 

0.6-:845 
0.67244 
0. 64514 
C. 59170 
C. 64396 
0.64918 
C. 70998 
0.81647 
C . 396- 1 
0.91130 
C.R9t05 
C.  329<^'0 
(.,92478 
C.9349Q 
C. 94257 


Nominal  Valu«^  ir  Degrees  for  Month 


Sf^ford  ATS3 

Stanford  ATSl 

Elevation 

31.3 

37.7 

Azimuth 

■  327.8 

221.2 

Eff.  latitude 

34.5 

34.'f 

Bff.  Longitude 

242.2 

234.H 

« 

S'  *. 

1 

1 

S 

1 

1 

TABIE  D.13 

k 

1 

Corr-ilation  Coefficients 

of  Vertical  Time  Delay  Pesiduals 

1 

i\ 

1 

t 

Evaluation  Condition  5 

1 

1 

► 

January  1969 

i 

1 

Corr 

elation  Pair:  Stanford  ATSl  -  Clark  Lake  AlTSl 

1 

i 

Number  of 

f 

( 

ITT  Hour  Interval 

6t^  Alg.  I  6t^  Alg.  II  6t 

y  Alg.  Ill 

Residual  Pairs 

1 

I  0,  1) 

0.95487 

0.95383 

-  14 

9 

i  It  2) 

0.95721 

0.82320 

0. 

54169 

15 

1 

» 

s 

\  2t  3) 

0.S1997 

0.51454 

0.55292 

15 

1 

1 

1  3t  4) 

0.95544 

0.91123 

0. 

94047 

15 

1 

i 

1  4,  5) 

0.96208 

0.7)575 

0.3713? 

15 

i 

1 

1  5t  6) 

0.83233 

0.79735 

0. 

50697 

15 

m 

0 

1  6,  7) 

0.65185 

0.57^43 

0. 

43250 

15 

§: 

r 

! 

1  7,  8) 

0.R2399 

0. 50789 

0. 

50146 

12 

1 

1  8»  9) 

0.55161 

0.57926 

0. 

40511 

18 

« 

1 

1  9fl0) 

0.55321 

0.43017 

0. 

50087 

18 

9 

1 

fin,  11) 

0.59098 

0.51395 

0. 

53253 

!B 

s 

1 

111 f  12) 

0.69354 

0.5755'^ 

0. 

59558 

18 

9 

1 

|]2tl3) 

0.84145 

0.91247 

0. 

80513 

IP 

9 

1 

113,14) 

0.93350 

0.89441 

0. 

86053 

IP 

.m 

114,15) 

0.601  0“^ 

0.72230 

0. 

55218 

IB 

fl 

1 

115,  1ft) 

0.80921 

0.9?939 

0. 

5725P 

1« 

1 

1 1  '=^,  17  ) 

0.39955 

0.53902 

0. 

52548 

16 

9 

1 

117, IP) 

0.93030 

O.Vll’O 

0. 

8103.? 

15 

9 

1 

1 18t  19) 

O.oiO^O 

0.95353 

0. 

91373 

14 

n 

1 1 9,20) 

0.96540 

0.9756? 

0. 

95166 

17 

1; 

1 

1 ’0, 21 ) 

0.32216 

0.P7569 

0. 

73085 

IP 

n 

1 

1 '1  ,22) 

0.81373 

0.917.54 

0. 

54315 

IP 

K- 

1 

1  ?2, 23) 

0.91236 

0.93173 

0, 

PP471 

16 

1 

1 

I23,2M 

0.^^513 

0.9‘917 

0. 

9519? 

14 

9 

1 

i 

• 

9 

1 

■Nominal  Value  in  Degrees  foi'  Month 

3 

1 

1 

1 

_ 

Stanford  ATSl 

Clark  Lake  ATSl 

• 

3 

'f 

\ 

Elevation 

37.7 

37.3 

S' 

* 

Azimuth 

221.1 

23c.  3 

1 

Eff.  Latitude 

34.6 

30.9 

* 

• 

r ' 

d 

Eff.  Longitude 

234.8 

240.1 

11 


I 
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TABIE  D,14 

Cori'elatioa  Coefficients-  of  Vertical  Time  Delay  Eesidvials 
Evaluation  Condition  5 
January  1969 

Correlation  Pair:  Stanford  ATSl  -  Fort  CoUihs  ATSl 


Ifl?  Hour  Interval  6t^  Alg.  I  6t^  Alg.  II 

t  Ot  1) 

0.45386 

0,91245 

0, 

!  1.  ?.) 

0.64558 

0, 76800 

0 

1  2»  3) 

0.59731 

0.52916 

0. 

1  3,  4) 

0.54796 

0.44613 

c. 

1  4,  5) 

0.6<775 

0,57516 

0. 

1  5,  81 

0.55332 

0.50636 

0, 

1  6,  7) 

0.59405 

0.52726 

0 

1  7,  9) 

O.'BOPO 

0.56131 

0 

1  8»  91 

0.20929 

0.75091 

0 

1  9,10) 

0.20074 

0.13759 

0 

1 10, 11) 

0,17136 

0.15264 

0 

tn,i?) 

0.27676 

0.29771 

0 

(12,  13.) 

0.26709 

0.33416 

0 

in,  14) 

0.29339 

0.42861 

0 

114,15) 

-0,07934 

0.39112 

0 

in, 15) 

0.32901 

0.37119 

0 

116,17) 

0.21631 

0.55015 

0 

1)7,18) 

0.42*02 

0.74145 

0 

119,19) 

0.95141 

0.93991 

0 

119,20) 

0,75196 

0.93536 

0 

120,21) 

0.67397 

0.901  O-' 

0 

121  ,22) 

0.74781 

0.8471? 

0 

|?7,2-») 

0,01513 

0.77069 

0 

173,7’-) 

0.90‘^'»2 

0.67085 

0 

nominal  Value  in  Degrees  for  Month 

*, 

Stanford  ATSl 

Fort  Co3.3.ins  ATSl 

Elevation 

37.7 

25.«o 

Azimuth 

221.1 

236.8 

Eff,  Latitude 

3^.6 

37.0 

Eff.  Longitude 

234.8 

248.:? 

Nunibe?'  of 
Residual  Pairs 


80 

94 

90 
84 
94 
93 
81 

91 
91 

78 

70 

PI 

91 

84 

84 

81 

p«. 

76 

79 
9! 
90 
7P 

77 


I 

i* 

n 

I 

4 


I 
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TABLE  D.  15 

Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 
January  I969 

Correlation  Pair:  Rosman  ATS3  -  Urtana  ATS3 


UT  Hour  Interval 


6T  Alg.  I 


6t  Alg.  II 


6t  Alg.  Ill 


Number  of 
Residual  Bairs 


O.P9«69 

3.96173 

0.35770 

0.84107 

0.85292 

0.77327 

0.7700? 

0.70367 

0.69997 

0.5665? 

0.56254 

0.65617 

0.53651 

0.72156 

0.80759 

0.33777 

0.88290 
0.34229 
0.93416 
0.931?'. 
0.93729 
0.9?7?2 
0,?'f  3? 
0.93')82 


0.01204 
0.88737 
0.37 >93 
0.93793 
0.35990 
0.77747 
0.71246 
0.*«f  961 
0.69397 
0.57361 
0.65197 
0.72739 
0.27765 
0.71297 
0.66045 
0.34365 

0.94453 

0.9584? 

0.08307 

0.97977 

0.96051 

0.93457 

0.96869 

0.95427 


0.90733 
0.91360 
0.91716 
0.37710 
0.36443 
0.74593 
0.63635 
0.62349 
0.65664 
0.52383 
0.55012 
0.f53?6 
0.51046 
0.64535 
0.65505 
0.73024 
0.3624? 
p.341 34 

0.93455 

0.92705 

0.91335 

0.90219 

0.9303? 

0.91144 


TJominal  Value  in  Degrees  for  Month 


Rosman  ATS 

3 _ 

L  _ 

Urbana  ATS 

Elevation 
Azimuth 
Eff.  Latitude 
Eff.  Longitude 


46.7 

156.9 

32.6 

278.4 


39.6 

151.8 

37.C 

273.8 


1 

1 

TABLE  D.i6  •  I 

•  JV 

Correlation  Coefficients  of  Vertic  il  Time  Delay  Residuals  J 

Evaluation  Conditi.on  5  J 

Januaiy  1969  1 

Correlation  Pair:  Arecibo  ATS3  -  Cold  Bay  ATSl  J 


OT  Horn:’  Interval 

6t^  Alg.  I 

6t^  Alg.  II 

6T^Alg.  Ill 

Number  of 
Residual  Bairs 

1  0,  1) 

0.4511.3 

0.31124 

0.30969 

78 

t 

1  It  2) 

0.2 74  OS 

0.31 356 

0.338.57 

• 

1 

1  ?,  3) 

0.4Q717 

0.22783 

0.5S733 

n 

■> 

1  3t  4) 

0.33U4 

0.22524 

0.25045 

72 

1  4,  5) 

0.01135 

0.31448 

-0.31244 

70 

1 

1  5,  ^>) 

0.01183 

0.20774 

0.05704 

69 

i 

1  6t  7) 

-0.33431 

-0.02233 

-0.18494 

72 

} 

1  7,  8) 

-0.50112 

-0.21173 

-0.39233 

»1 

1  8t 

-0. <48907 

-0.34402 

-0.4314? 

71 

i 

1  ‘>tl0) 

-0.31 5?5 

-0.37238 

-3.27297 

73 

jk 

■f 

no, in 

-0.33737 

-0.44353 

-0.23413 

7? 

1 

in  ,1?) 

-0.23321 

-0.35775 

-0.19825 

78 

i 

I!2t!3» 

-0.1'i 27 

-0.00141 

-0.22-23 

76 

i 

113,14) 

-0.33034 

0.30235 

-0.1292? 

70 

114,18) 

O.OS383 

0.3317? 

-0.0145O 

92 

i 

m,  u) 

0.14394 

0.317^1 

0.34231 

79 

5 

116,17) 

0,1 8345 

0.?H340 

0.02242 

79 

4 

.|17,1P) 

0.32086 

0.??001 

•  C.1»344 

73 

'? 

118,  n) 

0.23633 

-0.294^7 

3.13712 

71 

i 

11", ?0) 

0.2514? 

-0.’4574 

0.17484 

74 

l?0,?l) 

0.00496 

0.18)  13 

0.12735 

7/* 

i 

12!,??  ) 

-3.13123 

0.71490 

0.14447 

70 

; 

12?, ?3) 

-0.1 9303 

0.12487 

0.17117 

77 

• 

1  ?3,  ?4) 

0.1 ?106 

0.38150 

0.38700 

7B 

nominal  Value  in  Degrees  for  Month 


TABIE  D  .18 

Corr*ilation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 
April  1969 

Correlation  Pair:  Stanford  ATSl  -  Fort  Collins  ATSl 


UP  Hour  Interval 


6t^  Alg.  I 


6%Alg.  11 


6T^Alg.  Ill 


Number  of 
Residual  Bairs 


i  c>  n 

0.89298 

C.9C993 

G.  3‘9339 

1  It  2} 

0.86817 

0.9AA60 

•?.a7A74 

1  2t  3) 

0.80R8A 

0.89570 

C. 3^235 

1  3,  4) 

C.90A99 

C.7588A 

0.38612 

I  At  5) 

0.913A6 

0.8A66C 

0.92745 

I  6t  6} 

C. 81762 

C.C35A9 

0.95C46 

1  6t  71 

0.67195 

G.82C32 

0.69120 

1  7t  8) 

0.7C652 

0.87799 

0.72269 

1  8t  9) 

0.6667A 

0.89649 

0.69995 

1  9tlC) 

C.512A6 

0.83524 

0.55634 

1 ICt 11) 

C. 56589 

0.92239 

0.63157 

Illtl2) 

C.7C3A5 

C. 93451 

''.75566 

|12t 13) 

C.5137A 

0.9^198 

0.57786 

( 13t lA) 

C. 05592 

-0.26118 

0. 15395 

|IA,15) 

-0.03259 

0.53086 

30  25 

1 15tl6) 

C.7Afc7C 

0.63436 

0.5457) 

|16tl7) 

0.81179 

C.6126C 

c.7:>)2e 

117,18) 

C.856A3 

0.60730 

?. 82553 

1  18,19) 

C.67C87 

C. 77391 

0.8*088 

119,20 

0.9C217 

C. 72881 

0.67906 

120,21) 

C.9C299 

0.76915 

0.85544 

121,22) 

0.9A3A5 

0.85)14 

9*396 

122, 23) 

C. 95589 

0.95664 

n. 05159 

123, 2A) 

0.95O35 

C. 89*94 

C.^523? 

Nominal  Vttlue  in  Degrees  for  Month 

' 

Stanford  ATSl 

Port  Collins  ATSl 

Elevation  37,9 
Azimuth  220.8 
Eff.  Latitude  34.6 
Eff.  Longitude  234.9 


25.2 

236.6 

3?.0 

248.6 


TABIE  D.I9 

Correlation  Coefficients-  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 
Apr'l  1969 

Correlation  Pair:  Rosman  ATS3  -  Urbana  ATS3 


Number  of 


1  0 

1) 

i  1 

2) 

1  2 

3) 

1  3 

4) 

1  ^ 

5) 

i  3 

6) 

1  6 

7) 

1  7 

8) 

1  3 

9) 

1  9 

10) 

1 1C 

11) 

til 

12) 

1  12 

13) 

113 

14) 

lU 

15) 

1  13 

16) 

116 

17) 

117 

18) 

1 16 

19) 

119 

20 

120 

21) 

|21 

22) 

122 

23) 

123 

24) 

6t^  Alg.  I 

Alg.  II 

6t^  Alg.  Ill 

Residua. 

C. 92578 

0.95864 

0. 

39 

0.9C2C7 

0.94884 

0.  9‘NlOq 

39 

C. 87555 

C. 94141 

0.  8 7*^93 

38 

0.65331 

0.913C7 

9.8,5335 

36 

C. 83596 

0.90163 

•'>.8351 5 

38 

C. 72912 

0.87103 

0.72569 

39 

C.6489^ 

0.85945 

''.6616^ 

37 

0.6054C 

0.86059 

r.6iy3-‘' 

37 

0.66268 

O.f  3C7C 

r.6-327'-- 

38 

C. 66634 

0.72591 

0.65-'48 

36 

0.65231 

0.77423 

8.65263 

34 

C. 76804 

0.82481 

0.75391 

38 

C.6565C 

0.76415 

r..  61 154 

36 

C.7C117 

0.72275 

0.65931 

36 

0.04803 

0.84553 

3.31078 

36 

C.37C75 

C. 85654 

?. 34^50 

36 

C. 90559 

C. 91545 

094?? 

36 

C. 89946 

0.87922 

17043 

34 

C. 93089 

C. 92320 

.).9?'»27 

35 

C. 93473 

C. 94429 

'•.03255 

33 

C.945C4 

C. 95605 

94418 

34 

0.93299 

C. 95179 

C. 932?? 

32 

C. 92778 

C. 96928 

•'♦''277^ 

36 

0.91t5C 

C. 97916 

''.•>1618 

36 

T\’onj3nal  Value  in  Degrees  for  Month 


Elevation 
Azimuth 
Eff,  Latitude 
Eff,  Longitudo 


Rosman  ATS3 

47.8 

162.3 

32.6 

278.1 


Urbana  ATS3 
40.9 

156.4 

37. '> 

273.4 
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TABIE  D20 

Conelation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  ^ 

April  1969 

Correlation  Pair:  Stanford  ATSl  -  Stanford  ATS3 


ITT  Hour  Interval 


6\  Al€.  I 


6T^Alg.  II 


6T^Alg.  Ill 


Number  of 
Residual  Bairs 


c 

1) 

1 

2) 

2 

3) 

3 

4) 

4 

5) 

5 

6) 

6 

7) 

7 

8) 

8 

9» 

9 

10) 

1C 

11) 

11 

12) 

12 

13) 

13 

141 

14 

15) 

15 

16) 

16 

17) 

17 

13) 

le 

19) 

19 

20) 

2C 

21 ) 

21 

2?) 

22 

23) 

23 

24) 

O.S3?77 
0.94376 
C.94CC6 
0.92604 
0.9C655 
0.87271 
C. 86461 
C.86C02 
0.83395 
0.83743 
C. 81012 
0.7853C 
0.79214 
C.792B0 
C. 92229 
€.90046 
C.333d2 
C. 82523 
0.83890 
0.8430C 
C.B56-37 
C. 67061 
C. 89216 
0.90549 


0.85380 
0.91631 
0.95538 
0.90499 
0.90382 
0.89920 
0.89228 
C..93435 
0.939P3 
C. 88724 
C. 88855 
C. 88637 
0.91823 
0.33766 
0,96216 
0.92643 
0,34677 
C.774Ce 
0.77754 
0.73964 
0.72900 
0,70635 
C. 71029 
0, 78109 


^.94229 

'‘,94545 
'*.93074 
?. 92692 

Q'Mr? 

f. 83635 
83746 
34618 
''.83194 
r-.  79567 
73958 
'‘.69818 

'‘.75977 

1',  q7nS5 

•'‘.3472'^ 

343 5''' 
8  '  4  08 
0,87326 
87643 
(.07008 
87816 
99'»4C 
9f'90f- 

“.92744 


Nominal  Value  in  Degrees  for  Month 


Stanford  ATSl 

Stanford  ATS 

Elevation 

37.9 

22.0 

Azimuth 

220.8 

n7.3 

Eff.  Latitude 

3^.6 

34.3 

Eff.  Longitude 

234.9 

244.5 

>250 


TABIE  D.P1 

Corielation  Ctjefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 
July  1969 

Correlation  Pair;  Stanford  ATSl  -  Clark  Lake  ATSl 


UP  Hour  Interval 


6t^  Alg.  I 


6T^Alg.  II 


6T^Alg.  Ill 


Nvuriber  of 
Residual  Pairs 


0* 

1) 

It 

2) 

2t 

3) 

3t 

4) 

4t 

5) 

5f 

6) 

6t 

7) 

7, 

81 

8t 

9) 

9t 

10) 

lOt 

111 

lit 

12) 

12t 

13) 

I3t 

14) 

14t 

15) 

15* 

16) 

I6t 

17) 

17, 

18) 

18, 

19) 

19, 

20) 

20, 

21) 

21, 

22) 

22, 

23) 

23, 

24) 

0.91517 

0.88735 

0.90767 

O.B6P35 

0.62018 

0.46142 

0.38519 

0.53334 

0.92500 

0.91454 

0.95127 

0.95584 

0.67821 

0.85085 

0.81672 

0.90773 

0.91189 

0.93285 

0.94716 

0.93953 

0.96158 

0.96350 

0.96722 

0.94535 


0,‘>0396 

0.90097 

0.93134 

0.88358 

0.81251 

0.72631 

0.79539 

0.73388 

0.89654 

0.82443 

0.93007 

0.93666 

0.77900 

0.37531 

0.78747 

0.37937 

0.80475 

0.91032 

0.92529 

0.90317 

0.94593 

0.94562 

0.95799 

0.94041 


nominal  Value  in  Degrees  for  Month 


D.B6139 

0.97940 

0.95393 

0.76910 

0.54842 

0.37431 

0.65363 

0.67395 

0.92563 

3,99903 

0,''4674 

C. 05299 

0.5^55? 

0.8Q‘^11 

0.62312 

3,93934 

0.96143 

0.91551 

0,91967 

0,97050 

0.9?096 

0,91930 

0.97961 

0.90994 


Elevatio. 
Azimuth 
Eff.  Latitude 
Eff.  Longitude 


Stanford  ATSl 

38.1 

220.U 

34.6 

234.9 


Clark  Lake  ATS] 

37.6 

2':9.8 

30.9 

240.2 
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TABIE  D.2P. 

Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 

July  1969 

Correlation  Pair;  Stanford  ATSl  -  Fort  Collins  ATSl 


Hour  Interval 

6t^  Alg.  I 

6t^  Alg.  n 

6t^  Alg.  I 

i  Of 

1) 

0.53855 

0.60406 

''.407^;) 

1  li 

2) 

0.21069 

0.37332 

M5722 

1  2t 

3) 

0.50904 

0.54935 

24^77 

i  3f 

4) 

0.53682 

0.52630 

?7?3f 

i  4t 

5) 

0.40843 

0.45586 

f'.  14S74 

1  5f 

6) 

0.10644 

0.32095 

-".15675 

1  6f 

7) 

-0.01535 

0.35473 

-C. 19470 

1  7f 

8) 

0.09337 

0.41491 

i  a« 

9) 

0.11446 

0.45010 

-'■.'•■7569 

1  9f 

10) 

0.07267 

0.45835 

-0.  12285 

|10« 

11) 

•  0.08203 

0.51787 

*.36  93 

lilt 

12) 

0.06797 

0.52006 

-  '».f5'‘14 

il2f 

13) 

0.05473 

0.36485 

-•).  12477 

Il3t 

14) 

0.10830 

0.23369 

-".19971 

114, 

15) 

0.46883 

0.43398 

".191  59 

115, 

16) 

0.54207 

0.41452 

••*,43584 

1  16, 

17) 

0.64704 

0.54986 

".63';  16 

117, 

18) 

0.71092 

0.64967 

9,72859 

118, 

19) 

0.77294 

0.69369 

<•.72775 

119, 

20) 

0.7347B 

0.73081 

•',73044 

120, 

21) 

0.77919 

0.75350 

'•.77148 

i  21, 

2?) 

0.77737 

0.77293 

•'.91357 

122, 

23) 

0.78466 

0.31374 

■*.  91  ■’55 

123, 

24) 

0.74600 

0.51565 

746 2 8 

T^ominal  Value  in  Doprees  for  Month 


Stanford  ATSl 

Fort  Collins  ATSl 

Elevation 

Azimuth 

Eff,  Latitude 

Eff.  Longitude 

38.1 

220.  U 

34.6 

234.9 

2S.4 

236.3 

37.0 

248.7 

Number  of 
Residual  Pairs 


TABIE  D.  23 

Correlation  Coefficients  of  Verticil  Time  Delay  Residuals 
Evaluation  Conditj.on  5 
July  1969 

Correlation  P&ir:  Roaman  ATS3  -  Urbana  ATS3 


Number  of 


UP  Hour  Interval  6t^  Alg.  I  6t^  Alg.  II 

6t^  Alg.  Ill 

Residual  I^irs 

*- 

1 

•s 

i  0, 

11 

0.76674 

0.81496 

0,71^65 

80 

-1 

i  1* 

2) 

0.77996 

0.83158 

'■.73:75 

80 

3 

3 

1  2, 

3) 

0.83121 

0.84192 

'54?6 

80 

> 

i  3, 

4) 

0.73323 

0.71431 

r,6''687 

79 

i 

i  4, 

5) 

0.70158 

0.60920 

r. 54673 

83 

1  5, 

6) 

0.65786 

0.53779 

'■.r 

55766 

S3 

1 

* 

1  &t 

7) 

0.66887 

0.66  .'la 

3,6156^ 

80 

1  7, 

8) 

0.60994 

0.49554 

C. 4511  7 

87 

1 

t  8, 

9) 

0.55450 

0.54453 

r. 37719 

82 

1  9t 

13) 

0.43936 

0.54702 

•3,  376'^1 

82 

110, 

11) 

0.53341 

0.63172 

3'’95^> 

88 

1 

ill, 

12) 

0.69170 

0.72514 

53576 

87 

112, 

13) 

0.73627 

0.75396 

''.674  7? 

87 

*i 

113, 

14) 

0.83672 

0.32539 

I?.  81144 

S3 

1 

114, 

15) 

0.82867 

0.3300^ 

Q 

/  •  C 

'•994 

81 

(15, 

16) 

0.84175 

0.82656 

»»1568 

8? 

1 

-e 

116, 

17) 

0.64423 

0.87845 

'".82139 

88 

! 

i 

117, 

18) 

0.85320 

0.91092 

.r.f 

>3840 

86 

1 

|ia. 

19) 

0.35393 

0.39224 

''.83^65 

87 

■i 

119, 

20) 

0.65460 

0.39606 

83791 

S3 

j 

l?0. 

21  ) 

0.84119 

0.83781 

0.  =<1863 

82 

4 

121, 

22) 

0.83337 

0.3P402 

-Kgo 

79 

* 

122, 

23) 

0.64610 

0.87317 

'••.•^??94 

67 

123, 

2V) 

0.83346 

0,38593 

'..79  ".7  C 

79 

1 

6 

i 

V 

3 

nominal  Value  in  Degrees  for  Month 

• 

4' 

3 

Rosman  ATS3 

Urbana  ATS3 

't 

% 

i 

Elevation 

33.6 

26.7 

4 

Azimuth 

127.3 

125.3 

i 

Eff.  Latitude 

32.5 

36.8 

1 

Eff.  Longitude 

281.1 

277.3 

'A 

% 

■3 

.  .. 

3 

I 

I 

a 


"i 


'i 

.# 
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TAEIJi:  D.1’4 

Correlation  Coefficients  of  Vertical  Time  Delay  Residuals 

Evaluation  Condition  5 
October  I969 

Correlation  Pair;  Stanford  ATSl-Clark  Lake  ATSl 


m  Hour  Interval  6t^  Alg.  I  6t^  Ale.  II 


Alg.  Ill 


Number  of 
Residual  Pairs 


I  Of  II 
I  If  2) 
I  2f  3) 
I  3f  4) 
I  4f  51 

(  5f  6) 

I  Af  71 
I  7f  6) 
I  8f  9} 
I  9fl0) 
110,111 
lilt  12) 
112, 131 
113,141 

114.15) 

115.16) 

116.17) 

117.18) 
118,  19) 

119.23) 
120,21) 
121,22) 

122.23) 

123.24) 


0. 99007 

0.99236 

0,95949 

0.92699 

0.95085 

0.88822 

0.84978 

0.80639 

0.78266 

0.73569 

0.72884 

0.30016 

0,85673 

0.80117 

0,85890 

0.68C05 

0.82293 

0,94326 

0.95259 

0,96410 

0,92583 

0,96552 

0.98134 

0.98417 


0.97365 

0.98869 

0.98403 

0.96073 

0,97495 

0,97094 

0.94750 

0.94695 

0.94551 

0,94211 

0.90213 

0.93064 

0.92519 

0,91806 

0.762 

0,92351 

0,85311 

0,90738 

0.97299 

0,97487 

0,93506 

0,96524 

0,95002 

0,91113 


r. 95’ 87 

0,9585’ 

0.93556 

*.95935 

0,95189 

''.9368J 

O,  92767 
''.94760 
■).  94566 

C. 9 0721 

P,  01 7-2 

'',96201 

0.9*]  26 
0. 9365* 
0.95flfi 
”,95'‘** 
0,93857 
'*.510  62 
9287? 
^.94*^66 
'■,95805 


42 

42 

42 

39 

37 

36 

39 

36 
26 
24 
21 
23 
29 

37 
42 
42 

38 

38 

39 
42 
42 
42 
37 
39 


Nominal  Value  in  Degrees  for  Month 


Stanford  ATSl 

Clark  Lake  ATSl 

Elevation 

38,3 

37.9 

Azimuth 

219.9 

229.4 

Eff.  Latitude 

34.6 

30.9 

Eff,  Longitiaie 

235.0 

240.3 

-254- 
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TABLE  D.  '■^5 

Conflation  Coefficients  of  Vertical  Time  Delay  Residuals 

Evaluation  Condition  5 
October  1969 

Correlation  Pair:  Stanford  ATSl  -  Port  Collins  ATSl 


OT  Hour  Interval 


6t^  Alg.  I 


6T^Alg.  11 


6t^  Alg.  Ill 


Number  of 
Residual  Pairs 


p 

5^ 

1  0,  1) 

0.94390 

0.86504 

/\ 

86526 

74 

i  It  2) 

0.93363 

0.89578 

C 

87357 

72 

i 

1  2i  3) 

0.82730 

0.85757 

r> 

55411 

72 

1  3t  4} 

0.82968 

0.85308 

r- 

P6P65 

72 

1 

1  At  5) 

0.84588 

0.B9923 

9'>?54 

72 

1' 

1  5t  6) 

0.75661 

0.83981 

0 

69 

1  6t  71 

0.58350 

0.66985 

69142 

69 

1  7t  8) 

0.55273 

0.75278 

7H24 

75 

1 

i  St  9) 

0.51570 

0.76716 

f*. 

73215 

75 

.§ 

1  9tl0) 

0.35815 

0.72709 

0 

742  66 

75 

llOtll) 

0.24700 

C. 71047 

p 

68837 

71 

I 

Illtl2) 

0.35459 

0.72343 

n 

67?''S 

69 

1 

Il2tl3) 

0.37106 

0.73437 

0 

71860 

69 

1 

ll3tlA) 

0.12932 

0.44684 

57293 

72 

llAtl5) 

0.28574 

0.30563 

n 

67278 

76 

.1 15t  16) 

0.73495 

0.66471 

6??4'' 

80 

116,17) 

0.88755 

0.64031 

/* 

7713? 

77 

1 

117,18) 

0.92485 

0.80416 

(, 

y3P02 

77 

118,19) 

0.92948 

0.84866 

0 

70444 

78 

119,20 

0.93336 

0.86555 

0 

73J  35 

78 

1 

120,21) 

0.93073 

0.597.35 

?( 

73431  ' 

78 

121,22) 

0.Q4559 

0.89165 

r. 

81260 

74 

122,23) 

0.95458 

0.02487 

f 

0403: 

77 

f 

1 

123,24) 

0.96186 

0.89289 

30(f,4  7 

71 

■ 

nominal  ydue  in  Degrees  for  Month 

Stanford  ATSl 

_ Pbrt  Collins  ATSl 

Elevation 

38.3 

25.6 

Azimuth 

219.9 

2J5.9 

1 

Eff,  Latitude 

34.6 

37.0 

Eff.  Longitude 

235.0 

248.8 
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TABLE  D.26 


Correlation  Coefficients-  of  Vertical  Time  Delay  Residuals 
Evaluation  Condition  5 
October  I969 

Correlation  Pair:  Rosman  AT33  -  Urbana  ATS3 


UT  Hour  Interval  6t^  Alg.  I  Alg.  II 


Number  of 

Alg.  Ill  Residual  Pairs 


Of  1) 

0.93251 

0.94330 

•3.  90439 

32 

If  2) 

0.95587 

0.95528 

0,93703 

27 

2f  3) 

0.97053 

0.96401 

0,9*^525 

27 

3f  4) 

0.97412 

0.96473 

0.964  51 

32 

4f  5) 

0.98680 

0.98949 

0.97616 

29 

5f  6) 

0.97239 

C. 97936 

0,97799 

29 

6f  7) 

0.95915 

0.97678 

0,  Q#'-47'' 

29 

7f  8) 

0.96246 

0.98829 

C. 99067 

30 

8f  9) 

0.91202 

0.97934 

o,97?35 

29 

1 

9f  10) 

0.37333 

0.96369 

C. 96191 

30 

j 

lOf 11) 

0.9C226 

0.97053 

■'.9607! 

30 

! 

Ilf  12) 

0.85170 

0.Q9442 

o,04?41 

.32 

■ 

I2f 13) 

0.90566 

0,96366 

0,91 ?69 

35 

1 

13f 14) 

0*86341 

0.91505 

3.935O0 

33 

1 

14f 15) 

0.91689 

0,91772 

0 ,  967''5 

33 

i 

15f 16) 

0.86003 

0.90511 

0,03037 

32 

) 

1 

16f 17) 

0,91574 

0,94521 

7,92707 

36 

{ 

17f 18) 

0,91679 

0,94502 

0,91719 

36 

{ 

I8f 19) 

0.92138 

0.94161 

0.8  32  *'>3 

36 

1 

19t  20) 

0.95179 

C. 97099 

?, 93273 

32 

\ 

f 

20f 21 ) 

0.94722 

0,97937 

0.,  03149 

36 

21f 22) 

0,88891 

0.96455 

0, ?7t  16 

53 

22f  23) 

0,82.323 

0.95479 

n, 

35 

; 

23.24) 

0,93477 

0.90505 

79529 

55 

1 

"Nominal  Value  in  Degrees  for  Month 


Elevation 
Azimuth 
Eff.  Latitiide 
Eff.  LongitMe 


Rosman  ATS3 

33.0 

126.6 

32.4 

281.2 


Urbana  ATS; 

26.4 

12'n9 

56.8 

27^.4 
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APPENDIX  E 

TABIES  OF  CORREIATION  COEFFICIENTS  OF  VERTICAL  TIME  DELAY  - 
OVER  ENTIRE  EVAIAJATION  INTERVAL 


TADU3E.1 

Correlation  CoGlficiontc  of  Vorticnl  Time  Delay  Residuals 

Over  Entire 

Evaluation  Condition  1 


Number  of 

6t  Alg.  I  6t  Alg.  II  6t  Alg.  Ill  Residual  Pairs 

V  V _ _ V 

Correlation  Pair:  Stanford  ATS3  -  Stanfoi*d  ATSl 

O.S7915  0.89689  0.88513-  5334 


Correlation  Pair:  Stanford  ATS3  -  Urbana  ATS3 

0.73917  0.68118  0.57075  4657 

Correlation  Pair:  Urbana  ATS3  -  Sagamore  Hill  ATS3 

0.83953  0.79823  .O.806C6 
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TABIE  E.2 


Correlation  Coefficients-  of  Vertical  Time  Delay  Residuals 

Over  Entire 

Evaluation  Condition  2 


Numbo"’  of 

6t^  Alg.  I  Alg.  II  6t^  Alg.  Ill  Residual  Pairs 


Correlaticxi  Pair;  Stanford  ATS3  -  Stanford  ATSl 


TABIi:  E.3 

Correlation  CoefiicientG  of  Verticil  'Ji'jmc  l»c].ay  HocjAiuula 

Over  Entire 

Evaluation  Condition  5 


6T„  Alg.  I  6t  Alg.  II  6t  Alg.  Ill 


Nuniber  of 
Residual  Pairs 


Correlation  Pair:  Stanford  ATSl  -  Stanford  ATS3 


0,90284 


0.82314 


0.90038 


Correlation  Pair;  Stanford  ATSl  -  Clark  Lake  ATSl 


0.93295 


0. 9013Y 


0. 09543 


Correlation  Pair:  Stanford  ATSl  -  Port  Collins  ATSl 


0. 80925 


0. 71691 


0. 74504 


Correlation  Pair:  Rosraan  ATS3  -  Urbana  ATS3 


0.86105 


0.83425 


p.  87130 


Correlation  Pair;  Arecibo  ATS3  -  Cold  Bay  ATSl 


0.03723 


0.00581 


-0.08283 


TABLE  F.l 


"Uunuiativo  Frequency  Distribution  cf  the  Daily  Correlation 
Coefficients  of  the  Vortical  Time-  Delay  Residuals  over 
Eval.uatlon  Condit.’  on  1 


Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 
Total  Number  of  Correlation  Coefficients:  85 


Correlation  Coefficient  Cumulative  Frequency  Distribution  (Percent) 

Interval  6t  Alg.  I  6t,,  Alg.  Z1  6t  Alg.  Ill 


TABIE'P.2 

"TumuJative  Frequency  Distributioi  of  the  Daily  Correlation 
Coefficients  of  the  Vertical  Ti.rc  Delay  Residuals  over 
Evaluation  Condition  1 

Correlation  Pair:  Stanford  ATS3  -  Urbar^-.  fV?S3 
Total  Number  of  Correlation  Coefficients:  76 

Correlation  Coefficient  . Cumulative  Frequency  Distribution  (Percent) 

Interval  6t^  Alg.  I  6Ty  Alg.  II  6t^  Alg.  Ill 


0 

.0 

0 

0 

1.32 

0 

• 

1 

H 

1 

0 

1.32 

0 

• 

I 

H 

1 

0 

1.32 

0 

0 

2.63 

1.32 

0 

3.95 

1.32 

•0 

3.95 

2.63 

0 

5.26 

2.63 

1.*’' 

7.89 

3.95 

l"l,0) 

3.95 

7.89 

5.26 

3.95 

11.84  •• 

10.53 

6.58 

13.16 

17.11 

9.21 

17.11 

23.68 

9.21 

31.58 

31.58 

14.1i7 

39.47 

38.16 

23.63 

46.05 

55.26 

32.89 

61.84 

65.79 

53.95 

78.95 

81. 

69.47 

92.11 

94.74 

(-1,  11 

100.00 

100.00 

100.00 
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TABIE  F.3 

Cumulative  Frequency  Distribution  of  the  Daily  Correlation 
Cofcfficients  of  the  Vertical  Time.  Delay  Residuals  over 
Evaluation  Condition  1 

Correlation  Fair:  Urbana  ATS3  -  Sagamore  Hill  ATS3 
Total  Number  of  Correlation  Coefficients:  68 

Correlation  Coefficient  •  Cumulative  Frequency  Distribution  (Percent) 

Interval  Alg.  I  6t^  Alg.  II  Alg.  Ill 


0 

0 

0 

0 

0 

0 

0 

0 

0 

• 

1 

H 

1 

0 

0 

0 

[-l,>.5) 

0 

0 

0 

[.l,-.4) 

•  0 

0 

0 

[ -1,~.3) 

0 

0 

0 

CVJ 

• 

1 

•H 

1 

0 

0 

0 

0 

1.47 

0 

1-1,0) 

0 

1.47 

0 

[-l,.l) 

0 

1.47  ■ 

0 

l-l,*2) 

0  ■ 

1.47 

0 

l-l,.3) 

1.4? 

5.88 

1.47 

[-1,.4) 

1.47 

8.82 

2.94 

[-1,.5) 

1.J|7 

13.24 

5.88 

[-1,.6) 

4.41 

19.12 

8.82 

(-l,.7) 

10.29 

29.41 

20.59 

[-1,.8) 

16.18 

45.59 

38.24 

l-l,.9) 

5'<.4l 

77.94 

73.53 

[-1,  1] 

JO'J.OO 

JUO.UU 

lOli.OO 
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TABIE  p,4 

Cumulative  frequency  Distribution  of  the  Daily  Correlation 
®o®^icients  of  the  Vertical  Time  Delay  Residuals  over 
Evaluation  Condition  2 


Correlation  Pair:  Stanford  ATS3  -  Stanford  ATSl 
Total  Number  of  Correlation  Coefficients:  88 


Correlation  Coefficient 
Interval 


•  Cumulative  Frequency  Distribution  (Percent) 


6t  Alg.  I 


6Ty  Ale.  II 


6TyAl€.  Ill 


TABIE  P.5 

CuB'ulative  J’requency  Distribution  of  the  Daily  Correlation 
Coefficients  of  the  Vertical  Time  Delay  Residuals  over 
Evaluation  Condition  5 

Correlation  Pair:  Stanford  ATSi  -  Clark  Lake  ATSl 
Total  Nuniber  of  Correlation  Coefficients:  69 

Correlation  Coefficient  •  Cumulative  Frequency  Distribution  (Percent) 

Interval  6t^  Alg.  I  6^^  Alg.  II  Alg.  Ill 


0 

0 

0 

CO 

• 

1 

H 

1 

0 

0 

0 

• 

1 

H 

1 

0 

0 

0 

0 

0 

0 

t 

H 

t 

0 

0 

0 

1-1,-.!^) 

0 

0 

0 

0 

0 

0 

l-l,-.2) 

0 

0 

0 

1 

•  0 

0 

0 

1-1,0) 

0 

0 

0 

t-l,.l) 

1.45 

0 

1.45 

[-1,.2) 

1.45 

0 

1.45 

[-1,.3) 

1.45 

0 

1.45 

l-l,.4) 

2.90 

1.45 

1.45 

t-1,.5) 

4.35 

2.90 

2.90 

[-1,.6) 

5.80 

5.80 

5.80 

[-1,.?) 

10.14 

8.70 

7.25 

l-l,.8) 

15.94 

i5.r'i 

21.74 

(-5,. 9) 

r.o 

.>  1  •  1 

'’iO.:/'. 

t-1,  n 

100.00 

-266- 

100.00 

100.00 

TABI£  F.6 


Cumulative  Frequency  Distribution  of  the  Daily  Correlation 
Coefficients  of  the  Vertical  Time  Delay  Residuals  over 
Evaluation  Condition’  5 


Correlation  Pair:  Stanford  ATSl  -  Fort  Collins  ATSl 


Total  Number  of  Correlation  Coefficients:  95 


Correlation  Coefficient 
Interval 

t-l,-.9) 

t-l,-.8) 

[-1,-7) 

[-l,-.6) 

t-l,".5) 

l-l,-.3) 

t-l,-.2) 

1-1,0) 

l-l,.l) 

l-l,.2) 

[-1,.3) 

t-1,.4) 

t-1,.5) 

[-l,.6) 

t-l,-.8) 


.  Cumulative  Frequency  Distribution  (Percent) 


Alg.  I 


6Ty  Alg.  II 


Alg.  Ill 


2.11 


10.53 

13.68 

17.89 

26.32 
32.63 
36.8J^ 

46.32 

55.79 

81.05 


12.63 

15.79 


20.00 


30.53 

45.26 

58.95 

87.37 


5.26 


14.74 

25.26 

42.11 

67.37 

88.42 


t-1,  .1) 


100.00 


100.00 


100.00 


TABUS  F.7 


Ctimulative  Frequency  Distribution  of  the  Daily  Correlation 
Coefficients  of  the  Vertical  Time.  Delay  Residuals  over 
Evalxoation  Condit-*on  5 


Correlation  Pair:  Stanford  ATSl  -  Stanford  ATS3 
Total  Number  of  Correlation  Coefficients:  29 


Correlation  Coefficient 

.  Cumulative  Frequency  Distribution  (Percent) 

J 

Interval 

6t^  Alg.  I 

Alg.  II 

6Ty  Alg.  Ill 

A 

% 

0 

•  0 

0 

CO 

• 

1 

1 

0 

Q 

0 

i 

1 

0 

0 

0 

% 

1 

• 

1 

H 

1 

0 

■  0 

0 

I 

1 

1 

[-l,-.5) 

0 

0 

C 

1 

1 

0 

0 

0 

1 

0 

0 

0 

'4 

1 

3.45 

0 

0 

3.45 

0 

3.45 

4 

1-1,0) 

3.45 

3.4p 

6.90 

'1 

1 

l-l,.l) 

3.45 

3.45  ” 

6.90 

f 

'1 

.  « 

-S 

6,90 

3.45 

6.90 

t 

H 

• 

U> 

6,90 

3.45 

6.90 

4 

-J 

l-l,.4) 

6.90 

3.45 

6.90 

i 

6.90 

3.45 

6.90 

6.90 

10.34 

10.34 

it 

l-l,.7) 

10.34 

17.24 

10.34 

"i 

[-1,.8) 

20.69 

31.03 

17.24 

1 

[-1,-9)  • 

31.03 

68.97 

48.28 

'S 

4 

l-l,  1] 

100.00 

100.00 

100.00 

% 

1 

wwtiar^wi^ 


.t»CaK8»««^W»»J 


TABIB  F.8 


Cumulative  Frequency  Distribution  of  the  Daily  Correlation 
Coefficients  of  the  Vertical  Tinx2  Delay  Residuals  over 
Evaluation  Condition  5 


Correlation  Fiiri  Rosman  ATS3  -  Urbana  ATS3 
Total  Nvimber  of  Correlation  Coefficients:  85 


Correlation  Coefficient 
Interval 


.  Cumulative  Frequency  Distribution  (Percent) 


6Ty  Alg.  I 


6Ty  Alg.  II 
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6t^  Alg.  Ill 


0 

■  0 

0 

1 

H 

1 

• 

OD 

0 

0 

0 

0 

0 

0 

0 

« 

0 

0 

0 

.  0 

0 

[-1,-4) 

0 

0 

0 

[-l,-.3) 

0 

0 

0 

l-l,-.2) 

0 

0 

0 

l-l,-.l) 

0 

0 

0 

[-1,0) 

0 

0 

0 

[-l,.l) 

0 

0 

0 

[-l,.2) 

0 

0  .  . 

0 

[-1,.3) 

p 

0 

0 

[-1,.4) 

1.18 

0 

1.18 

[-1,-5) 

2.53 

0 

1.18 

[-1,.6) 

5.88 

3.53 

1.18 

[-1,-7) 

11.76 

8.24 

4.71 

[-1,.8) 

23.53 

18.82 

10.59 

[-1,.9)  ■ 

''3.53 

!:0.00 

?6.)-7 

[-1,  11 

100.00 

100.00 

100.00 

I  M 


J  1 


rS 


L 


4 


•4 

1 


1  a 

I  S 

i  I 

j  2 


ii 


TABUS  P.9 

Cumulative  I’requency  Distribution  of  the  Daily  Correlation 
Coefficients’  of  the  Vertical  Time  Delay  Residuals  over 
Evaluation  Condition  5 


Corre nation  Pair:  Arecibo  ATS3  -  Cold  Bay  ATSl 
Total .Nvunber  of  Correlation  Coefficients:  30 


Correlation  Coefficient 
Interval 

•  Cumulative  Frequency  Distribution  (Percent) 

6t^  Alg.  1  6Ty  Alg.  II  6t^  Alg.  Ill 

% 

Oi 

0 

0 

i 

1 

F 

0 

0 

0 

[-1,-7) 

0 

0 

0 

i 

fi 

[-l,-.6) 

6.67 

3.33 

6.67 

1 

3 

s 

1 

M 

[-l,-.5) 

16.67 

10.00 

16.67 

20.00 

13.33 

23.33 

1 

3 

1 

l-l,-.3) 

23.33 

30.00 

36.67 

? 

\ 

1 

26.67 

36.67 

50.00 

[-1,-.!) 

36.67 

46.67 

60.00 

; 

I 

i 

1 

t-1,0) 

40.00 

66.67 

70.00 

50.00 

73.33 

76.67 

56.67 

83.33 

83.33 

.1 

[-1,  *3) 

66.67 

93.33 

90.0c 

l-l,.4) 

73.33 

96.67 

93.33 

[-l,.5) 

80.00 

100.00 

96.67 

i 

i 

[-1,.6) 

93.33 

100.00 

100.00 

t 

l-l,.?) 

96.67 

100.00 

100.00 

i 

i 

[-l,-.8) 

100.00 

100.00 

100.00 

1 

[-l,.9) 

300.00 

3,00.00 

300,00 

•  1 

% 

4 

[-1,  1] 

100.00 

100; 00 

100.00 

1 

4 

i 
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APPENDIX  G 

BASIC  lONOSHBRIC  FOHMUIAE 


Under  free  space  (vacuum)  conditions  the  satellite 
ranging  signal  would  propagate  at  the  speed  of  light,  c  , 
so  that  the  time  of  propagation  of  the  signal  would  he  a 
direct  measure  of  the  satellite -to-navigator  geometric  range 
as 


T(t) 


(G.l) 


where 

p(t) 


satellite-navigator  geometric  range  at 
time  t 


T(t)  =  time  of  propagation  of  the  ranging  signal 
between  satellite  and  navigator  at  time  t 


However,  in  the  presence  of  the  ionosphere,  the  interaction 
of  the  ranging  signal  with  the  medium  resulto  in  variations 
of  Iho  sa^nal  velocity  along  the  propagation  path.  In  this 
case 


T(t) 


where 


P(t) 


v(?,t) 


ds 


J  -^1--  (G.2) 

P(t)  v(r,t) 


phase  (optical)  path  between  satellite 
and  nav'gator  at  time  t 

position  vector  to  a  point  on  ?(t) 
from  a  geocentric  coordinate  system 

magnitude  of  the  velocity  of  the  ranging 
signal  at  ?  and  t 

differential  element  of  arc  along  P(t) 


In  oraer  to  determine  expressions  for  the  ionospheric 
effects  on  the  ranging  signals,  it  is  necessary  to  specify 
v(?,t)  in  the  ionosphere.  There  are  three  velocities  of 
propagation  associated  with  the  propagation  of  electro¬ 
magnetic  energy:  jhase,  group,  and  signal.  The  velocity 
of  interest  in  any  particular  system  is  a  function  of  the 
natures  of  both  the  measurement  and  the  mediiun.  In  the 
present  study,  the  velocity  of  interest  is,  strictly  speaking. 
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the  signal  velocity  [Brillouin,  i960].  However,  at  I6OO  MHz, 
the  absorption  and  dispersion  of  the  signal  should  be  small 
so  that  the  group  velocity  would  be  a  good  approximation  to 
the  signal  velocity.  In  the  present  study  we  will  ignore 
absorptive  effects  and  assume  that  the  ranging  signal  prop¬ 
agates  at  the  group  velocity.  The  error  involved  in  this 
assumption  is  small  but  is  an  open  question  that  is  in  need 
of  further  study. 

The  group  velocity  of  the  ranging  signal  witl; 
carrier  frequency  f  is  derived  from  the  phase  refractive 
index,  n  ,  of  the  ionosphere  as 


v(?,t)  =  4.  f 

ki(r,t) 


(G.3) 


where  n(r,t)  is  given  by  the  Appleton-Hartree  equation 
[Kelso,  196k],  Ignoring  absorptive  effects 


^o,e(?,fc)  =  1  - 


fN(r,t)  ^ 


where  in  rationalized  ml's  units 

H  (r,t)  s  incLex  of  refraction  at  r  and  time  t  for 
either  the  ordinary  or  extraordinary  wave 
for  the  two  senses  of  circular  polarization 


s  carrier  frequency,  Hz 


Un^me 


N(r,t)  =  electron  density  at  r  and  t 


s  electron  charge 


H  electron  mass 


s  permittivity  j''  free  space 
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n  a.', 


(yf)^ 

2(i-(yf)^) 


•  (yf)^ 

Mi-(yf)‘^)2 


+ 


{Tjtf 


1/2 


=  fjj(r,t)  cos  0 
sin  0 

fjj(r,t)  =  '2n~^~  H(r,t)  ,  electron  gyro-frequency 

/“♦  *  -♦ 

H(r,t)  s  magnetic  field,  strength  at  r  and  t 

0(r,t)  =  angle  between  propagation  vector  and 

magnetic  field  vector  at  ?  and  t 

=  permeability  of  free  space 


At  ionospheric  heights,  the  electron  gyro-frequency, 
fjj  ,  is  of  the  order  of  1,4  MIz,  thx:is  for  carrier  frequencies 
at  L-band  the  effect  of  the  magnetic  field  on  the  signal 
would  be  small.  Neglecting  the  magnetic  field  effects, 

M(r,t)  simplifies  to 


where 


|j(r,t) 


(G.5) 


K  =80,6  in  rationalized  inks  units. 


Then  from  Eqs,  (G,3)  and  (G,5)  the  group  velocity  of  the 
ranging  signal  is  given  approximately  as 


v(r,t)  «  c  1  -  N(r,t) 


1/2 


(G,6) 


In 

are  of  the  order 
of  the  order  of  I6OO  MHz 


the  ionosAere,  the  maximum  values  of  N(r,t) 
•der  of  10^  electrons/m3.  Thus  for  frequencies 


f^  »  K  N(r,t) 


-273- 


